clear;
n=10000;
A=[-0.4 0.75 ; 0.2 -0.3 ];
a=[0 ; 0.4];
m=55;
R=[cos(2*%pi/m) -sin(2*%pi/m) ; sin(2*%pi/m) cos(2*%pi/m)];

G=zeros(2,2,m);
G(:,:,1)=diag([1,1]);

for i=2:m

     G(:,:,i)=R*G(:,:,i-1);
 
end


x=zeros(2,1);

xset("window",1);
xbasc();

for i=1:n

x=G(:,:,floor(1+m*grand(1,1,'def')))*(A*x+a);
        plot2d(x(1),x(2),0);
end,