n = 13 (m = 6): trajectories parallel to the side X_0 X_(n-1); exact arithmetic in Q(zeta_26)
j= 1  simple preclosed piece: 2 segments, length 3.7709 = lambda sin(2j pi/n)/sin(pi/n) (exact); x_2 = r^b x_0 with b =  9 = -4j mod n; factor 13 = n/gcd(n,j); closed trajectory: 26 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 6.6780 = lambda sin(2j pi/n)/sin(pi/n) (exact); x_2 = r^b x_0 with b =  5 = -4j mod n; factor 13 = n/gcd(n,j); closed trajectory: 26 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 8.0552 = lambda sin(2j pi/n)/sin(pi/n) (exact); x_2 = r^b x_0 with b =  1 = -4j mod n; factor 13 = n/gcd(n,j); closed trajectory: 26 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 7.5870 = lambda sin(2j pi/n)/sin(pi/n) (exact); x_2 = r^b x_0 with b = 10 = -4j mod n; factor 13 = n/gcd(n,j); closed trajectory: 26 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 5.3808 = lambda sin(2j pi/n)/sin(pi/n) (exact); x_2 = r^b x_0 with b =  6 = -4j mod n; factor 13 = n/gcd(n,j); closed trajectory: 26 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 1.9419 = lambda sin(2j pi/n)/sin(pi/n) (exact); x_2 = r^b x_0 with b =  2 = -4j mod n; factor 13 = n/gcd(n,j); closed trajectory: 26 segments   [t = 3/10, 1/7, 5/9; midpoint t = 1/2: 1 segment, half the length, same factor]
factors d_1..d_m: [13, 13, 13, 13, 13, 13] ; gcd = 13 ; a quotient equals n: False
DONE [0.1s]
