RM = reactantMatrix N
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.