import math
def meteor(n):
pi = math.pi
ro = 1.2
v = 0
r = .01
x = 1
g = 9.8
m = 21. * .001
t_max = -(2 * h/g) ** 1/2 * 10.
dt = t_max/n
for i in range(n + 1):
dx = v * dt
x = x + dx
if x < 0.:
return v
F = pi/4 * ro * v ** 2 * r ** 2 - m * g
a = F/m
dv = a * dt
v = v + dv
return -999.
print(meteor(100000))
aW1wb3J0IG1hdGgKZGVmIG1ldGVvcihuKToKCXBpID0gbWF0aC5waQoJcm8gPSAxLjIKCXYgPSAwCglyID0gLjAxCgl4ID0gMQoJZyA9IDkuOAoJbSA9IDIxLiAqIC4wMDEKICAgICAgICB0X21heCA9IC0oMiAqIGgvZykgKiogMS8yICogMTAuCiAgICAgICAgZHQgPSB0X21heC9uCiAgICAgICAgZm9yIGkgaW4gcmFuZ2UobiArIDEpOgogICAgICAgICAgICAgICAgZHggPSB2ICogZHQKICAgICAgICAgICAgICAgIHggPSB4ICsgZHgKICAgICAgICAgICAgICAgIGlmIHggPCAwLjoKICAgICAgICAgICAgICAgICAgICAgICAgcmV0dXJuIHYKICAgICAgICAgICAgICAgIEYgPSBwaS80ICogcm8gKiB2ICoqIDIgKiByICoqIDIgLSBtICogZwogICAgICAgICAgICAgICAgYSA9IEYvbQogICAgICAgICAgICAgICAgZHYgPSBhICogZHQKICAgICAgICAgICAgICAgIHYgPSB2ICsgZHYKICAgICAgICByZXR1cm4gLTk5OS4KCnByaW50KG1ldGVvcigxMDAwMDApKQ==