Solving Ordinary Differential Equations in Python