Helge von Koch 雪花（MATLAB）

%Helge von Koch
clear all
clc
new=[0.5+(sqrt(3)/2)*i,-0.5+(sqrt(3)/2)*i,...
    0,0.5+(sqrt(3)/2)*i];
plot(new);
pause(0.5);
for k=1:4
    old=new;
    [m,n]=size(old);
    n=n-1;
    for j=0:n-1
        diff=(old(j+2)-old(j+1))/3;
        new(4*j+1)=old(j+1);
        new(4*j+2)=old(j+1)+diff;
        new(4*j+3)=new(4*j+2)+diff*((1-sqrt(3)*i)/2);
        new(4*j+4)=old(j+1)+2*diff;
    end
    new(4*n+1)=old(n+1);
    plot(new);
    pause(0.5);
end