|
208 | 208 | end
|
209 | 209 |
|
210 | 210 | (p, dX), 0 <--> X
|
211 |
| - (p, dY), 0 <--> Y |
| 211 | + (1.1*p, dY), 0 <--> Y |
212 | 212 | end
|
213 | 213 |
|
214 | 214 | # Creates model programmatically.
|
|
218 | 218 | rxs = [
|
219 | 219 | Reaction(p, nothing, [X], nothing, [1])
|
220 | 220 | Reaction(dX, [X], nothing, [1], nothing)
|
221 |
| - Reaction(p, nothing, [Y], nothing, [1]) |
| 221 | + Reaction(1.1*p, nothing, [Y], nothing, [1]) |
222 | 222 | Reaction(dY, [Y], nothing, [1], nothing)
|
223 | 223 | ]
|
224 | 224 | continuous_events = [
|
|
236 | 236 | # Tests that approaches yield identical results.
|
237 | 237 | @test isequal(rn_dsl, rn_prog)
|
238 | 238 |
|
239 |
| - u0 = [X => 5.0, Y => 3.0, Z => 3.5] |
| 239 | + u0 = [X => 6.0, Y => 4.5, Z => 5.5] |
240 | 240 | tspan = (0.0, 20.0)
|
241 |
| - ps = [p => 1.0, dX => 0.05, dY => 0.05, dY_up => 0.01] |
| 241 | + ps = [p => 0.5, dX => 0.025, dY => 0.025, dY_up => 0.01] |
242 | 242 |
|
243 | 243 | sol_dsl = solve(ODEProblem(rn_dsl, u0, tspan, ps), Tsit5())
|
244 | 244 | sol_prog = solve(ODEProblem(rn_prog, u0, tspan, ps), Tsit5())
|
|
0 commit comments