%%%%%%%%%%%%%% L9 ti 21.11.06 B-sali%%%%%%%%%%%%%%%
%A=[.5 -.6;.75 1.1]
phi=pi/4; 
r=1.2;  %repellori 
r=0.8   %attraktori
A=r*[cos(phi) -sin(phi);sin(phi) cos(phi)];

%x0=[1 ;1.5]
x0=[0.1; 0.1]
x0=[1;1]
close all
x=x0;
n=15;
for k=1:n
    y=A*x(:,end);
    x=[x y];
end;   

plot(x(1,:),x(2,:),'o')
hold on
plot(x(1,:),x(2,:))
plot(x(1,1),x(2,1),'rx')
