For purpose of solving a general first-order network via matrix exponential it is straightforward to handle the homogeneous system ⃗ x = A⃗ x directly and compute ⃗ x(t) = e A(t-t0) ⃗ x 0 . The reaction rate constant matrix A can be found as follows. For a network consisting of only species i and j are coupled reversibly i