###########################################################
#
#  Duffing oscilator
#
#  parametry a,f,w > 0 
#
#
# 
#
###########################################################
#
#       parameter values defining model
#
par a=0.5,f=0,w=0.2
#
#      
#      diferencialni rovnice
#
x'=v
v'=-a*v+x-x^3+f*cos(z)
z'=2*Pi*w
#
#      počáteční podmínky
#
x(0)=0.1
v(0)=0.1
z(0)=0
#
#
# Numerické parametry:
#
#
@ total=300,xlo=-2,xhi=2,ylo=-2,yhi=2,dt=0.01,maxstor=20000,xp=x,yp=v, bounds=10000, maxstor=10000000
#
done
