function tapis(x,y,l,n)

% La fonction TAPIS(X,Y,L,N) représente le tapis de Sierpinsky
% X,Y sont les coordonnées inférieures gauches du carré à
% enlever. L est la longueur du carré intérieur et N est
% le niveau ou le nombre d'itérations.

if n==0
	break
end

l=l/3;				% longueur du carré intérieur

xc=[x x+l x+l x];	% coordonnées des sommets du carré
yc=[y y y+l y+l];

fill(xc,yc,'k')		% remplissage en noir = effacement

for i=1:3
	tapis(x-2*l/3+l*(i-1),y-2*l/3,l,n-1)
end
for i=1:3
	tapis(x-2*l/3+l*(i-1),y+4*l/3,l,n-1)
end
for i=1:3
	tapis(x-2*l/3+l*(i-1)*2,y+l/3,l,n-1)
end
