くだすれFORTRAN（超初心者用）

くだすれFORTRAN（超初心者用） at TECH

185:デフォルトの名無しさん
06/02/20 02:13:21
>>184　
tan(X/2)=1/Yとおくと積分範囲は∞～１

dX=-2/(1+Y^2)dY

cos(X)=(Y^2-1)/(Y^2+1)

1/(1+cos(X)^4) = (Y^2+1)^4/(Y^8+6Y^4+1)

代入して
(1/p)∫exp(-p/tan(X/2)) / (1+cos(X)^4) dX = -(2/p)∫(Y^2+1)^3 exp(-pY) / (Y^8+6Y^4+1) dY
積分範囲を１～∞に反転させれば頭のマイナスが消える。

とりあえず原点付近の１／ｔａｎの発散を避けるようにしてみた積分範囲が無限大にまでになるが、
被積分項は急速に小さくなるので適当なところで打ち切れるだろう。

適当なので間違ってたらスマソ。

次ページ

続きを表示

1を表示