% 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