雙伽瑪函數
程式更新日期: 2014年5月23日
程式由網友 Skwai提供。
這個程使用Lanczos approximation的方法去計算雙伽瑪函數的值，能夠準確至有效數字九至十位，另外這個版本較特別之處是可以計算複數因變數的伽瑪函數。
注意: e^( 是按shift ex。
程式 (166 bytes)
?→X: X%%%→M: For X - M→X To X + 2 M Step M: C→B: D→C:
Abs( X→A: A + 3: ( Ans + . 85 )^( A - . 5 ) ÷ e^( Ans )( 1.071370357 +
39.41066934 ÷ A - 35.55136409 ÷ (A + 1) - . 0943944629 ÷ Ans +
6.326951402 ÷ ( A + 2→D: 0>X => - π ÷ ( D X sin( X
πr→D: Next:
( D - B ) ÷ ( 2 M C
限制：只支持約 | X | < 56 . 5（ - 56 . 5 < X < 56 . 5） 的情況。準確度只有五至六位有效數字。
例題: 計算雙伽瑪函數Ψ(0.35)的值。
按 Prog 1 再按 0.35 EXE (顯示答案為-2.971070875)
返回 CASIO fx-50FH 及 fx-50F PLUS 程式集