Could odeint be useful in this context? http://headmyshoulder.github.io/odeint-v2/ "Odeint is a modern C++ library for numerically solving Ordinary Differential Equations. It is developed in a generic way using Template Metaprogramming which leads to extraordinary high flexibility at top performance." Best, Arnd