x = 0.2;
y = 0.1;
plot(x,y);
hold on;
omega = sqrt(2);
mu = -0.1;
for i = 1:3000
   a = [cos(omega) -1*sin(omega); sin(omega) cos(omega)]*[(1+mu)*x-x*(x^2+y^2); (1+mu)*y-y*(x^2+y^2)];
   x = a(1);
   y = a(2);
   plot(x,y);
end