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 = -(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))