stetter_ode

stetter_ode, a Python code which sets up and solves an ordinary differential equation (ODE) for which a specific time step sequence causes the implicit trapezoidal method to be unstable, while the implicit midpoint method will be stable. The right hand side function f(t,y) is periodic, discontinuous, and piecewise linear.

Licensing:

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

Languages:

stetter_ode is available in a MATLAB version and an Octave version and a Python version.

Reference:

• Hans Joerg Stetter,
  Analysis of discretization methods for ordinary differential equations,
  Springer, 1973,
  ISBN: 9783642654718.

Related Data and codes:

python_ode, Python codes which sets up various systems of ordinary differential equations (ODE).

Source Code:

• stetter_ode.py, the source code.
• stetter_ode.sh, runs all the tests.
• stetter_ode.txt, the output file.

• stetter_backward_euler.png,
• stetter_euler.png,
• stetter_lambda.png,
• stetter_midpoint.png,
• stetter_trapezoidal.png,