¿Cómo dibujar la pirámide de Sierpinski usando MATLAB?

n=8;?% número de capas de iteración, la capa 1 tiene solo 1 triángulo, la capa 2 tiene 4 triángulos

x=cosd([90?210?330]') ;

y=sind([90?210?330]');

tri=[1?2?3];

tnum=1; dnum=3;

for?ii=2:n

tri2=zeros(3*tnum,3);

nn=1;

for?jj=1:tnum

xnew=(x(tri(jj,[1?2?3]))+x(tri(jj,[2?3?1] )))/2;

ynew=(y(tri(jj,[1?2?3]))+y(tri(jj,[2?3?1])))/2 ;

x=[x;xnew];

y=[y;ynew];

tri2(nn,:)=[tri(jj, 1),dnum+1,dnum+3];

tri2(nn+1,:)=[dnum+1,tri(jj,2),dnum+2];

tri2(nn+2,:)=[dnum+3,dnum+2,tri(jj,3)];

nn=nn+3;

dnum= dnum+3;

fin

tri=tri2;

tnum=3*tnum;

fin

triplot(tri,x,y);

eje?imagen;

eje?apagado;