[Maxima] Phase space trajectories
j.cubizolles at free.fr
Wed Nov 25 17:52:55 CST 2009
Hi, I'm trying to display several trajectories in phase space
corresponding to the same ODEs but with different initial conditions.
I'd also like to use sliders.
I tried to use the "parameters" options in plotdf but it seems you can
only use it for the ODE expressions, not the initial conditions.
So I turned to "manually" do it. I'm using a list of rk to numerically
solve the ODE system and with_slider to draw the trajectories. The code
is like :
sol : makelist(rk([-x,v],[v,x],[0,i],[t,0,10,.5]),i,1,3)$
but I'm running into the error :
Subscript must be an integer:
-- an error. To debug this try debugmode(true);
I guess the makelist(i,i,1,2) does not create a list of integers as i
thought it would. How could I fix it ?
Btw if anyone has a simpler way to draw such phase-space trajectories,
I'd be happy to try it.
Thanks for your help.
More information about the Maxima