A simple convergent solver for initial value problems

We present a stable and convergent method for solving initial value problems based on the use of differentiation matrices obtained by Lagrange interpolation. This implicit multistep-like method is easy-to-use and performs pretty well in the solution of mildly stiff problems and it can also be applied directly to differential problems in the complex plane.