NP=1000; % antal punkter
clear x; x(2,NP)=0; % allokera resultatvektorn
x(:,1)=[0;0]; % startpunkt
NS=10; % skippa lite punkter i början

A1=[0 0;0 0.16]; T1=[0;0];  p1=0.01;
A2=[0.84 0.04;-0.04 0.85]; 	T2=[0;1.6];
p2=0.85; p12=p1+p2;
A3=[0.2 -0.26;0.23 0.22]; 	T3=[0;1.6];
p3=0.07; p123=p12+p3;
A4=[-0.15 0.28;0.26 0.24]; 	T4=[0;0.44];
p4=0.07;
for j=1:NP-1,
 r=rand;
 if(r<=p1)
  x(:,j+1)=A1*x(:,j)+T1;
 elseif((p1 < r) & (r<=p12))
  x(:,j+1)=A2*x(:,j)+T2;
 elseif((p12 < r) & (r<=p123))
  x(:,j+1)=A3*x(:,j)+T3;
 else
  x(:,j+1)=A4*x(:,j)+T4;
end
end
plot(x(1,NS:NP),x(2,NS:NP),'g.')
axis('equal')