**trapezoidal_fixed**,
a Python code which
solves one or more ordinary differential equations (ODE)
using the (implicit) trapezoidal method, using a fixed point method
for the implicit equation.

Unless the right hand side of the ODE is linear in the dependent variable, each trapezoidal step requires the solution of an implicit nonlinear equation. Such equations can be approximately solved using methods such as fixed point iteration, or an implicit equation solver like fsolve().

The computer code and data files described and made available on this web page are distributed under the GNU LGPL license.

**trapezoidal_fixed** is available in
a MATLAB version and
an Octave version and
a Python version and
an R version.

backward_euler, a Python code which solves one or more ordinary differential equations (ODE) using the backward Euler method.

euler, a Python code which solves one or more ordinary differential equations (ODE) using the forward Euler method.

leapfrog, a Python code which uses the leapfrog method to solve a second order ordinary differential equation (ODE) of the form y''=f(t,y).

midpoint, a Python code which solves one or more ordinary differential equations (ODE) using the midpoint method.

rk4, a Python code which applies the fourth order Runge-Kutta (RK) algorithm to estimate the solution of an ordinary differential equation (ODE) at the next time step.

rkf45, a Python code which implements the Runge-Kutta-Fehlberg (rkf) solver for the solution of a system of ordinary differential equations (ODE).

trapezoidal, a Python code which solves one or more ordinary differential equations (ODE) using the (implicit) trapezoidal method, using fsolve() to handle the implicit equation.

- trapezoidal_fixed.py, the source code.
- trapezoidal_fixed.sh, runs all the tests.
- trapezoidal_fixed.txt, the output file.