Odeint - Solving ordinary differential equations in C++