noah/CS-7863-Sci-Stat-Proj-3
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

Adding extra plot features for ODE prey solution comparisons

Browse Source
This commit is contained in:
Noah L. Schrick 2023-02-15 23:21:05 -06:00
parent 21c4dfe69c
commit 7222500638
1 changed files with 10 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

13
Schrick-Noah_Homework-3.R
View File

@ -268,10 +268,17 @@ plot.pred.prey(pp.rk4.sol, "RK4")
# plot comparing Prey solutions to ode45
plot(pp.sol$t, pp.sol$y[,1], type="l", col="blue")
plot(pp.sol$t, pp.sol$y[,1], type="l", col="green", lwd=3, ylim=c(20,120),
xlim=c(0,200), main="Comparing Prey solutions")
par(new=T)
lines(pp.euler.sol$t,pp.euler.sol$y[,1], col="red")
lines(pp.euler.sol$t,pp.euler.sol$y[,1], col="red", lwd=2)
par(new=T)
lines(pp.rk4.sol$t,pp.rk4.sol$y[,1], col="blue", lty="dashed", lwd=2)
legend(x = "topleft",
legend = c("ODE45", "Euler", "RK4"),
lty = c(1, 1, 2),
col = c("green", "red", "blue"),
lwd = 2)
# c) Use k3=0.02
pp.params.c <- c(.01, .1, .02, .05)