biochemical reaction networks. We are interested in finding polynomial

conditions on the parameters such that the steady state ideal of a chemical

reaction network is binomial under every specialisation of the parameters if

the conditions on the parameters hold. We approach the binomiality problem

using Comprehensive Gr\\"obner systems. Considering rate constants as

parameters, we compute comprehensive Gr\\"obner systems for various reactions.

In particular, we make automatic computations on n-site phosphorylations and

biomodels from the Biomodels repository using the grobcov library of the

computer algebra system Singular.

},\n}\n'