Displaying 1 result from an estimated 1 matches for "zout".
Did you mean:
out
2010 Dec 11
2
Predator Prey Models
...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=100), times, LVmod0D, pars))
plot(zout[,1],zout[,2], type="l",col="black", ylim=range(0,25),
xlab="time",ylab="population")
points(zout[,1],zout[,3],type="l",col="red")
title(main="prey (black...