「コンピュータ演習」の授業,今年最後の「13日」を記念した特別企画。
本日,2023年12月13日は(水曜日だけど)今年最後の「13日」。また,今年最初の「13日」,2023年1月13日は金曜日であった。いわゆる「13日の金曜日」。
なぜ「13日の金曜日」が特別に扱われるのかについては,素養が必要であるので,各自調査すること。ここでは,Wikipedia をあげておく。
「コンピュータ演習」の授業,今年最後の「13日」を記念した特別企画。
本日,2023年12月13日は(水曜日だけど)今年最後の「13日」。また,今年最初の「13日」,2023年1月13日は金曜日であった。いわゆる「13日の金曜日」。
なぜ「13日の金曜日」が特別に扱われるのかについては,素養が必要であるので,各自調査すること。ここでは,Wikipedia をあげておく。
「万有引力の2体問題の運動方程式を数値的に解く前の下ごしらえ ver. 2」で下ごしらえした式を Python で数値的に解く。
「万有引力の2体問題の運動方程式を数値的に解く前の下ごしらえ ver. 2」で下ごしらえした式を Maxima で数値的に解く。 続きを読む
「万有引力の2体問題の運動方程式を数値的に解く前の下ごしらえ」の改良版。変数の規格化をもう少しシンプルに。
万有引力の2体問題(は結局1体問題に帰着するんだけど)の運動方程式を数値的に解く前の下ごしらえ。闇雲な初期条件からはじめるのではなく,そもそも数値計算しなくても楕円になることはわかっているのだから,ルンゲ・クッタ法なりを使って数値的に解く前に,それなりの下ごしらえをしておこう。
「楕円軌道上の時刻ごとの位置を求めるための下ごしらえ」で下ごしらえした式を Python で数値的に解く。
続きを読む
「楕円軌道上の時刻ごとの位置を求めるための下ごしらえ」で下ごしらえした式を Maxima で数値的に解く。
続きを読む
楕円軌道は $r(\phi)$ で表されるので,$\phi$ を $t$ の関数として数値的に求めるための下ごしらえ。 続きを読む
授業で「Matplotlib でアニメーション」をやってて気づいたのですが,作成した mp4 ファイルを Microsoft Edge で見ようとしてもできないようです。Mac の Edge では見えるし,同じ Windows でも Firefox なら見える。 続きを読む
正負の点電荷がつくる電場の向きを表す方向場を Python の SymPy Plotting Backends (SPB) で描く。