円盤領域で熱方程式や波動方程式を考えるとき、 Bessel 関数が活躍する。 プログラムの正しさの検証 (状況証拠集め) をするためには、 Bessel 関数を計算したくなる。 ところが三角関数のような初等関数と違って、 Bessel 関数は標準のライブラリィには含まれていないことが多い。 Bessel 関数をどうすれば計算できるかについては、
『Bessel 関数の数値計算』を見よ。 中西君の使用言語は Java ということで、 Colt Project のライブラリィを利用した。
http://nalab.mind.meiji.ac.jp/~mk/labo/2004/computing-bessel-function/
…実は最初は Numerical Recipes in C のコードを Java に書き換えて使っていたのだが、 それでは公開が難しいので、探した末に Colt Project に行き着いた。
余談になるが (卒研でもないのだが)、 C で Bessel 関数が必要になった人がいて、 その人には、 大浦卓哉氏の「ベッセル関数 - 整数次」 (http://momonga.t.u-tokyo.ac.jp/~ooura/bessel-j.html) を 勧めておいた。 これ便利ですね。