### Simulation Freier Fall

## Konstanten:
g =	       9.81 # m/s² (Fallbeschleunigung Erde)

## Anfangsbedingungen:
h_init = 100    # m    (Anfangshöhe)
v_init =   0    # m/s  (Anfangsgeschwindigkeit)

## Genauigkeit:
dt =       0.2  # s    (Zeitschritt)

list_t = []
list_y = []
y  =  h_init
vy =  v_init
ay = -g
t  =  0
while y > 0:
    list_t.append(t)
    list_y.append(y)
    ay = -g
    y  += vy * dt # Neues y aus dem alten Wert von vy berechnen.
    vy += ay * dt # Neues vy aus dem alten Wert von ay berechnen.
    t  += dt

## Textausgabe:
print("Falldauer: ", round(t,3), " s")
print("Endgeschwindigkeit:", round(vy,3), "m/s")

## Diagramm:
import matplotlib.pyplot as plt
plt.plot(list_t, list_y, 'ro')
plt.ylabel('Höhe (m)')
plt.xlabel('Zeit (s)')
plt.title('Simulation: Freier Fall aus '+str(h_init)+'m')
plt.show()