powered by
ode
deSolve
predpreyLV(t, y, params)
Stevens. M.H.H. 2009. A Primer of Ecology with R. Use R! Series. Springer.
predpreyRM
params1 <- c(b=.5, a=.01, s=.2, e=.1) Time <- seq(0,100, by=.1) # Set time here LV.out <- ode(c(H0=25,P0=5), Time, predpreyLV, params1) matplot(Time, (LV.out[,2:3]), type="l", ylab="Population Size")
Run the code above in your browser using DataLab