# NOT RUN {
library(deSolve)
# MUST use the 'euler' integration method with integer time steps
p <- c(rd=1, alpha=.01)
time <- 0:10
initialN <- 10
out <- ode(y=initialN, times=time,
func=dlogistic, parms=p, method='euler')
plot(time, out[,-1], type='l')
# }
Run the code above in your browser using DataLab