koch曲線的matlab程序

function kk_plot(n)
if nargin==0
   n=3;
end
a=0;b=10;
y=[a,b];
while n>0
h=[];m=length(y)-1;
for i=1:m
x=k(y(i),y(i+1));h=[h x];
end
h(end+1)=y(i+1);y=h;n=n-1;
end
plot(y),axis equal off
title('fractal curve')
function y=k(a,b)
c=a+(b-a)/3;
e=a+2*(b-a)/3;
d=c+(e-c)*exp(i*pi/3);
y=[a c d e];