In [1]:
%matplotlib notebook
from libsignetsim import Model, TimeseriesSimulation
In [2]:
m = Model()
m.setName("Enzymatic Reaction")
e = m.listOfSpecies.new("E")
s = m.listOfSpecies.new("S")
p = m.listOfSpecies.new("P")
vmax = m.listOfParameters.new("vmax")
km = m.listOfParameters.new("km")
r = m.listOfReactions.new("Enzymatic reaction")
r.listOfReactants.add(s)
r.listOfModifiers.add(e)
r.listOfProducts.add(p)
r.kineticLaw.setPrettyPrintMathFormula("vmax*E*S/(km+S)", forcedConcentration=True)
e.setValue(10)
s.setValue(12)
p.setValue(0)
vmax.setValue(0.211)
km.setValue(1.233)
In [3]:
sim = TimeseriesSimulation([m], time_min=0, time_ech=1, time_max=20)
sim.run()
sim.plot()
Out[3]: