I don't propose to go into the maths in my case, unless someone wants to point me towards it, but I'd like to ask a question that's on my mind: I know that it is possible to write a precise differential equation for any 1 equilibrium which (after integration) will yield a result of being able to calculate the concentration of any species involved in the equilibrium at any time, given the forward and backward rate constants, equilibrium constant, and initial concentrations of all species.

Is it possible, in theory, to write a differential equation for *any* large *system* of equilibria, which can then be solved to calculate the concentration of any species involved in the system at a certain time, given the rate constants for every equilibrium, every equilibrium constant, and the initial concentrations of all species involved?