def p(self, b): return mars_density(b.radius) * b.speed2 / 2
def q(self, b): return self.k * sqrt(mars_density(b.radius) / b.r_nose) * b.speed ** 3