e64a27d4-cfe0-46e1-86d6-2feb306d8656 193 12 6 decimal

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

Elliptic integral of the 1-st kind (LegandreKc(k)) k Kc a 1 : g 1 k 2^- sqrt : 12 mat j 125 range a1 a : a 0.5 a g + (* : g a1 g * sqrt : 13 mat eval a g - 252 -^≤ break 0 if 21 line for π 2 a * / eval 31 line :

Elliptic functions sn,cn,dn x k scdn a 1 el 1 : g 1 el 1 k 2^- sqrt : c 1 el k : 13 mat n 125 range a n 1 + el 0.5 a n el g n el + (* : g n 1 + el a n el g n el * sqrt : c n 1 + el 0.5 a n el g n el - (* : c n 1 + el 252 -^≤ break 0 if 41 line for φ 2 n 1 -^a n el * x * : m 1 n 1 - (range j n 1 + m - : φ 0.5 φ c j el a j el / φ sin * asin + (* eval : 21 line for φ sin φ cos 1 k φ sin * (2^- sqrt stack eval 51 line :

n ξ EMQF nz n n 2 mod - 2 / : j 1 nz range n 2 mod 1 ≡ (j n 1 + 2 / ≡ (& x j el 0 : v 2 j * 1 - n / ξ 1 -^Kc * ξ 1 -^scdn : x j el v 2 el v 3 el / - eval : 21 line if q j el ξ x j el (2^+ 21 x j el (2^- sqrt * ξ 2^x j el (2^- sqrt * / eval : 21 line for l ξ 2^x j el (2^- (2^(j 1 nz product ξ n 2 n 2 mod * -^1 x j el (2^- (2^(j 1 nz product * / eval : g l l 1 + / sqrt x j el (2^(j 1 nz product ξ n 2 n 2 mod * -^sqrt (/ * eval : x q g 31 sys 51 line :

n 7 :

nz n n 2 mod - 2 / :

ξ 25 :

t n ξ EMQF : X t 1 el : Q t 2 el : g t 3 el : 14 mat

s H s 2^ξ 2^X j el (2^/ + (j 1 nz product s 2^ξ sqrt Q j el / s * + ξ + (j 1 nz product s ξ sqrt + n 2 mod 1 cases * / eval :

f 2 π * 1 - sqrt * 11000 range * 100 / eval :

h 20 g f H vectorize abs vectorize eval * log10 vectorize * :

1 1000 range 100 / h augment