correct, working reaction networks. still missing equality on
multisets, and in particular MassActionReaction