¿Cuál fue la distancia más lejana a la Tierra que alcanzaron Armstrong, Aldrin y Collins en el Apolo 11?

Este es el comienzo de una respuesta parcial, estoy feliz si alguien se hace cargo o le agrega algo.

¿Por qué no hacer las seis misiones en lugar de solo una?

Usando scripts de estos 1 , 2 , tracé la distancia entre los centros de la Tierra y la Luna como la línea negra, luego las barras de colores gruesas son esas distancias menos el radio de la Tierra, más el radio de la Luna más una altitud del lado lejano de 183 kilómetros

Las longitudes de las barras son las longitudes totales aproximadas de la misión, por lo que tendremos que sumar los tiempos de cada entrada y salida de pasar detrás de la Luna.

No estoy seguro de dónde encontrarlos, pero si alguien puede indicarme la dirección correcta, puedo raspar los números y agregarlos aquí.

from skyfield.api import Topos
from skyfield.api import Loader
import numpy as np
import matplotlib.pyplot as plt

from skyfield.api import load
loaddata = Loader('~/Documents/fishing/SkyData')  # avoids multiple copies of large files

ts = loaddata.timescale() # include builtin=True if you want to use older files (you may miss some leap-seconds)
eph = loaddata('de421.bsp')

earth, moon, venus = [eph[x] for x in ('earth', 'moon', 'venus')]

apollos = [(10, 1969, 5, 18, 26), (11, 1969, 7, 16, 18),
           (12, 1969, 11, 14, 24), (13, 1970, 4, 11, 17),
           (14, 1971, 1, 31, 40), (15, 1971, 7, 26, 38),
           (16, 1972, 4, 16, 27), (17, 1972, 12, 7, 19)]
# https://en.wikipedia.org/wiki/Apollo_program

timez_apollo = []
for n, year, month, d_start, d_stop in apollos:
    times = ts.utc(year, month, range(d_start, d_stop+1))
    timez_apollo.append(times)

days = 1 + np.arange(5*365.2564+1)
times = ts.utc(1969, 1, days)
years = days/365.2564
t_1969 = times.tt[0]

earthpos, moonpos = [x.at(times).position.km for x in (earth, moon)]
r = np.sqrt(((moonpos - earthpos)**2).sum(axis=0))

Rmoon, Rearth = 1737., 6378.
alt = 183.

fig = plt.figure()
ax1 = fig.add_subplot(1, 1, 1)
ax1.plot(years, r, '-k', linewidth=0.5)

for timez in timez_apollo:
    
    yearz = (timez.tt - t_1969) / 365.2564
    
    earthpoz, moonpoz = [x.at(timez).position.km for x in (earth, moon)]
    r = np.sqrt(((moonpoz - earthpoz)**2).sum(axis=0))
    d = r + Rmoon + alt - Rearth

    ax1.plot(yearz, d, linewidth=2.5)
ax1.set_xlim(0.2, 4.0)
ax1.set_ylim(357500, 405500)
ax1.set_ylabel('km')
ax1.set_title('Earth Moon and Apollo far-side orbit to Earth surface distances')
ax1.set_xlabel('years since 1969-01-01')
plt.show()