Op code niga wott

from math import
masa = 75
gravedad = 9.8
alpha = 5500.0
beta = 0.0069
alturaInicial = 36068
velocidadInicial = 0
re = 6371009.0 # Radio de la tierra
nu = 0.4



def c(y, t):
    return (nu*t + (3-t)*beta*(e**(-y/alpha)))/3.0

def g(y):
    return gravedad * ((re/(re+y))**2)

def uNMasUno(un, vn, h):
    return un + h*vn

def vNMasUno(un, vn, h, t):
    return (1-(1-4*h*c(un,t)*(vn-h*gravedad))**0.5)/(2*(h*c(un,t)))


h = 0.25
t = 0
altura = alturaInicial
alturaEstrella = 0
velocidad = velocidadInicial
duracion = 0
vmax = 0
while altura > 0:
    duracion += h
    if altura <= 1500 and t < 3:
        t += h
    alturaEstrella = round(uNMasUno(altura, velocidad, h),5)
    velocidad = round(vNMasUno(alturaEstrella, velocidad, h, t), 5)
    altura = round(uNMasUno(altura, velocidad, h),5)
    pasoCritico = abs(1/(c(altura,t)*velocidad))
#    if h >= pasoCritico:
#        break
    if velocidad < vmax:
        vmax = velocidad
    print velocidad/1000*-3600, duracion
print duracion
print vmax/1000*3600
variable = raw_input('input something!: ')