language: Python (python 2.7.2)
date: 102 days 18 hours ago
link:
visibility: public 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23

import math
def meteor(n):
        r = .09
        m = 7.
        b = 1.5e-3
        x = 200. * 1000.
        v = -20. * 1000.
        H = 7.6 * 1000.
        g = 9.8
        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 = b * v ** 2 * math.exp(-x/H) - m * g
                a = F/m
                dv = a * dt
                v = v + dv
        return -999.
 
print(meteor(100000))