« moraダウンロード販売で悲愴感が1位に(8/27) | トップページ | 奈良に行ったけど、全く観光できなかったので地図だけでも撮ってきた。 »

2008年8月28日 (木)

ExcelでEmden方程式を解いてみる。

恒星の内部構造についてのエムデン方程式

Emdeneq1

http://scienceworld.wolfram.com/physics/Lane-EmdenEquation.html

http://mathworld.wolfram.com/Lane-EmdenDifferentialEquation.html

を考える。これが解きにくいのはξ=0のとき特異になっていること。いろいろ調べたらξが小さいときの級数解を求めて0から離れたところからスタートすればいいようだ。

久々に手計算した結果はこれ。

Emdeneq2

これを使って、ξ=0.001の値を求めて、あとはDormand-Prince(8次)でnをいろいろ変えて計算。

Emdeneq

n=1の厳密解と比べてもいい精度であっているんでまあまあか。

この前、トーマスフェルミ方程式を解いたけど、あれもこうやって解いたほうが精度がいいのかな?また検討しよう。

« moraダウンロード販売で悲愴感が1位に(8/27) | トップページ | 奈良に行ったけど、全く観光できなかったので地図だけでも撮ってきた。 »

学問・資格」カテゴリの記事

コメント

コメントを書く

(ウェブ上には掲載しません)

トラックバック

この記事のトラックバックURL:
http://app.cocolog-nifty.com/t/trackback/512682/42116670

この記事へのトラックバック一覧です: ExcelでEmden方程式を解いてみる。:

« moraダウンロード販売で悲愴感が1位に(8/27) | トップページ | 奈良に行ったけど、全く観光できなかったので地図だけでも撮ってきた。 »

最近のコメント

2018年11月
        1 2 3
4 5 6 7 8 9 10
11 12 13 14 15 16 17
18 19 20 21 22 23 24
25 26 27 28 29 30  
フォト
無料ブログはココログ