Adding correction for x0_orbiter

This commit is contained in:
Noah L. Schrick 2023-03-01 21:58:23 -06:00
parent a503e9fffe
commit 220038e604

View File

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