仮想計算機構

IT業界と無縁な派遣社員のブログ

物理学

PythonによるN体シミュレーション

はじめに N体シミュレーションを行いました。 初期値等は以下の情報を参考にしています。 Create Your Own N-body Simulation (With Python) | by Philip Mocz | The Startup | Medium nbody-python/nbody.py at master · pmocz/nbody-python · GitHub 上記…

電気力線みたいな線を描画する

SageMath を用いて電磁気学などで使われる概念についてプロットしてみます。物理学的な厳密性は投げ出しているので差し引いて読んでください。まず、電荷をA~Dと名付け、位置を xs に、電荷の値を qs に格納しておきます。 A = vector([0,0]) B = vector([5,…

SageMathで固有振動のアニメーションを作る

固有振動 n=1,2,...に対して基準振動を下記で定義します。 は弦の長さ、 は弦を伝わる横波の速さを表します。 を変化させた場合のアニメーションをSageMathで作ります。 SageMath - Open-Source Mathematical Software System プログラム main.sageを作成し…