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