n = 33 (m = 16): trajectories parallel to the side X_0 X_(n-1); exact arithmetic in Q(zeta_66)
j= 1  simple preclosed piece: 2 segments, length 3.9639 = lambda sin(2j pi/n)/sin(pi/n) (exact); x_2 = r^b x_0 with b = 29 = -4j mod n; factor 33 = n/gcd(n,j); closed trajectory: 66 segments   [t = 3/10, 1/7, 5/9; midpoint t = 1/2: 1 segment, half the length, same factor]
j= 2  simple preclosed piece: 2 segments, length 7.7845 = lambda sin(2j pi/n)/sin(pi/n) (exact); x_2 = r^b x_0 with b = 25 = -4j mod n; factor 33 = n/gcd(n,j); closed trajectory: 66 segments   [t = 3/10, 1/7, 5/9; midpoint t = 1/2: 1 segment, half the length, same factor]
j= 3  simple preclosed piece: 2 segments, length 11.3237 = lambda sin(2j pi/n)/sin(pi/n) (exact); x_2 = r^b x_0 with b = 21 = -4j mod n; factor 11 = n/gcd(n,j); closed trajectory: 22 segments   [t = 3/10, 1/7, 5/9; midpoint t = 1/2: 1 segment, half the length, same factor]
j= 4  simple preclosed piece: 2 segments, length 14.4537 = lambda sin(2j pi/n)/sin(pi/n) (exact); x_2 = r^b x_0 with b = 17 = -4j mod n; factor 33 = n/gcd(n,j); closed trajectory: 66 segments   [t = 3/10, 1/7, 5/9; midpoint t = 1/2: 1 segment, half the length, same factor]
j= 5  simple preclosed piece: 2 segments, length 17.0613 = lambda sin(2j pi/n)/sin(pi/n) (exact); x_2 = r^b x_0 with b = 13 = -4j mod n; factor 33 = n/gcd(n,j); closed trajectory: 66 segments   [t = 3/10, 1/7, 5/9; midpoint t = 1/2: 1 segment, half the length, same factor]
j= 6  simple preclosed piece: 2 segments, length 19.0522 = lambda sin(2j pi/n)/sin(pi/n) (exact); x_2 = r^b x_0 with b =  9 = -4j mod n; factor 11 = n/gcd(n,j); closed trajectory: 22 segments   [t = 3/10, 1/7, 5/9; midpoint t = 1/2: 1 segment, half the length, same factor]
j= 7  simple preclosed piece: 2 segments, length 20.3545 = lambda sin(2j pi/n)/sin(pi/n) (exact); x_2 = r^b x_0 with b =  5 = -4j mod n; factor 33 = n/gcd(n,j); closed trajectory: 66 segments   [t = 3/10, 1/7, 5/9; midpoint t = 1/2: 1 segment, half the length, same factor]
j= 8  simple preclosed piece: 2 segments, length 20.9212 = lambda sin(2j pi/n)/sin(pi/n) (exact); x_2 = r^b x_0 with b =  1 = -4j mod n; factor 33 = n/gcd(n,j); closed trajectory: 66 segments   [t = 3/10, 1/7, 5/9; midpoint t = 1/2: 1 segment, half the length, same factor]
j= 9  simple preclosed piece: 2 segments, length 20.7318 = lambda sin(2j pi/n)/sin(pi/n) (exact); x_2 = r^b x_0 with b = 30 = -4j mod n; factor 11 = n/gcd(n,j); closed trajectory: 22 segments   [t = 3/10, 1/7, 5/9; midpoint t = 1/2: 1 segment, half the length, same factor]
j=10  simple preclosed piece: 2 segments, length 19.7930 = lambda sin(2j pi/n)/sin(pi/n) (exact); x_2 = r^b x_0 with b = 26 = -4j mod n; factor 33 = n/gcd(n,j); closed trajectory: 66 segments   [t = 3/10, 1/7, 5/9; midpoint t = 1/2: 1 segment, half the length, same factor]
j=11  simple preclosed piece: 2 segments, length 18.1389 = lambda sin(2j pi/n)/sin(pi/n) (exact); x_2 = r^b x_0 with b = 22 = -4j mod n; factor  3 = n/gcd(n,j); closed trajectory: 6 segments   [t = 3/10, 1/7, 5/9; midpoint t = 1/2: 1 segment, half the length, same factor]
j=12  simple preclosed piece: 2 segments, length 15.8291 = lambda sin(2j pi/n)/sin(pi/n) (exact); x_2 = r^b x_0 with b = 18 = -4j mod n; factor 11 = n/gcd(n,j); closed trajectory: 22 segments   [t = 3/10, 1/7, 5/9; midpoint t = 1/2: 1 segment, half the length, same factor]
j=13  simple preclosed piece: 2 segments, length 12.9473 = lambda sin(2j pi/n)/sin(pi/n) (exact); x_2 = r^b x_0 with b = 14 = -4j mod n; factor 33 = n/gcd(n,j); closed trajectory: 66 segments   [t = 3/10, 1/7, 5/9; midpoint t = 1/2: 1 segment, half the length, same factor]
j=14  simple preclosed piece: 2 segments, length 9.5975 = lambda sin(2j pi/n)/sin(pi/n) (exact); x_2 = r^b x_0 with b = 10 = -4j mod n; factor 33 = n/gcd(n,j); closed trajectory: 66 segments   [t = 3/10, 1/7, 5/9; midpoint t = 1/2: 1 segment, half the length, same factor]
j=15  simple preclosed piece: 2 segments, length 5.9009 = lambda sin(2j pi/n)/sin(pi/n) (exact); x_2 = r^b x_0 with b =  6 = -4j mod n; factor 11 = n/gcd(n,j); closed trajectory: 22 segments   [t = 3/10, 1/7, 5/9; midpoint t = 1/2: 1 segment, half the length, same factor]
j=16  simple preclosed piece: 2 segments, length 1.9909 = lambda sin(2j pi/n)/sin(pi/n) (exact); x_2 = r^b x_0 with b =  2 = -4j mod n; factor 33 = n/gcd(n,j); closed trajectory: 66 segments   [t = 3/10, 1/7, 5/9; midpoint t = 1/2: 1 segment, half the length, same factor]
factors d_1..d_m: [33, 33, 11, 33, 33, 11, 33, 33, 11, 33, 3, 11, 33, 33, 11, 33] ; gcd = 1 ; a quotient equals n: True
DONE [1.6s]
