# the undamped pendulum
theta'=v
v'=-(g/l)*sin(theta)
par g=9.8,l=2,m=1
aux E=m*((l*v)^2/2+g*l*(1-cos(theta)))
done