% script principal

clear all
clf
rand('state',sum(100*clock));
[M,p]=initIFS;
P=[0,0];
for cpt=1:10000
    P=imageIFS(P,p,M);
    plot(P(1),P(2),'k.');
    hold on
end
axis([-7,7,0,10])
axis equal

function j=choix(p)
x=rand;
j=1;
s=p(1);
while s<x
    j=j+1;
    s=s+p(j);
end

function [M,p]=initIFS()
M=[0,0,0,0.16,0,0;...
   0.8492,0.0371,-0.0371,0.8492,0,1.6;...
   0.1968,-0.2264,0.2264,0.1968,0,1.6;...
   -0.15,0.2834,0.2598,0.2378,0,0.44];
p=[0.1,0.7,0.1,0.1];

function Q=imageIFS(P,p,M)
j=choix(p);
Q=[M(j,1)*P(1)+M(j,2)*P(2)+M(j,5),...
   M(j,3)*P(1)+M(j,4)*P(2)+M(j,6)];