T₁ - время полёта, расчитывается прибавлением dTime₁ к предыдущем значению.
M₁ - полная масса ракеты в начале итерации, берётся из данных или из M₂ предыдущей итерации (строки).
V₁ - скорость ракеты в начале, берётся из данных или из V₂ предыдущей итерации.
S₁ - высота полёта. берётся из данных или вычисляется путём прибавления к предыдущему значению S₁ скорости V₁ умноженной на dTime₁.
F_t1 - тяга на данной высоте ( S₁ ). Вычисляется путём вычитания из тяги в вакууме разницы между двумя тягами домноженной на процент поверхностной плотности воздуха (см. ниже таблицу плотности). F_t1 = F_t1v -(F_t1v -F_t1m) * Ro.
I₁ - удельный импульс на данной высоте ( S₁ ). Вычисляется путём вычитания из импульса в вакууме разницы между двумя импульсами домноженной на процент поверхностной плотности воздуха (см. ниже таблицу плотности). I₁ = I_1v -(I_1v -I_1m) * Ro.
a₁ - ускорение приобретаемое ракетой за счёт двигателей. Вычисляется делением тяги двигателей на массу ракеты.
a₂ - ускорение приобретаемое ракетой за счёт действия сил гравитации. Вычисляется по закону всемироного тяготения. Гравитационная постоянная умножается на массу планеты и делится на квадрат расстояния от ракеты до центра планеты: a₂ = GravPost*M_pl/(R_pl+S₁)².
a₃ - полное ускорение, Вычисляется путём сложения ускорений получаемых от двигателей и гравитации a₃ = a₁ + a₂.
v₂ - скорость в конце итерации. Вычисляется путём сложения скорости в начале итерации и полного ускорения умноженного на промежуток времени v₂ = v₁ + a₃ * dTime₁.
M_t - расход топлива. Вычисляется путём умножения тяги двигателя на промежуток времени и деления на удельный импульс: F_t1*dTime₁/I₁.
M₂ - полная масса ракеты в конце итерации, Вычисляется путём вычитания расхода топлива из массы ракеты в начале итерации. M₂ = M₁ - M_t.

Зависимость плотности воздуха от высоты. Таблица международной стандартной атмосферы (МСА) Высота над уровнем моря, км Плотность, кг/м3 Плотность, % от уровня моря 0 1.250 100% 1 1.134 90.7% 2 1.027 82.2% 3 0.927 74.2% 4 0.836 66.9% 5 0.751 60.1% 6 0.673 53.8% 7 0.601 48.1% 8 0.536 42.8% 9 0.475 38.0% 10 0.421 33.7% 11 0.371 29.7% 12 0.317 25.4% 13 0.271 21.7% 14 0.231 18.5% 15 0.197 15.8% 16 0.169 13.5% 17 0.144 11.5% 18 0.123 9.8% 19 0.105 8.4% 20 0.090 7.2% 21 0.077 6.1% 22 0.065 5.2% 23 0.056 4.5% 24 0.048 3.8% 25 0.041 3.3% 30 0.018 1.44% 35 0.008 0.67% 40 0.004 0.32% 45 0.002 0.16% 50 0.001 0.09% 60 0.0003097 0.02477% 70 0.00008285 0.006628% 80 0.00001846 0.0014768% 90 0.000003418 0.00027344% 100 0.0000005550 0.00004440% 120 0.00000002440 0.000001952%

Расчёт плотности для произвольной высоты Высота, км. : Результат:   ...

Серые клетки в таблице - это не данные МСА, а параметры взятые из книги "Краткий справочник по физике" 1976, Енохович А. С., локальную страницу из которого можно посмотреть по этой ссылке (страница 47). Они нужны, чтобы можно было вычислить плотность до высоты в 120 км.