This paper presents a simple method for solving multireaction chemical equilibrium problems via Gibbs Minimization through the use of a stochastic global optimizer. The method is explained at the level of upperundergraduate chemistry or chemical engineering students. Its matlab code is also presented and explained. Four examples with varying difficulties illustrate the versatility and robustness of the method, which is simple, general, and scalable and can be easily extended to simultaneous chemical and phase equilibrium problems. ß 2016 Wiley Periodicals, Inc. Comput Appl Eng Educ 24:899-904, 2016; View this article online at wileyonlinelibrary.com/journal/ cae;