Advanced Search

Question: plotting the phase portrait of a hamiltonian

Question:plotting the phase portrait of a hamiltonian

meyestani 0 Maple
phase
December 25 2011
0

hi

I need to plot the phase portrait of a hamiltonian of 4-dimensional phase-space, and I have an adiabatic invariant, i.e. one of the momenta is conserverd. 

I don't know how to replace a constant for the conserved momentum and plot the 2-dimensional phase space.

I am writiing 

sys := {diff(Y1(t), t) = -(diff(K, y1)), diff(Y2(t), t) = -(diff(K, y2)), diff(y1(t), t) = diff(K, Y1), diff(y2(t), t) = diff(K, Y2)}

and using the DEplot tool as 

> phaseportrait([Y1(t), y1(t)], t = -2 .. 2, [[y1(0) = 1, Y1(0) = 0]], stepsize = 0.5e-1, linecolour = sin((1/2)*t*Pi), method = classical[foreuler]);
 
but nothing happens! could anyone tell me how to do the replacement and plot the phase portrait?
Please Wait...