« FNNの夕方のニュースのナレーションって合成音声? | トップページ | RFワールドNo.4の表紙に推薦しておきます。 »

2008年8月31日 (日)

Excelでトーマス・フェルミ方程式を解いてみる(改善版)。

Thomas-Fermi方程式

y'' = t-1/2 y3/2

y(0)=1, y(∞)=0

前に解いたけど、低精度の差分法でちょっと気になっていた。でこの前Emden方程式を解いたときの方法が使えるかと思ってやり直し。

t=0付近での近似解を求めて、そこからスタートすることにしよう。単純な級数では計算できないので、Picard(ピカール。スタートレックのジャン・リュック・ピカードと同じ)の逐次近似でやってみた。手計算は…もう無理。Maximaで計算してみた。y'(0)=cとすると、

y(t)は以下のように展開できる(10次まででも簡単)。

Thomasfermi3

これで t=Δt の時の初期値を得てから、Runge-Kutta8次(Dormand-Prince)で計算した結果がこれ。

Thomasfermi2_2 

フィッティング関数により近くなったな。まあフィッティング関数が正しいというより、今回の方が正しいだろうな。ちょうどt=1くらいのところでフィッティングの分子分母が切り替わるイメージかな。

« FNNの夕方のニュースのナレーションって合成音声? | トップページ | RFワールドNo.4の表紙に推薦しておきます。 »

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

コメント

コメントを書く

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

トラックバック

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

この記事へのトラックバック一覧です: Excelでトーマス・フェルミ方程式を解いてみる(改善版)。:

« FNNの夕方のニュースのナレーションって合成音声? | トップページ | RFワールドNo.4の表紙に推薦しておきます。 »

最近のコメント

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  
フォト
無料ブログはココログ