ClearAll l = 1.45;(*Länge des Besen*) g = 9.81;(*Fallbeschleunigung*) stepsize = 0.01;(*Zeitschritt \[CapitalDelta]t*) zeit = 0;(*Anfangsbedingung Zeit*) startw = 0.017453;(*Anfangsbedingung Winkel*) geschw = 0;(*Anfangsbedingung Geschwindigkeit*) dd\[CurlyPhi][\[CurlyPhi]_] := (Sin[\[CurlyPhi]]*3*g)/(2*l)(*Zu lösende Differentialgleichung*) Zeitschritt := Module[(*In der Module- Umgebung können wir nun die DGL mit ihren Anfangsbedingungen lösen. Wir erhalten alle Winkel zu allen Zeiten die in dem Bereich unserer \ Anfangsbedingungen liegen in den entsprechenden Zeitschritten.*) {ans, t, \[CurlyPhi], d\[CurlyPhi]}, ans = {{zeit, startw}}; \[CurlyPhi] = startw; t = zeit; d\[CurlyPhi] = geschw; While[ \[CurlyPhi] < 1.5707, d\[CurlyPhi] = d\[CurlyPhi] + dd\[CurlyPhi][\[CurlyPhi]]*stepsize;(*\[CurlyPhi]< 1.5707 muss gelten, da wir uns nur in einen viertel Kreis bewegen, somit 90 Grad.*) \[CurlyPhi] = \[CurlyPhi] + d\[CurlyPhi]*stepsize; t = t + stepsize; ans = Append[ ans, {t, \[CurlyPhi]} ] ]; ans ] loes = Drop[ Zeitschritt, -1];(*Der letzte Eintrag der Liste muss gelöscht \ werden, da er sonst über 90 Grad liegt.*) Last[loes](*Hier geben wir uns den letzten Wert unserer Liste aus. Dies ist der entscheidende Wert für unsere Auswertung, da er die \ Aufprallzeit bei dem entsprechenden Startwinkel ausgibt.*)