odeint is a highly flexible library for solving ordinary differential equations.
odeint - solving ordinary differential equations in c++ v2