# cup_pend.ode
# integrate them and then color code according to velocity
# let XPP figure out the bounds
phi1'=phi1p
phi1p'=-sin(phi1)+d*(phi2-phi1)
phi2'=phi2p
phi2p'=-sin(phi2)+d*(phi1-phi2)
par d=.01
init phi1=2.5,phi1p=.4,phi2=2.9,phi2p=0
@ total=800,trans=200
@ bound=100000,maxstor=40000
@ xp=phi1,yp=phi2,xlo=20,xhi=100,ylo=20,yhi=100
done