function dydt=duffing(t,y)

w0sqr=0.1; eta=0.1; F0=12;
 
dydt=[y(2) ; -w0sqr*y(1)-eta*y(1).^3+F0*cos(t)];