import math
def meteor(n):
	pi = math.pi
	ro = 1.2
	t = 0
	v = 0
	r = .01
	x = 1
	g = 9.8
	m = 21. * .001
        t_max = -(2 * x/g) ** (1/2) * 10.
        dt = t_max/n
	to = (2 * x/g) ** (1/2)
        for i in range(n + 1):
                dx = v * dt
                x = x + dx
                if x < 0.:
			change = t - to
			print(change)
                        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))

aW1wb3J0IG1hdGgKZGVmIG1ldGVvcihuKToKCXBpID0gbWF0aC5waQoJcm8gPSAxLjIKCXQgPSAwCgl2ID0gMAoJciA9IC4wMQoJeCA9IDEKCWcgPSA5LjgKCW0gPSAyMS4gKiAuMDAxCiAgICAgICAgdF9tYXggPSAtKDIgKiB4L2cpICoqICgxLzIpICogMTAuCiAgICAgICAgZHQgPSB0X21heC9uCgl0byA9ICgyICogeC9nKSAqKiAoMS8yKQogICAgICAgIGZvciBpIGluIHJhbmdlKG4gKyAxKToKICAgICAgICAgICAgICAgIGR4ID0gdiAqIGR0CiAgICAgICAgICAgICAgICB4ID0geCArIGR4CiAgICAgICAgICAgICAgICBpZiB4IDwgMC46CgkJCWNoYW5nZSA9IHQgLSB0bwoJCQlwcmludChjaGFuZ2UpCiAgICAgICAgICAgICAgICAgICAgICAgIHJldHVybiB2CiAgICAgICAgICAgICAgICBGID0gcGkvNCAqIHJvICogdiAqKiAyICogciAqKiAyIC0gbSAqIGcKICAgICAgICAgICAgICAgIGEgPSBGL20KICAgICAgICAgICAgICAgIGR2ID0gYSAqIGR0CiAgICAgICAgICAgICAgICB2ID0gdiArIGR2CgkJdCA9IHQgKyBkdAogICAgICAgIHJldHVybiAtOTk5LgoKcHJpbnQobWV0ZW9yKDEwMDAwMCkp