%Plots the Sierpinski fractal
x=[0 0 1];
y=[1 0 0];
for i=1:7    x=[x/2 x/2 x/2 + 1/2];
    y=[y/2 + 1/2 y/2 y/2];
    clf; plot(x, y); axis equal; 
    pause(1/5);
end;

%M=[1 1/2; 0 sqrt(3)/2]*[x; y];
%plot(M(1,:), M(2,:)); axis equal;
