wikipediaのローレンツ方程式に書いてあるそのままのパラメータで解いてみた。画像はgnuplotで描画したもの。 同じものをいくつか他の解法で解いて比較すると面白そう。 # coding: utf-8 # lorenz attractor # scheme: euler fx = lambda x,y,z,r,p,b: -p*x …