@@ -39,11 +39,11 @@ dna_rs = @reaction_network begin
3939 k6, DNAR --> DNA + P
4040end
4141rates = [:k1 => 0.5 ,
42- :k2 => (20 * log (2.0 ) / 120.0 ),
43- :k3 => (log (2.0 ) / 120.0 ),
44- :k4 => (log (2.0 ) / 600.0 ),
45- :k5 => 0.025 ,
46- :k6 => 1.0 ]
42+ :k2 => (20 * log (2.0 ) / 120.0 ),
43+ :k3 => (log (2.0 ) / 120.0 ),
44+ :k4 => (log (2.0 ) / 600.0 ),
45+ :k5 => 0.025 ,
46+ :k6 => 1.0 ]
4747tf = 1000.0
4848u0 = [:DNA => 1 , :mRNA => 0 , :P => 0 , :DNAR => 0 ]
4949prob = DiscreteProblem (dna_rs, u0, (0.0 , tf), rates, eval_module = @__MODULE__ )
@@ -103,7 +103,7 @@ oscil_rs = @reaction_network begin
103103 0.05 , SP2 --> 0
104104end
105105u0 = [:X => 200.0 , :Y => 60.0 , :Z => 120.0 , :R => 100.0 , :S => 50.0 , :SP => 50.0 ,
106- :SP2 => 50.0 ] # Hill equations force use of floats!
106+ :SP2 => 50.0 ] # Hill equations force use of floats!
107107tf = 4000.0
108108prob = DiscreteProblem (oscil_rs, u0, (0.0 , tf), eval_module = @__MODULE__ )
109109"""
@@ -121,7 +121,7 @@ specs_sym_to_name = Dict(:S1 => "R(a,l)",
121121 :S8 => " A(Y~P,r!1).R(a!1,l)" ,
122122 :S9 => " A(Y~P,r)" )
123123rsi = Dict (:R0 => 1 , :L0 => 2 , :A0 => 3 , :kon => 4 , :koff => 5 ,
124- :kAon => 6 , :kAoff => 7 , :kAp => 8 , :kAdp => 9 )
124+ :kAon => 6 , :kAoff => 7 , :kAp => 8 , :kAdp => 9 )
125125params = (5360 , 1160 , 5360 , 0.01 , 0.1 , 0.01 , 0.1 , 0.01 , 0.1 )
126126rs = @reaction_network begin
127127 kon, S1 + S2 --> S4
146146rates = [:kon , :kAon , :koff , :kAoff , :kAp , :kAdp ] .=>
147147 params[[rsi[:kon ], rsi[:kAon ], rsi[:koff ], rsi[:kAoff ], rsi[:kAp ], rsi[:kAdp ]]]
148148u0 = [:S1 => params[1 ], :S2 => params[2 ], :S3 => params[3 ], :S4 => 0 , :S5 => 0 ,
149- :S6 => 0 , :S7 => 0 , :S8 => 0 , :S9 => 0 ]
149+ :S6 => 0 , :S7 => 0 , :S8 => 0 , :S9 => 0 ]
150150tf = 100.0
151151prob = DiscreteProblem (rs, u0, (0.0 , tf), rates, eval_module = @__MODULE__ )
152152"""
@@ -220,7 +220,7 @@ rn = @reaction_network begin
220220 c8, P2G --> P2 + G
221221end
222222rnpar = [:c1 => 0.09 , :c2 => 0.05 , :c3 => 0.001 , :c4 => 0.0009 , :c5 => 0.00001 ,
223- :c6 => 0.0005 , :c7 => 0.005 , :c8 => 0.9 ]
223+ :c6 => 0.0005 , :c7 => 0.005 , :c8 => 0.9 ]
224224varlabels = [" G" , " M" , " P" , " P2" , " P2G" ]
225225u0 = [:G => 1000 , :M => 0 , :P => 0 , :P2 => 0 , :P2G => 0 ]
226226tf = 4000.0
0 commit comments