2017分子動力学による粒子のシミュレーション粒子の動きを分子動力学で計算しmatplotlibで可視化してみる。各粒子間には相互作用は働かないシンプルな系として考える。速度の計算はベルレ法(参考)。各粒子はお互いが見えていないので、初速のまま真っ直ぐに飛び続ける。境界を越えたら反対側へ。dpi=50)実行の結果は下記。gifファイルサイズを小さくするため下記サイトを使ってみた。https://compressor.io/結果は下記。余白が多いようなgifアニメは結構圧縮が効く様子。2017cygwinでmatplotlibとscipyを使えるようにする↓こちらのサイトを参考に、そのままの手順でできました。pip共にpython用外部ライブラリのパッケージ管理ツールeasy_installの方が古いらしく、pipはこれらを置き換える目的で開発されている。2017任意の数のランダムなxy座標を生成するリスト形式でappendしてやることで、指定数を追加できる。2017関数に指定した引数は参照渡しPythonは関数に指定した引数がすべて参照渡しとなるらしい。指定した引数の型によっては変更不可となる。20id関数を使ってみると、関数内でイミュータブル変数を変更すると識別子が変わった。最近の投稿ある点(x,y)が多角形の内部にあるかどうかを判定するアルゴリズムボロノイ図を使いポリゴン図形の中に中心線を引くpolygonを形成する頂点列からボロノイ図を作成する頂点列で表すPolygonを構成する長いエッジに対しエッジ分割し、頂点を追加する。