CHW Example 2.3 (Example 1 of AS08)

Loads the chemJac command file, which contains all commands used below

\(\(\( (*\ \ SetDirectory["C:\CHEMmma"]\ \ *) \)\(
\)\(<< "chemJac.m"\)\)\)

The RHS of the ODE has the following form.
Here the functions k1,....,k6 are increasing in each argument when the other arguments are fixed.

f = {-k1[x1] + k2[-x1, -x3], -k3[x2] + 
          k4[-x2, x1], -k5[x3] + k6[-x3, x2]}

The list of species concentrations is denoted by ' vars'.

vars = {x1, x2, x3}

The Jacobian of the RHS of the ODE, and its Craciun-Feinberg determinant

j = jac[f,  vars]; 
j // MatrixForm

det = Det[j]

coeffs[det]