Octaveでの実装を参考にCommon Lispで正規分布の累積分布関数を計算 - 丸井綜研

（2023-01-23追記：正規分布の計算を見直してみました →正規分布の確率密度関数と累積分布関数の計算 - 丸井綜研） 正規分布の確率密度関数は定義通りに ;; Normal Distribution ;; (normpdf 2 0 1) => 0.05399096581690089 ;; (defun normpdf (x mu sigma) (* (/ 1 (* (sqrt (* 2 pi)) sigma)) (exp (- (/ (expt (- x mu) 2) (* 2 sigma sigma)))))) で計算できるのですが、累積分布関数は確率密度関数を定積分したものなので、やっかいでなかなかプログラムが作れずにいました。 Octaveではどういう実装になっているのかと、normcdfのMファイルを読み辿っていったら、誤差関数(erf)を使って計算していることが分かりました

[MagnesiumRibbon]

残念ながら誤差関数は標準ではないですねー/添削したくなるコードだ…オーバーヘッドすごい