Tak, to są wartości między etapem spadku a eksplozji. Kod poniżej obejmuje dwa pierwsze etapy (wznoszenie i opadanie). W każdym kroku albo (1) silnik jest włączony zużywając 0.01 masy paliwa albo (2) jest wyłączony nie zużywając paliwa. Błędy zmiennoprzecinkowe dają niedokładne wartości. m1 to masa paliwa, m0 to całkowita masa rakiety. n* to tymczasowe zmienne.
Czy to co piszę ma sens?
[bbone=python,2381]
for t in arange(t0, t2, step):
if m1 > 0:
nm0 = m0 - step
nm1 = m1 - step
f = grav + Fa*v + thrust
else:
nm0 = m0
nm1 = m1
f = grav + Fa*v
nv = v + f*(step/m0)
np = p + v*step
p,v,m0,m1 = np,nv,nm0,nm1[/bbone]