Displaying 1 result from an estimated 1 matches for "dceq1".
Did you mean:
ceq1
2012 Mar 02
0
?Syntax on Taking differential on both sides of the equation in 'R'
...ntial of both sides of the equation and then solve for the inverse of the first as follows (the parameter values are made 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...