Displaying 1 result from an estimated 1 matches for "ode_sys".
Did you mean:
de_sys
2006 Aug 26
1
problems with loop
...153846)
times3=c(0.115384615,0.125,0.134615385,0.166666667,0.25,0.416666667,0.5)
times4=c(0.58333333,0.666666667,0.75,0.83333333,0.916666667,1)
times5=c( 1.166666667,1.25,1.5,1.6,1.75,2,2.5)
times6=seq(3,20)
times = c(times1,times2,times3,times4,times5,times6)
ODE_sys = function(t, x, w) #w=c("beta0","mu","k")
{
I=x[1]
dI=-w["mu"]*I+w["beta0"]*exp(-w["k"]*t)*(1-I)
list(c(dI),c(S=1-I))...