「万有引力の2体問題の運動方程式を数値的に解く前の下ごしらえ」の改良版。変数の規格化をもう少しシンプルに。
万有引力の2体問題（は結局1体問題に帰着するんだけど）の運動方程式を数値的に解く前の下ごしらえ。闇雲な初期条件からはじめるのではなく，そもそも数値計算しなくても楕円になることはわかっているのだから，ルンゲ・クッタ法なりを使って数値的に解く前に，それなりの下ごしらえをしておこう。
楕円軌道上の時刻ごとの位置を Python で数値的に求める
「楕円軌道上の時刻ごとの位置を求めるための下ごしらえ」で下ごしらえした式を Python で数値的に解く。
続きを読む
楕円軌道上の時刻ごとの位置を Maxima で数値的に求める
「楕円軌道上の時刻ごとの位置を求めるための下ごしらえ」で下ごしらえした式を Maxima で数値的に解く。
続きを読む
楕円軌道上の時刻ごとの位置を数値的に求めるための下ごしらえ
楕円軌道は $r(\phi)$ で表されるので，$\phi$ を $t$ の関数として数値的に求めるための下ごしらえ。 続きを読む
Matplotlib で作成した mp4 ファイルは Windows の Edge では見ることができない
授業で「Matplotlib でアニメーション」をやってて気づいたのですが，作成した mp4 ファイルを Microsoft Edge で見ようとしてもできないようです。Mac の Edge では見えるし，同じ Windows でも Firefox なら見える。 続きを読む
SymPy Plotting Backends で電場ベクトルの方向場を描く
正負の点電荷がつくる電場の向きを表す方向場を Python の SymPy Plotting Backends (SPB) で描く。
Matplotlib で電場ベクトルの方向場を描く
正負の点電荷がつくる電場の向きを表す方向場を Python の Matplotlib で描く。
続きを読む
Python で作成したグラフの部分だけを pdf として保存する：アップデート版
以前，「Python で作成したグラフの部分だけを pdf として保存する」のアップデート版。最近の弘大 JupyterHub の環境では，（matplotlib のバージョンアップにより）特にエラーが出ることもなく，簡単に pdf として保存できる。
続きを読む
Weinberg の教科書にある CMB の「双極的異方性」の式の食い違い？
Weinberg の2つの教科書
- Gravitation and Cosmology (1972)
- Cosmology (2008) （日本語訳「ワインバーグの宇宙論」(2013)）
に書いてある宇宙背景放射 CMB の「双極的異方性」の式が異なっている件。