“…The open literature reveals that there are several numerical methods, both local and global, devoted to studying reaction networks, but usually with <10 reactions. ,,, These small systems do not exhibit the complexity of most practical chemical reaction systems, where closed analytical expressions are often found embedded in a convolved legacy code. Without the ability to manipulate the complex reaction equations symbolically, it is difficult to apply advanced global methods such as the global terrain method of Lucia and Feng, , the interval method of Kearfott, and the α-BB algorithm of Esposito and Floudas .…”