Question: how to solve the challenging non linear coupled 1st ode when bcs are given

G := 6.6743*10^(-8);

a := 1.9501*10^24;

b := .3306;

c := 2.99792458*10^10;

d := 2.035;

pi := 3.143;

eps := 3.8220*10^35;

g(r) = 1-s(r)/0.06123;

j(r) = e^(-(1/2)*w(r))*(1-2*G*v(r)/(r*c^2))^.5

sys := diff(v(r), r) = 4*pi*r^2*eps/c^2, ics=v(0)=0

diff(u(r), r) = -G*(eps+u(r))*(v(r)+4*Pi*r^3*u(r)/c^2)/(c^2*(r^2-2*G*r*v(r)/c^2)),u(0)=1.3668*10^34

diff(w(r), r) = 1.485232054*10^(-28)*(v(r)+4.450600224*10^(-21)*pi*r^3*u(r))/(r^2-2*G*r*v(r)/c^2), conditions: w(0)=0,iterate it to find w(688240)=-2.05684, it solve must satistfy the both conditions.

diff(r^4*j(r)*(diff(g(r), r)), r)+4*r^3*g(r)*(diff(j(r), r)) = 0, conditions dg(r)/dr =0  at r=0, g(688240) =0.87214

diff(J(r), r) = (8*pi*(1/3))*(eps/c^2+u(r)/c^2)*(g(r)*j(r).(r^4))/(1-2*G*v(r)/(r*c^2)) condition J(0)=0.