Automated symbolic computation of the ODE system's Jacobian

Expected to significantly improve solver's speed as compared to using a different, faster solver.

ODE are in CRN from: dx/dt = f(x) = N*v(x), which, for the Jacobian @f/@x requires first computation of derivatives of parametric rates functions @v/@x. These should be coded in C and interfaced from R, as the ODEs are.

Est: 4d