Compare commits
2 Commits
4f39845b57
...
717e29c9f7
| Author | SHA1 | Date | |
|---|---|---|---|
| 717e29c9f7 | |||
| 220038e604 |
@ -385,7 +385,7 @@ orbiter.estimate <- findZeroRelax(orbiter.f, 1000000)
|
|||||||
l <- orbiter.estimate[1]-r_earth
|
l <- orbiter.estimate[1]-r_earth
|
||||||
|
|
||||||
x0_orbiter_tmp <- c(l, orbiter.f(l), 0)/1e3 # km
|
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)
|
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)
|
sunorbit.sol <- rk4sys(f=sunearth.f, a=0, b=365.25, y0=y.init.orbiter, n=365)
|
||||||
|
|||||||
Loading…
x
Reference in New Issue
Block a user