function a = A363445( numIter )
    a = [3,-2, 3,-2, 3, 2, 0, 2, -3, 2, 3, 2,-3, 2];
    L = [0,14,56,202];
    r1 = 6; r2 = 6; r3 = 6; r4 = 8; r5 = 2; r6 = 2; r7 = 4;
    c1 = r2; c2 = r2 + r3; c3 = r2 + r4; c4 = r2 + r7; 
    for k = 1:numIter 
        if k > 3
            L(k+1) = 6*L(k) - 10*L(k-1) + 6*L(k-2) - L(k-3);
        end
        if k > 1
            r1(k) = r5(k-1) + r1(k-1) + r7(k-1);
            r2(k) = r2(k-1) + r7(k-1);
            r3(k) = 2*r6(k-1) + r3(k-1) + r4(k-1) + r7(k-1);
            r4(k) = 2*r6(k-1) + 3*r3(k-1) + r4(k-1);
            r5(k) = r6(k-1) + r3(k-1);
            r6(k) = r5(k-1) + r6(k-1) + r7(k-1);
            r7(k) = r6(k-1) + 2*r3(k-1);
            c1(k) = r2(k) + L(k);
            c2(k) = r2(k) + r3(k) + L(k);
            c3(k) = r2(k) + r4(k) + L(k);
            c4(k) = r2(k) + r7(k) + L(k);
        end
        a = [a a((1+L(k)):c1(k)-1) -a(c1(k)) a(c1(k)+1:c2(k)-1) -a(c1(k)), a(c1(k)+1:c2(k)-1) -a(c1(k)) a(c1(k)+1:c3(k)-1) -a(c4(k)) a(c4(k)+1:c2(k)-1) -a(c1(k)) a(c1(k)+1:c2(k)-1) -a(c1(k)), a(c1(k)+1:c3(k)-1) -a(c3(k)) a(c3(k)+1:L(k+1))];
    end
    drawSeq( a ); axis equal;
end
function drawSeq( s )
    phi = 0; x = 0; y = 0;
    l = sqrt(3);
    x = [x x(end)+l*cos(phi)];
    y = [y y(end)+l*sin(phi)];
    for k = 1:length(s)
        if abs(s(k)) == 3
            if l == sqrt(3)
                l = 1;
            else
                l = sqrt(3);
            end
        end
        phi = phi + (s(k)/6)*pi;
        x = [x x(end)+l*cos(phi)];
        y = [y y(end)+l*sin(phi)];
    end
    plot(x,y);
end