N=100000;
A=[-0.1 0.35 ; 0.2 0.5];
a=[0.5 ; 0.4];
g=zeros(2,2,6);
g(:,:,1)=[1 0 ; 0 1];
g(:,:,2)=[1/2 -sqrt(3)/2 ; sqrt(3)/2 1/2];
g(:,:,3)=g(:,:,2)*g(:,:,2);
g(:,:,4)=[1 0 ; 0 -1];
g(:,:,5)=g(:,:,4)*g(:,:,2);
g(:,:,6)=g(:,:,4)*g(:,:,3);

//x=zeros(2,1);
X=zeros(2,N);

for i=1:N-1
    X(:,i+1)=g(:,:,floor(1+6*rand(1,1,'def')))*(A*X(:,i)+a);
//    plot2d(x(1),x(2),0,rect=[-1.5,-1.5,1,1.5])
end

plot2d(X(1,:),X(2,:),0,rect=[-1.5,-1.5,1,1.5]);