  In this paper we introduce a new representation for the multistationarity
region of a reaction network, using polynomial superlevel sets. The advantages
of using this polynomial superlevel set representation over the already
existing representations (cylindrical algebraic decompositions, numeric
sampling, rectangular divisions) is discussed, and algorithms to compute this
new representation are provided. The results are given for the general
mathematical formalism of a parametric system of equations and so may be
applied to other application domains.
