diff --git a/README.md b/README.md index 3b74e54..72f01b7 100644 --- a/README.md +++ b/README.md @@ -54,8 +54,8 @@ tspan = (0, tf) ps = [1e-4, 1e-2] # Model: non-Markovian part (delay reactions) -τ = 20.0 delay_trigger_affect! = function (integrator, rng) + τ = 20.0 append!(integrator.de_chan[1], τ) # add a delay time τ to the first delay channel end delay_trigger = Dict(1 => delay_trigger_affect!) # the first reaction S+I -> E+I will trigger a delay reaction: E --> I after τ time. @@ -107,10 +107,11 @@ tf = 200.0 tspan = (0, tf) ps = [0.0282, 3.46] # Model: non-Markovian part -τ = 130.0 + delay_trigger_affect! = [] for i in 1:burst_sup push!(delay_trigger_affect!, function (integrator, rng) + τ = 130.0 append!(integrator.de_chan[1], fill(τ, i)) end) end