Advanced Search

Question: Numerical solution of ODE in polar coords with singularities

Question:Numerical solution of ODE in polar coords with singularities

TigraDance 0 Maple
plot student
July 01 2010
0

Please help plot solution of ODE in polar coords or tell how can i define a piecewise function that will make me able to plot the right and left side of singularities.

>sys1:=diff(r(t),t$2)=2*cos(teta(t))/r(t)^3,diff(teta(t),t$2)=sin(teta(t))/r(t)^3;funcs:=[r(t),teta(t)]:InC:=r(0)=1,teta(0)=15,D(r)(0)=0,D(teta)(0)=-2;

>FNUM:=dsolve({sys1,InC},numeric,funcs);

>odeplot(FNUM,funcs,t=0..230,numpoints=2000,coords=polar,view=[-20..20,-20..20]);

Thanks
Please Wait...