import mumot

mumot.about()

model1 = mumot.parseModel(r"""
U -> A : g_1
A -> U : a_1
A + U -> A + A : r_1
""")

model1.showODEs()

model2 = model1.substitute('U = N - A')

modelStreamPlot = model2.stream('A', showFixedPoints = True)