デバイ関数（Debye's function）をkeisan.casio.jpにUP!

先日に引き続き統計力学編。デバイ関数ですよ。

リンクはこちら↓。

デバイ関数(Debye's function)

デバイの比熱式に出てくるデバイ関数
fD(x) = (3/y^3)∫[0,y]exp(x)*x^4 / (exp(x)-1)^2 dx
を計算します。

計算にはシンプソンの積分公式を使っているんだけれど、二重積分公式のほうがいいかも。誰か書き直してくれれば（めんどくさい）。