  An algorithm is given in this paper for the computation of dynamically
equivalent weakly reversible realizations with the maximal number of reactions,
for chemical reaction networks (CRNs) with mass action kinetics.
