H. T. Davis Γa(x)

e90820e8-7be9-432a-bbee-06e7cb2c7a0c 22 4 5 decimal

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

H. T. Davis Γa(x)

Abramowitz & Stegun, p. 256 c 1 0.5772156649015329 0.6558780715202538 - 0.0420026350340952 - 0.1665386113822915 0.0421977345555443 - 0.0096219715278770 - 0.0072189432466630 0.0011651675918591 - 0.0002152416741149 - 0.0001280502823882 0.0000201348547807 - 0.0000012504934821 - 0.0000011330272320 0.0000002056338417 - 0.0000000061160950 0.0000000050020075 0.0000000011812746 - 0.0000000001043427 0.0000000000077823 0.0000000000036968 - 0.0000000000005100 0.0000000000000206 - 0.0000000000000054 - 0.0000000000000014 0.0000000000000001 261 mat :

1 z Γ / c k el z k^* k 1 ∞ sum ≡

z Γpoly c k el z k^* k 126 sum eval :

z Γa 1 z Γpoly / :

x c 10 x ≤ (x 1 ≤ (& cases :

x Γa x c * 1 1 + 15 red 1 5 mat 2 1 sys if/otherwise ... (x ≤ 0)&(x ≤ 15) suite IterativeProduct

Gamma(x)

constant 15 D accuracy @ integers ... can be made iterated product(s) x Γ x Γa x 0 ≤ x Γa 0 x ≤ (x 1 ≤ (& x 1 - (x 1 - Γa * 1 x ≤ (x 2 ≤ (& x 2 - (x 1 - (* x 2 - Γa * 2 x ≤ (x 3 ≤ (& x 3 - (x 2 - (* x 1 - (* x 3 - Γa * 3 x ≤ (x 4 ≤ (& x 4 - (x 3 - (* x 2 - (* x 1 - (* x 4 - Γa * 4 x ≤ (x 5 ≤ (& x 5 - (x 4 - (* x 3 - (* x 2 - (* x 1 - (* x 5 - Γa * 5 x ≤ (x 6 ≤ (& x 6 - (x 5 - (* x 4 - (* x 3 - (* x 2 - (* x 1 - (* x 6 - Γa * 6 x ≤ (x 7 ≤ (& x 7 - (x 6 - (* x 5 - (* x 4 - (* x 3 - (* x 2 - (* x 1 - (* x 7 - Γa * 7 x ≤ (x 8 ≤ (& x 8 - (x 7 - (* x 6 - (* x 5 - (* x 4 - (* x 3 - (* x 2 - (* x 1 - (* x 8 - Γa * 8 x ≤ (x 9 ≤ (& x 9 - (x 8 - (* x 7 - (* x 6 - (* x 5 - (* x 4 - (* x 3 - (* x 2 - (* x 1 - (* x 9 - Γa * 9 x ≤ (x 10 ≤ (& x 10 - (x 9 - (* x 8 - (* x 7 - (* x 6 - (* x 5 - (* x 4 - (* x 3 - (* x 2 - (* x 1 - (* x 10 - Γa * 10 x ≤ (x 11 ≤ (& x 11 - (x 10 - (* x 9 - (* x 8 - (* x 7 - (* x 6 - (* x 5 - (* x 4 - (* x 3 - (* x 2 - (* x 1 - (* x 11 - Γa * 11 x ≤ (x 12 ≤ (& x 12 - (x 11 - (* x 10 - (* x 9 - (* x 8 - (* x 7 - (* x 6 - (* x 5 - (* x 4 - (* x 3 - (* x 2 - (* x 1 - (* x 12 - Γa * 12 x ≤ (x 13 ≤ (& x 13 - (x 12 - (* x 11 - (* x 10 - (* x 9 - (* x 8 - (* x 7 - (* x 6 - (* x 5 - (* x 4 - (* x 3 - (* x 2 - (* x 1 - (* x 13 - Γa * 13 x ≤ (x 14 ≤ (& x 14 - (x 13 - (* x 12 - (* x 11 - (* x 10 - (* x 9 - (* x 8 - (* x 7 - (* x 6 - (* x 5 - (* x 4 - (* x 3 - (* x 2 - (* x 1 - (* x 14 - Γa * 14 x ≤ (x 15 ≤ (& x 15 - (x 14 - (* x 13 - (* x 12 - (* x 11 - (* x 10 - (* x 9 - (* x 8 - (* x 7 - (* x 6 - (* x 5 - (* x 4 - (* x 3 - (* x 2 - (* x 1 - (* x 15 - Γa * cases :

Γ(x)=> usable as is up to x = 15 Higher numerical: use iterated product

x Γ

1/1 Mathcad sanity up to 15 14.987654321 Γ 84347399235.0809

x Γ ext'd iter'd product n x trunc : x j - (j 1 n product x n - Γa * eval 31 line :

Numerical Γ(x) . . .

Sanity vs Mathcad => . . . not exhausted

]]>

29.6 Γ 2.2897321973268 1030^*

Γ(x) [15 D] Alvaro

Γ.κ 0.99999999999999709182 57.156235665862923517 59.597960355475491248 - 14.136097974741747174 0.49191381609762019978 - 0.33994649984811888699 104 -^* 0.46523628927048575665 104 -^* 0.98374475304879564677 104 -^* - 0.15808870322491248884 103 -^* 0.21026444172410488319 103 -^* - 0.21743961811521264320 103 -^* 0.16431810653676389022 103 -^* - 0.84418223983852743293 104 -^* 0.26190838401581408670 104 -^* - 0.36899182659531622704 105 -^* 151 mat :

ζ Γ z ζ 1 - : zh 1 - 2 ζ * + 2 / : zg zh 607128 / + : zp zg 0.5 zh *^: 14 mat 2 π * sqrt Γ.κ 1 el Γ.κ κ 1 + el z κ + / κ 114 sum + (* zp 2^* e zg -^* eval 21 line :

Valid for Re(ζ)>0

14.987654321 Γ 84347399235.0768

29.6 Γ 2.2897321973267 1030^*