tmax=100
omega=10
gamma=1
epsilon=1
delta=1
tmax=100
pts=100
#       Duffing equation
#        y dotdot + delta y dot + y + gamma y^3 = epsilon cos(omega t)
