This uses the gsl ODE solver.

Files:

• The definition of the ODE solver return value ode.h
• The implementation ode.cpp