The solution of steady state chemical reactor models, using classic numerical methods, (like Newton-Raphson, Quasi-Newton, etc.), is often very difficult. It is proposed to select constraints for the independent variables as the first step of the solution, then to use a combination of the Newton-Raphson and continuation methods, to obtain the solution, starting from a feasible initial guess. Several methods for the solution of steady reactor models are compared. It is shown that the proposed method is the most reliable and efficient for locating the solution.