Fourth order approximation of the time derivatives of metabolite concentrations and calculation of the Jacobian.