function a = A363348( numIter )
    a = [3,-2, 3,-2, 3, 2, 0, 2, -3, 2, 3, 2,-3, 2];
    L = [14, 140, 1588];
    r1 = 6; r2 = 6; r3 = 6; r4 = 8; r5 = 2; r6 = 2; r7 = 4;
    r8 = 1; r9 = 1; r10 = 1; r11 = 1; r12 = 1; r13 = 1; r14 = 1;
    c1 = r2; c2 = c1; c3 = r2 + r3; c4 = r2 + r4; c5 = r2 + r7; c6 = c4;
    for k = 1:numIter 
        if k > 3
            L(k) = 12*L(k-1) - 7*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);
            r8(k) = r12(k-1) + r8(k-1) + r14(k-1);
            r9(k) = r9(k-1) + r14(k-1);
            r10(k) = 2*r13(k-1) + r10(k-1) + r11(k-1) + r14(k-1);
            r11(k) = 2*r13(k-1) + 3*r10(k-1) + r11(k-1);
            r12(k) = r13(k-1) + r10(k-1);
            r13(k) = r12(k-1) + r13(k-1) + r14(k-1);
            r14(k) = r13(k-1) + 2*r10(k-1);
            s1 = 0; s2 = 0; s3 = 0; s4 = 0; s5 = 0; s6 = 0;
            for m = 1:k-1
                s1 = s1 + r9(k+1-m)*L(m);
                s2 = s2 + L(m);
                s3 = s3 + (r9(k+1-m) + r10(k+1-m) - 1)*L(m);
                s4 = s4 + (r9(k+1-m) + r11(k+1-m) - 1)*L(m);
                s5 = s5 + (r9(k+1-m) + r14(k+1-m) - 1)*L(m);
            end
            c1(k) = r2(k) + s1;
            c2(k) = c1(k) - s2;
            c3(k) = r2(k) + r3(k) + s3;
            c4(k) = r2(k) + r4(k) + s4;
            c5(k) = r2(k) + r7(k) + s5 - s2;
            c6(k) = c4(k) - s2;
        end
        a = [a a(1:c1(k)-1) -a(c1(k)) a(c2(k)+1:L(k)) a(1:c3(k)-1) -a(c3(k)), a(c2(k)+1:L(k)) a(1:c3(k)-1) -a(c3(k)) a(c2(k)+1:L(k)) a(1:c4(k)-1) -a(c4(k)) a(c5(k)+1:L(k)) a(1:c3(k)-1) -a(c3(k)) a(c2(k)+1:L(k)) a(1:c3(k)-1) -a(c3(k)), a(c2(k)+1:L(k)) a(1:c4(k)-1) -a(c4(k)) a(c6(k)+1:L(k))]; 
    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