Solve ODEs with Python