method: ‘adams’ or ‘bdf’
Which solver to use, Adams (non-stiff) or BDF (stiff)

with_jacobian : bool
This option is only considered when the user has not supplied a
Jacobian function and has not indicated (by setting either band)
that the Jacobian is banded. In this case, with_jacobian specifies
whether the iteration method of the ODE solver’s correction step is
chord iteration with an internally generated full Jacobian or
functional iteration with no Jacobian.

nsteps : int
Maximum number of (internally defined) steps allowed during one
call to the solver.

first_step : float

min_step : float

max_step : float
Limits for the step sizes used by the integrator.

order : int
Maximum order used by the integrator,
order <= 12 for Adams, <= 5 for BDF.

This integrator is not re-entrant. You cannot have two ode
instances using the “zvode” integrator at the same time.

This integrator accepts the same parameters in set_integrator
as the “vode” solver.

Note

When using ZVODE for a stiff system, it should only be used for
the case in which the function f is analytic, that is, when each f(i)
is an analytic function of each y(j). Analyticity means that the
partial derivative df(i)/dy(j) is a unique complex number, and this
fact is critical in the way ZVODE solves the dense or banded linear
systems that arise in the stiff case. For a complex stiff ODE system
in which f is not analytic, ZVODE is likely to have convergence
failures, and for this problem one should instead use DVODE on the
equivalent real system (in the real and imaginary parts of y).