🚀 go-pugleaf

小口と申します。

<4t7cbd$iks@cmnews.nda.ac.jp>の記事において
asanoa@cc.nda.ac.jpさんは書きました。

>ガウスの誤差関数erf(x)=2/(pi^1/2)*integral(0->x)exp(-x^2)dxを
>数値解析するときのアルゴリズムを教えてください。
>nuemerical recipieという本に載っていると思われますが、日本語のは
>Ｃ言語で、良く理解できません。

ガンマ関数から書いてあるから分からない、のでは？
求めるだけならば、erfcc という関数
（erf の補関数 erfc = 1 - erf を Chebyshev 多項式を当てはめたもの）
を見れば、c の知識がほとんどなくても理解できると思います。
return の所は推測できるでしょう。
# 日本語版の初版で、p.176 の下の方

# 今、 nuemerical recipie を見て、こういう方法があるのか、と知った ^_^;
# 以前に自分で作ったときは、Gause-Helmite 積分で作ったけど...
# 何度も使う、精度があまり要らない、てのなら、
# プログラムの初めに、台形法で計算、配列に次々突っ込む、てのもあるかも
# というか、やったことがある ^_^;;
# （0.01 刻みで、4 位までで十分だったので、400 の配列 ^_^）
   _________________________________________
  /|   SEIKO EPSON CO.  Watch Design Dept.  |\ /|
 /@|<]          Keiichiro Oguchi            | X |
`>_|__oguchi.keiichiro@senet.ha.epson.co.jp_|/ \|