function dy = fduffeq(t,y)
dy = zeros(2,1);    
delta=0.25;
gamma=0.30;
om=1.0;
dy(1) = y(2);
dy(2) = y(1)-y(1)^3-delta*y(2)+gamma*cos(om*t);
