cbb556fa-0db7-499c-a47a-ecb6999bcd39 64 3 5 decimal

&[DATE] &[TIME] - &[FILENAME]

Intersection of two cylinders

X0 tmin tmax N D n f rows 1 + : f n el 0 : j 1 n range k 1 n range m j k el f j el x k el diff : for for x Jacobi m : b 1 - (x Jacobi 1 n 1 - n n submatrix * : A x Jacobi 1 n 1 - 1 n 1 - submatrix : n 2 > Δ 1 el b A 1 n 1 - 2 n 1 - submatrix augment det - : Δ n 1 - el A 1 n 1 - 1 n 2 - submatrix b augment det - : Δ n el A det - : k 2 : k n 2 - ≤ k k 1 + : Δ k el A 1 n 1 - 1 k 1 - submatrix b A 1 n 1 - k 1 + n 1 - submatrix augment det - : 11 line for 41 line Δ 1 el b : Δ n el A det : 21 line if i 1 n range D2 i el Δ i el Δ norme / : for t x D1 D2 : r X0 tmin tmax N t x D1 rkfixed : 101 line :

f x 1 el (2^x 2 el (2^+ 2 - x 2 el (2^x 3 el (2^+ 2 - 21 mat : B1 02 sqrt - 1013 -^13 mat transpose 012 - 120 D 112024 submatrix eval : γ 0.8232 - 0.4194 - 0.3827 0.5677 0.6187 - 0.543 0.009 0.6643 0.7474 33 mat : angle Ω2 angle cos 0 angle sin 010 angle sin - 0 angle cos 33 mat : 22 mat

α F B0 B1 1 col B1 2 col B1 rows 1 matrix 2 + augment (: k 1 B1 rows range L1 k el B0 k row B1 k row stack α Ω2 * γ * eval : L2 k el L1 k el 1 col L1 k el 3 col augment : S k el L2 k el 1010^1010^1010^1010^1010^16 mat transpose eval : k 1 ≡ F1 S k el eval : F1 F1 S k el stack eval : if 41 line for F1 31 line :

B11 B1 γ * 1 col B1 γ * 3 col augment eval 11 mat :

π 2 / F 3 2 / π * F B11 π F stack eval mat2sys.1 π 2 / F 3 2 / π * F B11 stack eval mat2sys.1