stoichiometricMatrix -- Computes the stoichiometric matrix of a Reaction Network.

i1 : N = reactionNetwork "A <--> B"

o1 = A-->B
     B-->A

o1 : ReactionNetwork

i2 : S = stoichiometricMatrix N

o2 = | -1 1  |
     | 1  -1 |

              2       2
o2 : Matrix ZZ  <-- ZZ

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 : S = stoichiometricMatrix N

o4 = | -1 1  0  0  0  1  |
     | -1 1  1  0  0  0  |
     | 1  -1 -1 0  0  0  |
     | 0  0  1  -1 1  0  |
     | 0  0  0  -1 1  1  |
     | 0  0  0  1  -1 -1 |

              6       6
o4 : Matrix ZZ  <-- ZZ