Chemical engineering problems often involve the solution of systems of nonlinear algebraic equations (NLE). There are several software packages that can be used for solving NLE systems, but they may occasionally fail, especially in cases where the mathematical model contains discontinuities and/or regions where some of the functions are undefined. Here we describe and demonstrate the use of a systematic procedure for reformulating the NLE system so as to increase the probability of finding a solution and at the same time speeding up the solution process.