This method, which is an indirect method, involves first processing the rules of the model specification (S1 File) to obtain the reaction network implied by the rules, as well as the corresponding ordinary differential equations (ODEs) for mass-action kinetics.