stiff_ode

stiff_ode, an Octave code which sets up an ordinary differential equation (ODE) y'=lambda*(cos(t)-y), which exhibits stiffness, making accurate solution difficult.

Licensing:

The information on this web page is distributed under the MIT license.

Languages:

stiff_ode is available in a C version and a C++ version and a Fortran90 version and a MATLAB version and an Octave version and a Python version and an R version.

Related Data and codes:

stiff_ode_test

octave_ode, an Octave code which sets up various ordinary differential equations (ODE).

stiff_exact, an Octave code which returns an exact solution y(t) to a stiff ordinary differential equation (ODE) of the form y' = lambda * ( cos(t-t0) - y).

Reference:

• John D Cook,
  Stiff differential equations,
  https://www.johndcook.com/blog/
  Posted 02 February 2020.

Source Code:

• stiff_euler_backward.m, uses the backward Euler method to solve the stiff ODE.
• stiff_euler.m, uses the Euler method to solve the stiff ODE.
• stiff_exact.m, evaluates the exact solution of the stiff ODE.
• stiff_midpoint.m, uses the midpoint method to solve the stiff ODE.
• stiff_deriv.m, evaluates the right hand side of the stiff ODE.
• timestamp.m, prints the YMDHMS date as a timestamp.