% R11

constants

%---
a1 = 6880e3;  % periapsis
a2 = 42250e3; % apoapsis

m0 = 2000;
ve = 3e3;
f = 10; % mass flow.
%---

a = (a1+a2)/2

v1 = sqrt(mu/(a1)) % parking circ. orbit
vGTO1 = sqrt(mu*(2/a1 - 1/a)) % GTO at periapsis.
vGTO2 = sqrt(mu*(2/a2 - 1/a)) % GTO at apoapsis
v2 = sqrt(mu/(a2)) % Geo circ. orbit.

deltaV1 = vGTO1 - v1

m1 = m0/exp(deltaV1/ve) 

deltam1 = m0-m1
tb1 = deltam1/f

deltaV2 = v2 - vGTO2

m2 = m1/exp(deltaV2/ve)
deltam2 = m1 - m2

tb2 = deltam2/f