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