Hossa 
Oder anders gefragt, wie würdest du das zweite Axiom auf
dieses Besipiel ummünzen!
Das zweite Newton’sche Axiom wird häufig als „Kraft gleich Masse mal Beschleunigung“ wiedergegeben:
\vec F=m\cdot\vec a
Das findet man auch so in vielen Büchern. Jedoch gilt diese Gleichung nur für den häufig vorkommenden Spezialfall, dass die Masse m konstant ist. Newton selbst hat die Gleichheit von Kraft F und Änderung des Impulses p=m*v mit der Zeit gefordert:
\vec F=\frac{d\vec p}{dt}\quad=\frac{d}{dt}\left(m\vec v\right)=\frac{dm}{dt},\vec v+m,\frac{d\vec v}{dt}=\dot m,\vec v+m,\vec a
Zum Zeitpunkt t habe die Rakete die Masse m(t). Sie erfährt dann die Kraft:
\vec F_r=m(t)\cdot\vec a
Das mit konstanter Geschwindigkeit vg ausgestoßene Gas erfährt die Kraft:
\vec F_g=-\dot m(t)\cdot\vec v_g
Dazu überlegt man sich, dass die Masse der Rakete von m(t) auf m(t+dt) abnimmt. Es wird also die Gasmenge [m(t)-m(t+dt)] ausgestoßen. Das ist gerade gleich -dm, wenn man dm wie üblich als [m(t+dt)-m(t)] definiert. Die Ausstoß-Geschwindigkeit ist vg. Der Impuls ist dann -dm*vg. Die Kraft ist die Zeitableitung des Impulses, also (-dm/dt)*vg.
Wegen des dritten Axioms „action = reactio“ sind beide Kräfte entgegengesetzt gleich:
\vec F_r=-\vec F_g\quad\left|,\mbox{einsetzen}\right.
m\cdot\vec a=-\left(-\dot m\cdot\vec v_g\right)\quad\left|,\mbox{umschreiben}\right.
m,\frac{d\vec v}{dt}=\frac{dm}{dt},\vec v_g\quad\left|,\cdot\frac{dt}{m}\right.
d\vec v=\vec v_g,\frac{dm}{m}\quad\left|,\int\cdots\right.
\int\vec dv=\vec v_g,\int\frac{1}{m},dm
\vec v-\vec v_0=\vec v_g,\left[\ln(m)+C\right]
Die Integrationskonstante v0 ist die Anfangsgeschwindigkeit der Rakete. Die Integrationskonstante C ergibt sich aus der Randbedingung, dass sich die Geschwinidgkeit der Rakete zu Beginn der Verbrennung noch nicht geändert hat:
\vec 0=\vec v\left(m_0\right)-\vec v_0=\vec v_g,\left[\ln\left(m_0\right)+C\right]\quad\Longrightarrow\quad C=-\ln\left(m_0\right)
Damit ergibt sich für die Geschwindigkeit der Rakete:
\vec v(m)=\vec v_0+\vec v_g,\left[\ln(m)-\ln\left(m_0\right)\right]=\vec v_0+\vec v_g,\ln\left(\frac{m}{m_0}\right)
Wegen m0 ist der Logarithmus negativ, um ihn positiv zu machen, kann man den Kehrwert nehmen und dafür ein Minuszeichen vor die Logarithmusfunktion setzen:
\vec v(m)=\vec v_0-\vec v_g,\ln\left(\frac{m_0}{m}\right)
Setzt man darin die Anfangsgeschwindigkeit der Rakete v0=0, erkennt man schön, dass vg und die Geschwinidigkeitsänderung der Rakete entgegengesetzt gerichtet sind.
Viele Grüße
Hasenfuß