Started on the integration of ODE solver.  Not done yet.
wanted to make the species parameterizations a separate function