The flux for the reaction A+B <-> C+D is given by:
J0(t) = k1*A(t)*B(t) - k2*C(t)*D(t)
Where k1 , k2 are rate constants for the forward and backward reactions, respectively.
Load the model below and then hit the 'Start' button to run the simulation. Adjust the parameter sliders (k1, k2, ) below the plot to see how they affect reactant and product concentrations as a function of time.