# forced duffing equation
x'=v
v'=a*x*(1-x^2) -f*v+c*cos(omega*t)
par a=1,f=.2,c=.3,omega=1
init x=.1,v=.1
done