Implement interpolation routines
Provide routines to interpolate the reaction rate coefficients and density to the current time step of the ODE solver.
The interpolation routines need to be called in Fcn() and Jacobian() routines.
User should be able to choose between no interpolation / nearest neighbour, linear or spline / polynomial interpolation.
Depends on #10 (closed)