Introduction - If you have any usage issues, please Google them yourself
To solve the coupled dynamical system, we use the initial condition and the 4th order Runge-Kutta method to find the first 3 values, and so to apply Adams-Moulton method to
obtain the approximation solution of the system