How to tabulate $(r\cos\theta,r\sin\theta)$ from $(t,r,\theta)$ obtained from solving a system of ODE using RK4?

How to tabulate

$(r\cos\theta,r\sin\theta)$ from

$(t,r,\theta)$ obtained from solving a system of ODE using RK4?

Abby11
53 ●1 ●3 ●6

I am trying to solve a system of first-order differential equations of the following form: $\frac{dr}{dt}=f(r,\theta), \quad\frac{d\theta}{dt}=g(r,\theta)$ I am using desolve_system_rk4 which gives me a list of triples $(t,r,\theta)$ . Now I need to plot the trajectory in the xy-plane. For this, I require to obtain the x and y coordinate using $x=r\cos\theta, \quad y=r\sin\theta$ and then obtain the xy-plot. How can I tabulate the pair $(r\cos\theta,r\sin\theta)$ from $(t,r,\theta)$ ?

Comments

Please add your code using desolve_system_rk4, so that it's easier to give an answer.

rburing ( 3 years ago )

Can a parametric plot solve your problem?

tolga ( 3 years ago )

add a comment

answered 3 years ago

Emmanuel Charpentier
7804 ●8 ●52 ●155

Starting from the example given in the documentation that you could have bothered to read by typing desolve_system_rk4? :

sage: from sage.calculus.desolvers import desolve_system_rk4
sage: x,y,t=var('x y t')
sage: P=desolve_system_rk4([x*(1-y),-y*(1-x)],[x,y],ics=[0,0.5,2],ivar=t,end_points=20)

you can see that the solution P is a list of [t, x, y] triplets. You can display the x and y functions by :

sage: line2d([[i, j] for i, j, k in P])
sage: line2d([[i, k] for i, j, k in P])

(line2d interpolates between the points that list_plot displays, which is aesthetically more pleasing).Your phase trajectory can therefore be displayed by :

sage: line2d([[j,k] for i, j, k in P])

HTH,

link

add a comment

Your Answer

Please start posting anonymously - your entry will be published after you log in or create a new account.

How to tabulate $(r\cos\theta,r\sin\theta)$ from $(t,r,\theta)$ obtained from solving a system of ODE using RK4?

Comments

1 Answer

Your Answer

Question Tools

Stats

How to tabulate (rcosθ,rsinθ)(r\cos\theta,r\sin\theta) from (t,r,θ)(t,r,\theta) obtained from solving a system of ODE using RK4? savecancel

Comments

1 Answer

Your Answer

Question Tools

Stats

How to tabulate $(r\cos\theta,r\sin\theta)$ from $(t,r,\theta)$ obtained from solving a system of ODE using RK4?