步驟：
1，將一個大正方形作為最基本的圖形；
2，將第一步得到的陣列縮小到三分之一；
3，將第二部得到的陣列複製八份，平移到相應的位置。
figure%畫圖
hold on
for k=1:8^n
fill(a(4*k-3) a(4*k-2) a(4*k-1) a(4*k)], [b(4*k-3) b(4*k-2) b(4*k-1) b(4*k),'b']
end
hold off
axis off%不要座標軸
axis equa%各座標軸同比例
set(findobj(gcf,'type','patch'),'edgecolor','none')%將線條隱去
步驟：
1，確定這個大正方形的四個角的座標，用兩個陣列分別表示橫座標和縱座標。可以把最左下角的座標設定為原點。
2，把這個大正方形填充成藍色。
3，將第一步得到的陣列縮小到原來的三分之一；
4，將第三步得到的陣列平移到原來大正方形的正中間位置。
function ditan(n)
a=[0 1 1 0];
b=[0 0 1 1];
fill(a,b,'b')%將正方形填充藍色
for k=1:n%迭代n次
for m=1:3^(k-1)%按行
for p=1:3^(k-1)%按列
fill(1/3^k+1/3^(k-1)*(p-1)+a/3^k,1/3^k+1/3^k+1/3^(k-1)*(m-1)+b/3^k,'w')
end
end
end
hold off
axis off%不要座標軸
axis equal%各座標軸同比例
set(findobj(gef,'type','patch'),'edgecolor','none')%將線條隱去