@@ -3,7 +3,7 @@ using DiffEqBiological, DiffEqJump, DiffEqBase, OrdinaryDiffEq, StochasticDiffEq
3
3
r1 = VariableRateReaction (1e-4 ,(1 ,2 ),((1 ,- 1 ),(2 ,1 )))
4
4
r2 = VariableRateReaction (0.01 ,[2 ],[(2 ,- 1 ),(3 ,1 )])
5
5
6
- prob = DiscreteProblem ([999.0 ,1 .0 ,0 ],(0.0 ,250.0 ))
6
+ prob = DiscreteProblem ([999.0 ,3 .0 ,0 ],(0.0 ,250.0 ))
7
7
jump_prob = GillespieProblem (prob,Direct (),r1,r2)
8
8
9
9
jump_prob = JumpProblem (prob,Direct (),build_jumps_from_reaction (r1),build_jumps_from_reaction (r2))
@@ -31,7 +31,7 @@ f = function (t,u,du)
31
31
du[4 ] = u[2 ]* u[3 ]/ 100000 - u[1 ]* u[2 ]/ 100000
32
32
end
33
33
34
- prob = ODEProblem (f,[999.0 ,1 .0 ,0.0 ,100.0 ],(0.0 ,250.0 ))
34
+ prob = ODEProblem (f,[999.0 ,3 .0 ,0.0 ,100.0 ],(0.0 ,250.0 ))
35
35
jump_prob = GillespieProblem (prob,Direct (),r1,r2)
36
36
37
37
srand (230 )
@@ -40,7 +40,7 @@ sol = solve(jump_prob,Tsit5())
40
40
@test 650 <= sol[end ][3 ] <= 900
41
41
42
42
r3 = VariableRateReaction (1e-2 ,[4 ],[(2 ,- 1 ),(3 ,1 )])
43
- prob = ODEProblem (f,[999.0 ,1 .0 ,0.0 ,1.0 ],(0.0 ,250.0 ))
43
+ prob = ODEProblem (f,[999.0 ,3 .0 ,0.0 ,1.0 ],(0.0 ,250.0 ))
44
44
jump_prob = GillespieProblem (prob,Direct (),r1,r2,r3)
45
45
sol = solve (jump_prob,Tsit5 ())
46
46
50
50
51
51
srand (330 )
52
52
println (" Turn Gillespie Problem into SDE" )
53
- prob = SDEProblem (f,g,[999.0 ,1 .0 ,0.0 ,1.0 ],(0.0 ,250.0 ))
53
+ prob = SDEProblem (f,g,[999.0 ,3 .0 ,0.0 ,1.0 ],(0.0 ,250.0 ))
54
54
jump_prob = GillespieProblem (prob,Direct (),r1,r2,r3)
55
55
sol = solve (jump_prob,SRIW1 ())
56
56
@test 650 <= sol[end ][3 ] <= 1500
0 commit comments