nbpoints=10000
P=zeros(2,nbpoints);
// P est la matrice nulle à 2 lignes et nbpoints colonnes
for i = 2:nbpoints
    tirage=rand()
    if tirage <1/2 then 
        A=[[0.5,-0.5] ;[0.5,0.5]];V=[0;0];
    else  A=[[-0.5,-0.5] ;[0.5,-0.5]];V=[1;0];
 end 
 P(:,i)=A*P(:,i-1)+V;
 //P(:,i) désigne la colonne i de la matrice P
end
plot(P(1,:),P(2,:),"*b");
//les : indiquent que toutes les colonnes sont retenues
    
    