For computation of the exact solution it is perhaps easiest to implement a function (in say MATLAB or EXCEL) that takes parameters a, b, c, d, F, G and time t as inputs and produces the concentration vector as the output by substituting the appropriate combination of rate constants for the inputs according to kinetic system (2).