哪位高手能教教我用matlab画koch雪花的具体步骤。先多谢拉

function plotkoch(k) %显示迭代k次后的Koch曲线图
p=[0,0;10,0]; %存放结点坐标，每行一个点，初始值为两结点的坐标
n=1; %存放线段的数量，初始值为1
A=[cos(pi/3),-sin(pi/3);sin(pi/3),cos(pi/3)]; %用于计算新的结点
for s=1:k %实现迭代过程，计算所有的结点的坐标
j=0; %
%以下根据线段两个结点的坐标，计算迭代后它们之间增加的三个
%结点的坐标，并且将这些点的坐标按次序存暂时放到r中
for i=1:n %每条边计算一次
q1=p(i,:); %目前线段的起点坐标
q2=p(i+1,:); %目前线段的终点坐标
d=(q2-q1)/3; %
j=j+1;r(j,:)=q1; %原起点存入r
j=j+1;r(j,:)=q1+d; %新1点存入r
j=j+1;r(j,:)=q1+d+d*A'; %新2点存入r
j=j+1;r(j,:)=q1+2*d; %新3点存入r
end %原终点作为下条线段的起点，在迭代下条线段时存入r
n=4*n; %全部线段迭代一次后，线段数量乘4
clear p %清空p ，注意：最后一个终点q2不在r中
p=[r;q2]; %重新装载本次迭代后的全部结点
end%
figure
plot(p(:,1),p(:,2)) %显示各结点的连线图
axis equal %各坐标轴同比例