# double well potential
x'=y
y'=(x-x^3)/m
par m=1
aux e=m*y^2/2+x^4/4-x^2/2
@ xp=x,yp=y,xlo=-1.5,ylo=-1.5,xhi=1.5,yhi=1.5
done