您可以使用或作为其参数
set_integrator在
ode类上调用方法。
'dopri5'``'dop853'
这是一个例子：
import numpy as npimport matplotlib.pyplot as pltfrom scipy.integrate import odedef fun(t, z, omega): """ Right hand side of the differential equations dx/dt = -omega * y dy/dt = omega * x """ x, y = z f = [-omega*y, omega*x] return f# Create an `ode` instance to solve the system of differential# equations defined by `fun`, and set the solver method to 'dop853'.solver = ode(fun)solver.set_integrator('dop853')# Give the value of omega to the solver. This is passed to# `fun` when the solver calls it.omega = 2 * np.pisolver.set_f_params(omega)# Set the initial value z(0) = z0.t0 = 0.0z0 = [1, -0.25]solver.set_initial_value(z0, t0)# Create the array `t` of time values at which to compute# the solution, and create an array to hold the solution.# Put the initial value in the solution array.t1 = 2.5N = 75t = np.linspace(t0, t1, N)sol = np.empty((N, 2))sol[0] = z0# Repeatedly call the `integrate` method to advance the# solution to time t[k], and save the solution in sol[k].k = 1while solver.successful() and solver.t < t1: solver.integrate(t[k]) sol[k] = solver.y k += 1# Plot the solution...plt.plot(t, sol[:,0], label='x')plt.plot(t, sol[:,1], label='y')plt.xlabel('t')plt.grid(True)plt.legend()plt.show()这将生成以下图：
