From a3068ff812b650792b60321170dc28595f8dbf15 Mon Sep 17 00:00:00 2001 From: noah Date: Thu, 16 Feb 2023 00:00:17 -0600 Subject: [PATCH] Final changes prior to submission --- Schrick-Noah_Homework-3.R | 8 +++++++- 1 file changed, 7 insertions(+), 1 deletion(-) 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]