MATLAB Implementation of Adams-Bashforth-Moulton ODE Solver

Adams-Bashforth-Moulton Variable-Step-Size Predictor-Corrector Numerical Integration of a System of Ordinary Differential Equations (ODEs)
This method solves the first-order system of ODE's of the following form:
Y' = F(t,Y(t)) with a <= t <= b and Y(a) = alpha where Y = mx1 vector and Y(a) = mx1 vector
The function "F" is evaluated using the user-defined function called "odefunc":
f(t,y(t)) = odefunc(t,y(t)) where y(t) = mx1
The user can also define a custom termination condition