@@ -19,14 +19,14 @@ function bc_model(du, u, h, p, t)
1919end
2020lags = [tau]
2121h (p, t) = ones (3 )
22- h2 (p, t) = ones (3 ) .- t * q1
22+ h2 (p, t) = ones (3 ) .- t * q1 * 10
2323tspan = (0.0 , 10.0 )
2424u0 = [1.0 , 1.0 , 1.0 ]
2525prob = DDEProblem (bc_model, u0, h, tspan, constant_lags = lags)
2626alg = MethodOfSteps (Vern9 ())
2727sol = solve (prob, alg, reltol = 1e-7 , abstol = 1e-10 )
2828prob2 = 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 )
3030
3131@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
3232@variables t x₀ (t) x₁ (t) x₂ (.. )
@@ -44,8 +44,8 @@ prob = DDEProblem(sys,
4444sol_mtk = solve (prob, alg, reltol = 1e-7 , abstol = 1e-10 )
4545@test sol_mtk. u[end ] ≈ sol. u[end ]
4646prob2 = 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 ],
4848 tspan,
4949 constant_lags = [tau])
5050sol2_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