はじめに 数値解析のレポートで、修正オイラー法や、4次のルンゲクッタ法の安定領域の図示にかなり時間食われたので、今後そんなことが無いようにここに書き残して置くことにしました。 やりたいこと 今回自分がやりたかったことは、修正オイラー法の安定条…

はじめに 学校のレポートをやるために、$ z^{ n } + a_{ n - 1 } z^{ n - 1 } + \cdots + a_0 = 0 \, ( z \in \mathbb{ C } ) $という多項式代数方程式の全ての根を数値的に求める必要が出たのでPythonで実装しました。 Newton法とか、1個だけ求めるアルゴリ…

はじめに 脳磁逆問題を解く際に複素線積分する必要があっただが、Pythonで複素線積分をしているものが見つからなかったので、記事を書こうと思いました。 まあ、複素積分って実質は線積分に過ぎないので、線積分が数値的に求められたら、それを複素数にして…

数値解析のレポートをやります。これが終われば春休みになったも同然です。はい。 ただ、授業聞いてなかったので、まずはInputからです。はい。 問題はこれ 問.調和振動子の初期値問題: $$ { \displaystyle \frac{d}{dt} \left( \begin{array}{c} p \\\ q \e…