Solving a one-dimensional ordinary differential equation
In the next section, we will look at how to solve one-dimensional ordinary differential equations.
Getting ready
We will need to follow some instructions and install the prerequisites.
How to do it...
- As with integration, SciPy has some extremely accurate general-purpose solvers for systems of ordinary differential equations of the first order:
- For real-valued functions, we have basically two flavors:
ode(with options passed with theset_integratormethod) andodeint(simpler interface). The syntax of ode is as follows:
ode(f,jac=None)
- 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 the system (set_integrator), the initial conditions (set_initial_value), and different parameters to be sent to the function or its Jacobian...
