6 | | - | let MULT18X = parseBigIntValue("1000000000000000000") |
---|
7 | | - | |
---|
8 | | - | let E18X = parseBigIntValue("2718281828459045235") |
---|
9 | | - | |
---|
10 | | - | let iters = [80, 79, 78, 77, 76, 75, 74, 73, 72, 71, 70, 69, 68, 67, 66, 65, 64, 63, 62, 61, 60, 59, 58, 57, 56, 55, 54, 53, 52, 51, 50, 49, 48, 47, 46, 45, 44, 43, 42, 41, 40, 39, 38, 37, 36, 35, 34, 33, 32, 31, 30, 29, 28, 27, 26, 25, 24, 23, 22, 21, 20, 19, 18, 17, 16, 15, 14, 13, 12, 11, 10, 9, 8, 7, 6, 5, 4, 3, 2, 1] |
---|
11 | | - | |
---|
12 | | - | @Callable(i) |
---|
13 | | - | func gamma (zStr18) = { |
---|
14 | | - | let z18 = parseBigIntValue(zStr18) |
---|
15 | | - | let z36 = (z18 * MULT18X) |
---|
16 | | - | func chain (acc,i) = ((z36 + (toBigInt(((i + i) - 1)) * MULT36X)) - fraction((toBigInt((i * i)) * MULT36X), MULT36X, acc)) |
---|
17 | | - | |
---|
18 | | - | let f = { |
---|
19 | | - | let $l = iters |
---|
20 | | - | let $s = size($l) |
---|
21 | | - | let $acc0 = MULT36X |
---|
22 | | - | func $f0_1 ($a,$i) = if (($i >= $s)) |
---|
23 | | - | then $a |
---|
24 | | - | else chain($a, $l[$i]) |
---|
25 | | - | |
---|
26 | | - | func $f0_2 ($a,$i) = if (($i >= $s)) |
---|
27 | | - | then $a |
---|
28 | | - | else throw("List size exceeds 80") |
---|
29 | | - | |
---|
30 | | - | $f0_2($f0_1($f0_1($f0_1($f0_1($f0_1($f0_1($f0_1($f0_1($f0_1($f0_1($f0_1($f0_1($f0_1($f0_1($f0_1($f0_1($f0_1($f0_1($f0_1($f0_1($f0_1($f0_1($f0_1($f0_1($f0_1($f0_1($f0_1($f0_1($f0_1($f0_1($f0_1($f0_1($f0_1($f0_1($f0_1($f0_1($f0_1($f0_1($f0_1($f0_1($f0_1($f0_1($f0_1($f0_1($f0_1($f0_1($f0_1($f0_1($f0_1($f0_1($f0_1($f0_1($f0_1($f0_1($f0_1($f0_1($f0_1($f0_1($f0_1($f0_1($f0_1($f0_1($f0_1($f0_1($f0_1($f0_1($f0_1($f0_1($f0_1($f0_1($f0_1($f0_1($f0_1($f0_1($f0_1($f0_1($f0_1($f0_1($f0_1($f0_1($acc0, 0), 1), 2), 3), 4), 5), 6), 7), 8), 9), 10), 11), 12), 13), 14), 15), 16), 17), 18), 19), 20), 21), 22), 23), 24), 25), 26), 27), 28), 29), 30), 31), 32), 33), 34), 35), 36), 37), 38), 39), 40), 41), 42), 43), 44), 45), 46), 47), 48), 49), 50), 51), 52), 53), 54), 55), 56), 57), 58), 59), 60), 61), 62), 63), 64), 65), 66), 67), 68), 69), 70), 71), 72), 73), 74), 75), 76), 77), 78), 79), 80) |
---|
31 | | - | } |
---|
32 | | - | let exp = pow(E18X, 18, -(z18), 18, 18, HALFUP) |
---|
33 | | - | $Tuple2(nil, ((toString(fraction(exp, MULT36X, f)) + " ") + toString(f))) |
---|
34 | | - | } |
---|
35 | | - | |
---|
36 | | - | |
---|
| 5 | + | unit |
---|