import math
def meteor(n):
pi = math.pi
ro = 1.2
v = 0
r = .01
h = 1
g = 9.8
m = 21. * .001
t_max = -x/v * 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))
aW1wb3J0IG1hdGgKZGVmIG1ldGVvcihuKToKCXBpID0gbWF0aC5waQoJcm8gPSAxLjIKCXYgPSAwCglyID0gLjAxCgloID0gMQoJZyA9IDkuOAoJbSA9IDIxLiAqIC4wMDEKICAgICAgICB0X21heCA9IC14L3YgKiAxMC4KICAgICAgICBkdCA9IHRfbWF4L24KICAgICAgICBmb3IgaSBpbiByYW5nZShuICsgMSk6CiAgICAgICAgICAgICAgICBkeCA9IHYgKiBkdAogICAgICAgICAgICAgICAgeCA9IHggKyBkeAogICAgICAgICAgICAgICAgaWYgeCA8IDAuOgogICAgICAgICAgICAgICAgICAgICAgICByZXR1cm4gdgogICAgICAgICAgICAgICAgRiA9IHBpLzQgKiBybyAqIHYgKiogMiAqIHIgKiogMiAtIG0gKiBnCiAgICAgICAgICAgICAgICBhID0gRi9tCiAgICAgICAgICAgICAgICBkdiA9IGEgKiBkdAogICAgICAgICAgICAgICAgdiA9IHYgKyBkdgogICAgICAgIHJldHVybiAtOTk5LgoKcHJpbnQobWV0ZW9yKDEwMDAwMCkp