楕円積分の解き方（マクローリン展開？）

解き方をすっかり忘れてしまいました…
下の楕円積分を解いて下さい。
よろしくお願いします。

No.1 ベストアンサー 回答者： info22_ 回答日時： 2010/01/15 16:06

kを変数のままマクローリン展開しても長ったらしいkの無限級数

展開式を求めても意味はありません。
E(k)=1.570796326794897-0.39269908169743*k^2-0.073631401025032*k^4
-0.030399798903041*k^6-0.038626564642975*k^8+0.41887581315451*k^10
-3.255408021537378*k^12- ...
K(k)=1.570796326794897+0.39269908169743*k^2+0.2208942030751*k^4
+0.1519989945152*k^6+0.27038595250082*k^8-3.76988231839055*k^10
+35.80948823691115*k^12- ...

最初の式は第2種完全楕円積分 E(k) の定義式そのものです。
なのでk(0≦k≦<1)を与えて数値積分しかないです。
下の積分サイトでkを与えて計算ボタンをおせば積分値をだしてくれます。
http://keisan.casio.jp/has10/SpecExec.cgi?path=0 …
例えば　k=0.5　→　E(k)=1.46746220933942715546

2番目の積分は第1種完全楕円積分 K(k) の定義式そのものです。
なのでk(0≦k≦<1)を与えて数値積分しかないです。
下の積分サイトでkを与えて計算ボタンをおせば積分値をだしてくれます。
http://keisan.casio.jp/has10/SpecExec.cgi?path=0 …
例えば　k=0.5　→　K(k)=1.685750354812596042871

この回答へのお礼

ご回答ありがとうございました。
ご紹介いただいたURLも参考になりました。

お礼日時：2010/01/18 10:03

No.3 回答者： Ae610 回答日時： 2010/01/16 06:34

「解け」・・・の意味がよく分からないので、取り敢えず

以下のように級数展開は出来る。

∫[0,π/2]｛1/√(1-(ksinθ)^2)｝dθ=π/2・F(1/2,1/2,1;k^2)

∫[0,π/2]｛√(1-(ksinθ)^2)｝dθ=π/2・F(-1/2,1/2,1;k^2)

(F(a,b,c;x)=(Γ(c)/(Γ(a)Γ(b)))・｛Σ(n=0～∞)(Γ(a+n)Γ(b+n)x^n)/(n!・Γ(c+n))｝なる超幾何関数)

この回答へのお礼

ご回答ありがとうございました。

お礼日時：2010/01/18 10:05

No.2 回答者： spring135 回答日時： 2010/01/15 23:13

楕円関数の本をみても級数展開のような式はなかなか出てきません。

私ならば数値積分をやります。シンプソンの台形公式を使って、0～π/2を10分割ぐらいすると、8ケタぐらいは簡単に出ます。

この回答へのお礼

ご回答ありがとうございました。
本が無いため、なかなか難しい所ではありますが、何とか調達してみます。

お礼日時：2010/01/18 10:04