Warren modeled a two-step electron transfer kinetic system as a reversibly coupled chain of three species. The resulting system of linear differential equations was then solved by Mathematica software to obtain an exact solution. Although effective and useful for investigation of closed-form (i.e., analytic) relationships, the method would be prohibitively tedious if general networks of varying complexity required quick analysis, especially if the network design had to be adjusted by varying a parameter. Automated matrix exponential solutions are particularly well-suited for overcoming this hurdle as e At in the reaction rate constant matrix A can be calculated in MATLAB by the command expm(At) and postmultiplication by a column vector of initial concentrations produces the concentration vector at time t for the network (assuming the process starts at t = 0).