ExcelでEmden方程式を解いてみる。
恒星の内部構造についてのエムデン方程式
http://scienceworld.wolfram.com/physics/Lane-EmdenEquation.html
http://mathworld.wolfram.com/Lane-EmdenDifferentialEquation.html
を考える。これが解きにくいのはξ=0のとき特異になっていること。いろいろ調べたらξが小さいときの級数解を求めて0から離れたところからスタートすればいいようだ。
久々に手計算した結果はこれ。
これを使って、ξ＝0.001の値を求めて、あとはDormand-Prince（8次）でnをいろいろ変えて計算。
n=1の厳密解と比べてもいい精度であっているんでまあまあか。
この前、トーマスフェルミ方程式を解いたけど、あれもこうやって解いたほうが精度がいいのかな？また検討しよう。
« moraダウンロード販売で悲愴感が１位に（8/27） | トップページ | 奈良に行ったけど、全く観光できなかったので地図だけでも撮ってきた。 »
「学問・資格」カテゴリの記事
- 高周波（RF・マイクロ波・ミリ波・5G）関連ニュース2020年1月16日 5G、レーダで使われるミリ波フェイズドアレイアンテナ各社まとめ、HuaweiはSAWまで取り込む、iPhone 11 ProがFBIにロック解除される、など。(2020.01.16)
- 高周波（RF・マイクロ波・ミリ波・5G）関連ニュース2020年1月15日 iPhone 5G mmWのAiP、TSMCとともにASEにも発注と報道、Microwave Journal 1月号の特集はレーダ＆アンテナ、ResonantのWhitepaperなど。(2020.01.15)
- ”楕円の円周と一方の径からもう一方の径を求める。”をカシオの高精度計算サイトkeisan.casio.jpに自作式としてUP.(2020.01.13)
- 高周波（RF・マイクロ波・ミリ波・5G）関連ニュース2020年1月13日 IMECがbeyond 5GのRF半導体技術（III-V/CMOS)について解説、ミンチークオはiPhone 5GモデルはSub6GHzもミリ波も同時に出ると語る、など。(2020.01.13)
- 高周波（RF・マイクロ波・ミリ波・5G）関連ニュース2020年1月4日 Qorvoが6GHz対応のWiFiを支持＆5G RF for dummies 2nd edを公開、TSMCがiPhone mmW向けにInFO AiPの注文を受ける、Akouticsが軍用レーダ向けBAW出荷など。(2020.01.04)
« moraダウンロード販売で悲愴感が１位に（8/27） | トップページ | 奈良に行ったけど、全く観光できなかったので地図だけでも撮ってきた。 »
コメント