N=10000;
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);
for i=1:N
    x=g(:,:,floor(1+6*rand(1,1,'def')))*(A*x+a);
    plot2d(x(1),x(2),0,rect=[-1.5,-1.5,1,1.5])
end