stetter_ode_test, a MATLAB code which calls stetter_ode(), which defines an ordinary differential equation (ODE) for which a specific time step sequence causes the implicit trapezoidal rule to be unstable, while the implicit midpoint rule 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, a MATLAB code which defines an ordinary differential equation (ODE) for which a specific time step sequence causes the implicit trapezoidal rule to be unstable, while the implicit midpoint rule will be stable. Note that the right hand side function f(t,y) is periodic, discontinuous, and piecewise linear.