diff --git a/Schrick-Noah_Homework-3.R b/Schrick-Noah_Homework-3.R index 4260361..021eeda 100644 --- a/Schrick-Noah_Homework-3.R +++ b/Schrick-Noah_Homework-3.R @@ -282,11 +282,17 @@ legend(x = "topleft", # c) Use k3=0.02 pp.params.c <- c(.01, .1, .02, .05) -pp.sol.c <- ode45(f = function(t,y){pp.f(t,y,k=pp.params)}, +pp.sol.c <- ode45(f = function(t,y){pp.f(t,y,k=pp.params.c)}, y = c(prey0, pred0), t0 = tmin, tfinal = tmax) plot.pred.prey(pp.sol.c, "ODE45 with k3=0.02") +pp.params.c2 <- c(.01, .1, .02, .75) +pp.sol.c2 <- ode45(f = function(t,y){pp.f(t,y,k=pp.params.c2)}, + y = c(prey0, pred0), + t0 = tmin, tfinal = tmax) +plot.pred.prey(pp.sol.c2, "ODE45 with k3=0.02, k4=0.75") + ## 5. Solve SIR model numerically from t=0 -> 20 sir.f <- function(t, y, k) { dS <- -k[1] * y[1] * y[2]