function dy = duffeq(t,y)
dy = zeros(2,1);    
delta=0.1;
dy(1) = y(2);
dy(2) = y(1)-y(1)^3-delta*y(2);
