Displaying 1 result from an estimated 1 matches for "dpred".
Did you mean:
pred
2010 Dec 11
2
Predator Prey Models
...0) and
then slowly work its way back to equilibrium. Does anybody know if this
could be done?
LVmod0D <- function(Time, State, Pars) {
with(as.list(c(State, Pars)), {
IngestPred <- rI * N * Pred
GrowthN <- rG * N * (1 - N/K)
MortPred <- rM * Pred
dN <- GrowthN - IngestPred
dPred <- IngestPred * AE - MortPred
return(list(c(dN, dPred)))
})
}
pars <- c(rI = 0.1, rG = 0.9, rM = 0.8, AE = 0.9, K = 20)
yini <- c(N = 20, Pred = 20)
times <- seq(0, 200, by = 1)
y(time=100)<-c(N=10, Pred=5)
print(system.time)
zout <- as.data.frame(lsoda(yini, y(time=1...