function ret=koch(p0,p4,n)
    if n==0
        ret=[p0,p4];
    else
        p1=p0+1/3*(p4-p0);
        p2=p1+(p1-p0)*exp(i*pi/3);
        p3=p0+2/3*(p4-p0);
        ret=[koch(p0,p1,n-1),koch(p1,p2,n-1),koch(p2,p3,n-1),koch(p3,p4,n-1)];
    end