yy1=1;
yy2=1;
plot(yy1,yy2,'g.')
hold on
alpha=1;
beta=5;
delta=0.02;
K=8;
omega=0.5;
tstep=(2*pi)/omega;
tstart=-tstep;
tfinal=0;
for jstep=1:3200
    tstart=tstart+tstep;
    tfinal=tfinal+tstep;
    y0=[yy1;yy2;1;0];
    tspan=[tstart tfinal];
    options=[];
    [t,y]=ode45(@duffn4,tspan,y0,options,alpha,beta,delta,K,omega);
    s=size(y);
    yy1=y(s(1),1);
    yy2=y(s(1),2);
    if jstep > 200
        plot(yy1,yy2,'r.')
        hold on
    end
end
hold off