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")