# NOT RUN {
library(deSolve)
p <- c(r=1,alpha=.01, a=10, tau=.1)
time <- 0:10
initialN <- 11
out <- ode(y=initialN, times=time,
func=alogistic, parms=p)
plot(time, out[,-1], type='l')
# }
Run the code above in your browser using DataLab