trapezoidal_fixed


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

Unless the right hand side of the ODE is linear in the dependent variable, each trapezoid 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().

Licensing:

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

Languages:

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

Related Data and codes:

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

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

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

rk12, an Octave code code which implements Runge-Kutta solvers of orders 1 and 2 for a system of ordinary differential equations (ODE).

rk23, an Octave code code which implements Runge-Kutta ODE solvers of orders 2 and 3.

rk34, an Octave code code which implements Runge-Kutta ODE solvers of orders 3 and 4.

rk4, an Octave code code which applies the fourth order Runge-Kutta (RK) algorithm to estimate the solution of an ordinary differential equation (ODE).

rk45, an Octave code code which implements Runge-Kutta ODE solvers of orders 4 and 5.

rkf45, an Octave code code which implements the Runge-Kutta-Fehlberg ODE solver.

trapezoidal_fixed_test

trapezoidal, an Octave code which solves one or more ordinary differential equations (ODE) using the (implicit) trapezoidal method, and the fsolve() function to handle the implicit equation.

Source Code:


Last revised on 26 April 2021.