Im Prinzip hast du schon die Lösung für die ganze Facharbeit aufgezeigt
Für alle die sich Tipparbeit sparen wollen:
kg = 1; m = 1; sek = 1; grad = Pi/180; d = 45*10^-4;
M = 547/1000000 kg; (*Masse Projektil*)
g = 981/100 m/sek^-2 ; (*Erdbeschleunigung*)
p = 1293/ 1000 kg m^-3 ; (*Luftdichte*)
J = 75/10 kg m^2 sek^-2; (*Energie in Joule*)
v0 = Sqrt[2 J/M]; (*Mündungs-Geschwindigkeit*)
a = 0.0 grad ; (*Winkel*)
h = 2 m; (*Höhe*)
ch = -2/100; cv = ch ; (*negativer Widerstandsbeiwert *)
Ah = Pi*d^2/4 m^2; Av = Ah; (*Fläche*)
tm = 2 sek; (*Flugzeit*)
x1 = NDSolve[{ch*Ah*p*rx'[t]^2/2 == M*rx''[t], rx'[0] == v0*Cos[a],
rx[0] == 0}, rx[t], {t, 0, tm}];
x2 = NDSolve[{-M*g + Sign[ry'[t]]* cv*Av*p*0.5 *ry'[t]^2 ==
M*ry''[t], ry'[0] == v0*Sin[a], ry[0] == h}, ry[t], {t, 0, tm}];
ListPlot[Table[
Flatten[Evaluate[{rx[t] /. x1, ry[t] /. x2}]], {t, 0, tm/2,
0.01 sek}], AspectRatio -> 1/3, Frame -> True, GridLines -> Automatic,
FrameLabel -> {"range [m]", "height [m]", "bullet trajectory"} ]