In [1]:
%matplotlib inline
In [2]:
J0: $x0 -> L; k0*x0
# REACTIONS from Spiro, Table 3
# Methylation
#J1_{me}{p}R: T{me}{p}R -> T{me+1}{p} + R; k1c*T{me}{p}R
J1_2pR: T2pR -> T3p + R; k1c*T2pR
J1_3pR: T3pR -> T4p + R; k1c*T3pR
J1_2R: T2R -> T3 + R; k1c*T2R
J1_3R: T3R -> T4 + R; k1c*T3R
In [ ]:
import tellurium as te
model = """
model test
compartment C1;
C1 = 1.0;
species S1, S2;
S1 = 10.0;
S2 = 0.0;
S1 in C1; S2 in C1;
J1: S1 -> S2; k1*S1;
k1 = 1.0;
end
"""
# load models
r = te.loada(model)
res1 = r.simulate(start=0, end=10, points=6)
print res1
In [ ]:
import tellurium as te
myfile = open("chemotaxis.mdl", "r")
contents = myfile.read()
# load models
r = te.loada(contents)
res1 = r.simulate(start=0, end=1000)
print res1