Chemical reaction networks (CRNs) are susceptible to mathematical modelling. The dynamic behavior of CRNs can be investigated by solving the polynomial equations derived from its structure. However, simple CRN give rise to non-linear polynomials that are difficult to resolve. Here we propose a procedure to locate the steady states of CRNs from a formula derived through algebraic geometry methods. We have applied this procedure to define the steady states of a classic CRN that exhibits instability, and to a model of programmed cell death.