例題で学ぶ微分方程式：1.1 常微分方程式と相図

なんとなく「時代は微分方程式だよなあ」と思って先日オライリーの「例題で学ぶ微分方程式」を買いました。
www.oreilly.co.jp

本書ではMathematicaでプログラムが書かれているので、これをPythonでというかJupyter notebook (scipy+matplotlib) で置き換えながら読み進めようと思います。

というか、jupyter notebookをgist経由でブログに埋め込むテストです。

例題1.1.1の図版の再現

こういう図を再現します。

SciPyで常微分方程式の数値解法を行うには、scipy.integrate.ode または scipy.integrate.odeint を使います。