axis([-6 6 -6 6]);
mm=100;
n=50;
w=zeros(mm,2*n);
z=zeros(mm,2*n);
for i=1:n,
   t0 = 0;
   tfinal =20;
   y0 = [1+3*cos(0.5+2*i*pi/n) 3*sin(0.5+2*i*pi/n)];
   tol = 1.e-7;		% Accuracy
   trace = 0;
   [t,y] = myode45('peetwo',t0,tfinal,y0,tol,trace);
   [m1 m2]=size(y);
   w(:,i)=[y(:,1)' zeros(mm-m1,1)']';
   z(:,i)=[y(:,2)' zeros(mm-m1,1)']';
end
for i=1:n,
   t0 = 0;
   tfinal =20;
   y0 = [0.6*cos(0.5+2*i*pi/n) 0.6*sin(0.5+2*i*pi/n)];
   tol = 1.e-5;		% Accuracy
   trace = 0;
   [t,y] = myode45('peetwo',t0,tfinal,y0,tol,trace);
   [m1 m2]=size(y);
   w(:,n+i)=[y(:,1)' zeros(mm-m1,1)']';
   z(:,n+i)=[y(:,2)' zeros(mm-m1,1)']';
end
plot(w,z,'.');