给定以下公式,我需要绘制粒子在时间 t 的位置:s(t) = -0.5*g(s)*t^2+v0*t
,其中g(s) = G*M/(R+s(t))^2
(G
,M
, 和R
是常数,s
是一个值,而不是函数s(t)
)。粒子被垂直发射,我想每秒打印一次它的当前位置,直到它撞到地面。但是我无法弄清楚如何在定义一个函数之前不使用另一个函数来定义它。到目前为止,这是我的代码:
G = 6.6742*10^(-11)
M = 5.9736*10^24
R = 6371000
s0 = 0
v0 = 300
t = 0
dt = 0.005
def g(s):
def s(t):
s(t) = -0.5*g(s)*t^2+v0*t
g(s) = G*M/(R+s(t))^2
def v(t):
v(t) = v(t-dt)-g(s(t-dt))*dt
while s(t) >= 0:
s(t) = s(t-dt)+v(t)*dt
t = t+dt
if t == int(t):
print s(t)
当我运行该函数时,它说它无法分配函数调用。