Feature Request boost ode solving engines

[mattfidler]

This allows Adams method, which is used by monolix:

(Adams methods for non stiff ODEs, and Backward Differentiation Formulas methods for stiff ODEs)

The absolute and relative tolerance for the numerical resolutions are 1e-6 and 1e-3 respectively for nonStiff ODEs and 1e-9 and 1e-6 respectively for stiff ODEs .

[mattfidler]

Or odeint

http://headmyshoulder.github.io/odeint-v2/

[mattfidler]

For now, simply keep it and show the difference.