ODE system solver made simple. For IVPs (initial value problems).
soypat
Fixed and variable-step Runge-Kutta solvers in Modern Fortran
jacobwilliams