stetter_ode, an Octave 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. Note that the right hand side function f(t,y) is periodic, discontinuous, and piecewise linear.
The information on this web page is distributed under the MIT license.
stetter_ode is available in a MATLAB version and an Octave version and a Python version.
octave_ode, an Octave code which sets up various ordinary differential equations (ODE).