科赫曲線 MATLAB

% 科赫曲線是一種像雪花的幾何曲線，所以又稱爲雪花曲線，它是de Rham曲線的特例。
% 科赫曲線是出現在海里格·馮·科赫的論文中，是形曲線中的一種。

x=[-1,0,1];
y=[0,sqrt(3),0];
c=complex(x,y);
c(4)=c(1);
plot(real(c), imag(c));
axis([-1.5 1.5 -1 2]);
axis square;
pause(1);
for i=1:8
    n=size(c);
    n=n(2);
    n=n-1;
    c1=complex(zeros(1,4*n),zeros(1,4*n));
    for j=1:n
        c1(4*j-3)=c(j);
        c1(4*j-2)=(2*c(j)+c(j+1))/3;
        c1(4*j)=(c(j)+2*c(j+1))/3;
        c1(4*j-1)=c1(4*j-2)+(c1(4*j)-c1(4*j-2))*complex(0.5,sqrt(3)/2);
    end
    c=c1;
    c(4*n+1)=c(1);
    plot(real(c), imag(c));
    axis([-1.5 1.5 -1 2]);
    axis square;
    pause(1);
end