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 * x/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
t = t + dt
return -999.
print(meteor(100000))
aW1wb3J0IG1hdGgKZGVmIG1ldGVvcihuKToKCXBpID0gbWF0aC5waQoJcm8gPSAxLjIKCXYgPSAwCglyID0gLjAxCgl4ID0gMQoJZyA9IDkuOAoJbSA9IDIxLiAqIC4wMDEKICAgICAgICB0X21heCA9IC0oMiAqIHgvZykgKiogKDEvMikgKiAxMC4KICAgICAgICBkdCA9IHRfbWF4L24KICAgICAgICBmb3IgaSBpbiByYW5nZShuICsgMSk6CiAgICAgICAgICAgICAgICBkeCA9IHYgKiBkdAogICAgICAgICAgICAgICAgeCA9IHggKyBkeAogICAgICAgICAgICAgICAgaWYgeCA8IDAuOgogICAgICAgICAgICAgICAgICAgICAgICByZXR1cm4gdgogICAgICAgICAgICAgICAgRiA9IHBpLzQgKiBybyAqIHYgKiogMiAqIHIgKiogMiAtIG0gKiBnCiAgICAgICAgICAgICAgICBhID0gRi9tCiAgICAgICAgICAgICAgICBkdiA9IGEgKiBkdAogICAgICAgICAgICAgICAgdiA9IHYgKyBkdgoJCXQgPSB0ICsgZHQKICAgICAgICByZXR1cm4gLTk5OS4KCnByaW50KG1ldGVvcigxMDAwMDApKQ==