i1 : N = reactionNetwork "A + 2B <--> 3A + 4B"
o1 = A+2B-->3A+4B
3A+4B-->A+2B
o1 : ReactionNetwork
|
i2 : RM = reactantMatrix N
o2 = | 1 2 |
| 3 4 |
2 2
o2 : Matrix ZZ <--- ZZ
|
A bigger example:
i3 : N = oneSiteModificationA()
o3 = S_0+E-->X
X-->S_0+E
X-->E+S_1
S_1+F-->Y
Y-->S_1+F
Y-->S_0+F
o3 : ReactionNetwork
|
i4 : RM = reactantMatrix N
o4 = | 1 1 0 0 0 0 |
| 0 0 1 0 0 0 |
| 0 0 1 0 0 0 |
| 0 0 0 1 1 0 |
| 0 0 0 0 0 1 |
| 0 0 0 0 0 1 |
6 6
o4 : Matrix ZZ <--- ZZ
|
The object reactantMatrix is a method function.