function demoKoch
global x y
a=[0;0];
b=[1;0];
e=[b(1)*cos(pi/3)-b(2)*sin(pi/3)
   b(1)*sin(pi/3)+b(2)*cos(pi/3)];
f=[b(1)*cos(-pi/3)-b(2)*sin(-pi/3)
   b(1)*sin(-pi/3)+b(2)*cos(-pi/3)];
n=4;
x=[];y=[];
koch(a,b,n);
koch(b,f,n);
koch(f,a,n);
h=fill(x,y,'c');
axis equal
axis off
set(gcf,'color','b')
set(gcf,'doubleBuffer','on')
figure(gcf)
while 1
    rotate(h,[90 0],1,[0 0 0])
    drawnow
end
