File: mixed-volume-experiment.m2

package info (click to toggle)
macaulay2 1.21%2Bds-3
  • links: PTS, VCS
  • area: main
  • in suites: bookworm
  • size: 133,096 kB
  • sloc: cpp: 110,377; ansic: 16,306; javascript: 4,193; makefile: 3,821; sh: 3,580; lisp: 764; yacc: 590; xml: 177; python: 140; perl: 114; lex: 65; awk: 3
file content (13 lines) | stat: -rw-r--r-- 373 bytes parent folder | download | duplicates (4) 
1
2
3
4
5
6
7
8
9
10
11
12
13
 
restart
needsPackage "ReactionNetworks"
needsPackage "PHCpack"
CRN = reactionNetwork "A <--> 2B, A + C <--> D, B + E --> A + C, D --> B + E"
stoichiometricSubspace CRN
F = steadyStateEquations CRN
R = ring first F
describe R
R' = CC[gens R]
toR' = map(R',R, gens R')
monoms = F / monomials / toR'
F' = apply(monoms, m->(m*random(CC^(numcols m),CC^1))_(0,0))
mixedVolume F'