@@ -19,14 +19,14 @@ function bc_model(du, u, h, p, t)
19
19
end
20
20
lags = [tau]
21
21
h (p, t) = ones (3 )
22
- h2 (p, t) = ones (3 ) .- t * q1
22
+ h2 (p, t) = ones (3 ) .- t * q1 * 10
23
23
tspan = (0.0 , 10.0 )
24
24
u0 = [1.0 , 1.0 , 1.0 ]
25
25
prob = DDEProblem (bc_model, u0, h, tspan, constant_lags = lags)
26
26
alg = MethodOfSteps (Vern9 ())
27
27
sol = solve (prob, alg, reltol = 1e-7 , abstol = 1e-10 )
28
28
prob2 = DDEProblem (bc_model, u0, h2, tspan, constant_lags = lags)
29
- sol2 = solve (prob , alg, reltol = 1e-7 , abstol = 1e-10 )
29
+ sol2 = solve (prob2 , alg, reltol = 1e-7 , abstol = 1e-10 )
30
30
31
31
@parameters p0= 0.2 p1= 0.2 q0= 0.3 q1= 0.3 v0= 1 v1= 1 d0= 5 d1= 1 d2= 1 beta0= 1 beta1= 1
32
32
@variables t x₀ (t) x₁ (t) x₂ (.. )
@@ -44,8 +44,8 @@ prob = DDEProblem(sys,
44
44
sol_mtk = solve (prob, alg, reltol = 1e-7 , abstol = 1e-10 )
45
45
@test sol_mtk. u[end ] ≈ sol. u[end ]
46
46
prob2 = DDEProblem (sys,
47
- [x₀ => 1.0 - t * q1, x₁ => 1.0 - t * q1, x₂ (t) => 1.0 - t * q1],
47
+ [x₀ => 1.0 - t * q1 * 10 , x₁ => 1.0 - t * q1 * 10 , x₂ (t) => 1.0 - t * q1 * 10 ],
48
48
tspan,
49
49
constant_lags = [tau])
50
50
sol2_mtk = solve (prob2, alg, reltol = 1e-7 , abstol = 1e-10 )
51
- @test sol2_mtk. u[end ]≈ sol2. u[end ] atol = 1e-5
51
+ @test sol2_mtk. u[end ] ≈ sol2. u[end ]
0 commit comments