function gasket (A,B,C,level);
% Rekursiv erzeugtes Sierpinski Dreieck
% input: A,B,C kartesische Koordinaten der Eckpunkte
%        level Anzahl der Rekursionschritte
%

%level
if level > 0
    MAB = (A+B)/2;
    MBC = (B+C)/2;
    MAC = (A+C)/2;
    gasket(A,MAB,MAC,level-1);
    gasket(B,MAB,MBC,level-1);
    gasket(C,MAC,MBC,level-1);
else
    fill([A(1),B(1),C(1)],[A(2),B(2),C(2)],'k'),pause(.5)
    axis([0,1,0,1])
    hold on
end
return