A=[0.4 -0.3733 ; 0.06 0.6];
B=[-0.8 -0.1867 ; 0.1371 0.8];
a=[0.3533 ; 0];
b=[1.1 ; 0.1];
N=100000;
X=grand(2,N,'def');
p=0.5;

for i=1:N-1
         X(:,i+1)=bool2s(grand(1,1,'def')>p)*(A*X(:,i)+a-B*X(:,i)-b)+B*X(:,i)+b;
        //plot2d(x(1),x(2),0,rect=[0,0,1.4,0.9]);
end  


plot2d(X(1,:),X(2,:),0);