今回は、Pythonによる常微分方程式を解く際によく使われるscipyライブラリのメモである。 Scipyによる常微分方程式 Pythonを用いて常微分方程式を解く際は、ライブラリを用いるのが最も手軽で、かつ高精度な結果が得られる。 常微分方程式を解くライブラリと…

常微分方程式の数値計算例題-ロジスティック方程式 ここでは、常微分方程式の数値計算の例題として、ロジスティック方程式を扱ってみる。 ロジスティック方程式 ロジスティック方程式は、生物の個体数の増加を表すモデルであり、以下のような形の方程式とし…

本稿から具体的な数値計算の手法についての話に入る。 まずは、古典力学の基本方程式であるNewtonの運動方程式を数値的に解く手法から。 古典力学と常微分方程式 古典力学は、基本的には質点の運動はNewtonの運動方程式で記述される。 Newtonの運動方程は、…

Matplotlib基本事項メモ ここでは、数値計算結果をグラフなどの可視化の際に便利なMatplotlibライブラリの基本的事項を記しておく。 Matplotlib Matplotlibはグラフ描画用のライブラリである。数値計算の結果を簡易的に確認するだけではなく、論文やプレゼン…

現在、以下の書籍で数値計算（特に物理に関する）におけるPythonの活用を学習している。 Pythonによる計算物理作者:大槻純也森北出版Amazon まずその１として、Pythonで数値計算コードを書く際に威力を発揮するNumpyライブラリごく基本的な事項のメモである…

ここでは、M1 Macにhomebrew、pyenvを使ったpython環境を構築した際のメモを記す。 基本的には以下のサイトを参考にした。 https://zenn.dev/unsoluble_sugar/articles/283ebf698c307c homebrewのアップデート Terminal上で以下のコマンドを打ってhomebrewを…

はじめに マルチフィジックスシミュレーションのCAEソフトウェアとして代表的なものとしてCOMSOLなどがあるが、オープンソースのもので該当するものがないか探してみた。オープンソースとして比較的歴史があるものとして、Elmerがある。Elmerはフィンランド…

Python リスト内包表記 基本的な書き方 リスト内包表記は、リストを生成する時の一つの方法である。例えば、1～10までのリストを作成する基本的な書式は以下のようになる。 list = [] for i in range(1,11): list.append(i) list [1, 2, 3, 4, 5, 6, 7, 8, 9…

gist32e2f928c41146efff4265dda7f5f32a

Pythonの勉強として、数値計算のプログラムを題材にしているが、その時に2次元のタプルでの記述が出てきたので、備忘録としてメモ。 例えば、２次元平面内に二つの電荷の位置座標と電荷の値をセットで定義をするのに以下のように記述する。 具体的には、電荷…