Validate Stoichiometry in Mass Action Rate Laws

[longxf21]

Implement functionality to verify that mass action rate laws correctly reflect the stoichiometry of chemical reactions. For example, in the reaction

2A + 3B -> C,

the mass action rate law should be represented as

rate = k * [A]^2 * [B]^3,

where k is the rate constant, and [A] and [B] are the concentrations of reactants A and B, respectively. The exponents 2 and 3 directly correspond to the stoichiometric coefficients of A and B in the chemical equation. This enhancement aims to ensure the accuracy and reliability of simulations involving mass action kinetics by enforcing that the rate laws used in the simulations are consistent with the principles of chemical stoichiometry.