search for: dveq1

Displaying 1 result from an estimated 1 matches for "dveq1".

Did you mean: veq1
2012 Mar 02
0
?Syntax on Taking differential on both sides of the equation in 'R'
...ade up): library(deSolve) rkMethod("rk45dp7") CSIeq1<-function(t,yeq1,pars) { with (as.list(c(yeq1,pars)),{ Neq1<-Seq1+Ceq1+Ieq1 dSeq1<-d[a*(1-Neq1)/(f*Veq1+m+d)] dCeq1<-d[(f*Seq1*Veq1+g*Ieq1+r*(1-Neq1)-b1*Veq1*Ieq1)/(b2+m+d+g)] dIeq1<-d[(-b2*Ceq1)-r*(1-Neq1)/(b1*Veq1-g-u)] dVeq1<-d[o*(Ceq1+Ieq1)/e] return(list(c(Seq1,Ceq1,Ieq1,Veq1),Neq1)) }) }pars <- c(a=0.1, m=0.0005, u=0.00005, b1=0.7, b2=0.2, f=0.01, g=0.4, r=0.3, o=20, e=90, d=0.01) # initial conditions yeq1 <- c(Seq1=0.99,Ceq1=0.01,Ieq1=0,Veq1=1) t <-seq(0,365, by=50) o1 <- ode(yeq1, t,...