import numpy as np
import matplotlib.pyplot as plt
def F(P,deltaT):
    l=1.45
    T=0
    g=10
    PP=0
    PPP=np.sin(P)*3*g/(2*l)
    while P < np.pi/2:
        P=P+deltaT*PP
        PP=PP+deltaT*PPP
        PPP=np.sin(P)*3*g/(2*l)
        T=T+deltaT
    return T
Plist=np.linspace(0.1,1.25,50)
Tlist=[F(i,0.001) for i in Plist]
plt.plot(Plist,Tlist)
plt.xlabel ("Anfangswinkel")
plt.ylabel ("Fallzeit")
plt.title ("Fallzeit")
plt.show()