From 220038e6040c5e90c6c869457bdff148e9df4191 Mon Sep 17 00:00:00 2001 From: noah Date: Wed, 1 Mar 2023 21:58:23 -0600 Subject: [PATCH] Adding correction for x0_orbiter --- Schrick-Noah_Homework-4.R | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/Schrick-Noah_Homework-4.R b/Schrick-Noah_Homework-4.R index d41e96b..8b8d86b 100644 --- a/Schrick-Noah_Homework-4.R +++ b/Schrick-Noah_Homework-4.R @@ -385,7 +385,7 @@ orbiter.estimate <- findZeroRelax(orbiter.f, 1000000) l <- orbiter.estimate[1]-r_earth x0_orbiter_tmp <- c(l, orbiter.f(l), 0)/1e3 # km -x0_orbiter <- x0_orbiter + x0_earth +x0_orbiter <- x0_orbiter_tmp + x0_earth y.init.orbiter <- c(x0_orbiter, v0_earth) sunorbit.sol <- rk4sys(f=sunearth.f, a=0, b=365.25, y0=y.init.orbiter, n=365)