Question: suggestion on faster solving set of nonlinear ODEs

#first_question :how can i solve set of nonlinear ODEs,faster or using any packages ?
#second_question :what can be some boundary conditions for this type of nonlinear ODEs? how many BCs are required for this set of nonlinear ODEs? ( to use numeric solution)

 


 

restart:with(DEtools):with(DifferentialAlgebra):

eq[1]:=diff(N(r),r$2)+2/r*diff(N(r),r)+diff(phi(r),r)/phi(r)*diff(N(r),r)-mu^2/(32*phi(r))*N(r);

diff(diff(N(r), r), r)+2*(diff(N(r), r))/r+(diff(phi(r), r))*(diff(N(r), r))/phi(r)-(1/32)*mu^2*N(r)/phi(r)

(1)

eq[2]:=diff(phi(r),r$2)+2/r*diff(phi(r),r)-1/2*diff(phi(r),r)^2/phi(r)-8*diff(N(r),r)^2/(omega*(1-2*G*M/r))*phi(r);

diff(diff(phi(r), r), r)+2*(diff(phi(r), r))/r-(1/2)*(diff(phi(r), r))^2/phi(r)-8*(diff(N(r), r))^2*phi(r)/(omega*(1-2*G*M/r))

(2)

dsolve({eq[1],eq[2]});

``


 

Download nonlinear_ODE.mw

Please Wait...