Introduction - If you have any usage issues, please Google them yourself
fortran90 procedures, the use of Runge-Kutta method for numerical solution of ordinary differential equations of the process, is very applicable, including source code, a numerical example and output procedures