%O3

constants
mu = G*M_E;

T = 97.4*60

a = (T^2*mu/(4*pi^2))^(1/3)

h = a-R_E

v = visviva_circ(a) % Velocity of satellite.

alpha = asin(R_E/a*sin(deg2rad(95)))

beta = 180-rad2deg(alpha)-95 % See drawn diagram
beta = deg2rad(beta);

arclength = 2*beta*a % See drawn diagram

omega_s = v/a % Angular velocity of satellite.
% t = arclength/v  % Time for non-rotating earth.
t = arclength/(omega_s*a + omega_E*a)% E->W rotation.
t = arclength/(omega_s*a - omega_E*a)% W->E rotation.
