As with integration, SciPy has some extremely accurate general-purpose solvers for systems of ordinary differential equations of first order.
For the case of real-valued functions we have basically two flavors –
ode (with options passed with the
set_integrator method) and
odeint (simpler interface). The syntax of
ode is as follows:
The first parameter,
f, is the function to be integrated, and the second parameter,
jac, refers to the matrix of partial derivatives with respect to the dependent variables (the Jacobian). This creates an
ode object, with different methods to indicate the algorithm to solve ...