|
|
十進BASICがインストールされたフォルダにある
Libraryフォルダに1000桁モード用の超越関数計算ルーチンのサンプルをおいてあります。
計算結果の正確さの評価などはしてないので推奨できるレベルのものではありませんが，
とりあえず使ってみることはできると思います。
10進1000桁モードで実行してください。
DECLARE EXTERNAL FUNCTION EXP
PRINT EXP(1)
END
! 1000桁モードで利用する指数関数
EXTERNAL FUNCTION EXP(x)
FUNCTION s(y,n)
LET t=y*x/n
IF ABS(t)<=EPS(0) THEN
LET s=y+t
ELSE
LET s=y+s(t,n+1)
END IF
END FUNCTION
LET EXP=s(1,1)
END FUNCTION
補足
カイ２乗分布の定義式にΓ関数がでてきますが，nが自然数のときのΓ(n+1/2)は簡単に計算できるので，
あまり面倒なことをしなくてすむと思います。
|
|