# Constructing a chemical reaction network from an ode system

**URL:** <https://discourse.julialang.org/t/constructing-a-chemical-reaction-network-from-an-ode-system/132682>\
**Category:** Modelling & Simulations\
**Tags:** catalyst\
**Created:** [September 26, 2025, 5:13pm UTC](https://discourse.julialang.org/t/constructing-a-chemical-reaction-network-from-an-ode-system/132682 "2025-09-26T17:13:02Z")\
**Posts on this page:** 2\
**Page:** 1

<div class="post-metadata">

**Author:** ![adhalanay](https://sea2.discourse-cdn.com/julialang/user_avatar/discourse.julialang.org/adhalanay/32/7371_2.png) [@adhalanay](https://discourse.julialang.org/u/adhalanay)\
**Post date:** [September 26, 2025, 5:13pm UTC](https://discourse.julialang.org/t/constructing-a-chemical-reaction-network-from-an-ode-system/132682/1 "2025-09-26T17:13:02Z")

</div>

The so called Hungarian Lemma states that a system of ODEs with polynomial RHS

x'\_j =f\_j(x),j=\overline{1,n}

comes from a mass-action chemical reaction network if and only if all the negative monomials of f\_f have x\_j among the factors. Of course there are many Chemical Reaction Networks inducing the same system of ODEs. My question is : is there a way to construct one such CRN starting with a system of ODEs satisfying the Hungarian Lemma?

---

<div class="post-metadata">

**Author:** ![ChrisRackauckas](https://sea2.discourse-cdn.com/julialang/user_avatar/discourse.julialang.org/chrisrackauckas/32/77_2.png) [@ChrisRackauckas](https://discourse.julialang.org/u/ChrisRackauckas)\
**Post date:** [September 26, 2025, 8:02pm UTC](https://discourse.julialang.org/t/constructing-a-chemical-reaction-network-from-an-ode-system/132682/2 "2025-09-26T20:02:36Z")

</div>

We don’t have tooling that does this right now, but it would be an interesting symbolic pass that can be written heuristically.