From 0eb34e59ac0fbbe93ea3830c8f6e6e0110826ce8 Mon Sep 17 00:00:00 2001 From: noah Date: Wed, 15 Feb 2023 14:27:04 -0600 Subject: [PATCH] Adding dt=h for Euler decay solution --- Schrick-Noah_Homework-3.R | 9 ++++++--- 1 file changed, 6 insertions(+), 3 deletions(-) diff --git a/Schrick-Noah_Homework-3.R b/Schrick-Noah_Homework-3.R index 16aa751..65dc567 100644 --- a/Schrick-Noah_Homework-3.R +++ b/Schrick-Noah_Homework-3.R @@ -89,15 +89,18 @@ k <- 0.03 y0 <- 100 decay.euler.sol <- my_euler(f=function(t,y){decay.f(t,y,k=k)}, - y0, t0, tfinal, dt=1) -plot(decay.euler.sol$t,decay.euler.sol$y) + y0, t0, tfinal, dt=h) +plot(decay.euler.sol$t,decay.euler.sol$y, xlab="t", ylab="y", + main="Euler Numerical Solution for the decay ode dy/dt=-ky") + par(new=T) lines(seq(t0,tfinal,len=50), y0*exp(-k*seq(t0,tfinal,len=50)), col="red") decay.rk4.sol <- my_rk4(f=function(t,y){decay.f(t,y,k=k)}, y0, t0, tfinal) -plot(decay.rk4.sol$t,decay.rk4.sol$y) +plot(decay.rk4.sol$t,decay.rk4.sol$y, xlab="t", ylab="y", + main="RK4 Numerical Solution for the decay ode dy/dt=-ky") par(new=T) lines(seq(t0,tfinal,len=50), y0*exp(-k*seq(t0,tfinal,len=50)), col="red")