Question: how do i solve the following system of differential equation

dsys := {1/2*diff(theta(t),`$`(t,2))*cos(phi(t))+5/4*diff(theta(t),`$`(t,2))+5/4*diff(phi(t),`$`(t,2))+1/2*sin(phi(t))*diff(theta(t),t)^2, -1/2*mu(t)*sin(phi(t))*diff(theta(t),`$`(t,2))-diff(phi(t),`$`(t,2))+1/2*mu(t)*cos(phi(t))*diff(theta(t),t)^2+5/4*diff(mu(t),`$`(t,2))-150000000000/(1+exp(27.0-15*sin(theta(t))-15*sin(theta(t)+phi(t))))^2*cos(theta(t)+phi(t))*exp(27.0-15*sin(theta(t))-15*sin(theta(t)+phi(t))), (-1-mu(t)*sin(phi(t)))*diff(theta(t),`$`(t,2))-1/2*mu(t)*sin(phi(t))*diff(phi(t),`$`(t,2))-1/2*mu(t)*cos(phi(t))*diff(phi(t),t)^2-mu(t)*cos(phi(t))*diff(phi(t),t)*diff(theta(t),t)-diff(mu(t),t)*sin(phi(t))*diff(theta(t),t)-diff(mu(t),t)*sin(phi(t))*diff(phi(t),t)+diff(mu(t),`$`(t,2))*(1/2*cos(phi(t))+5/4)+10000000000/(1+exp(27.0-15*sin(theta(t))-15*sin(theta(t)+phi(t))))^2*(-15*cos(theta(t))-15*cos(theta(t)+phi(t)))*exp(27.0-15*sin(theta(t))-15*sin(theta(t)+phi(t)))};

boundary conditions

bcs:={theta(0)=Pi/6,phi(0)=-Pi/6,phi(30)=Pi/6,theta(30)=5*Pi/6,D(phi)(0)=0,D(theta)(30)=0};

i am trying to solve the above system of differential equation with the given B.C. for theta, phi and mu, all are function of time (t), with mu as a parameter. 

i tried using continuation as well as using a higher value of initmesh and maxmesh, but "initial newton iteration not converging" always creeps in.

i would be thankful if anyone could help me in solving the above system.

thanks in advance

anshul

Please Wait...