有解析解

这个方程换一下元就很容易（是不是我想简单了）？
GM/(R+z)^2 = d^2z/dt^2
 
设 y = R+z  则
d^2y/dt^2 =  GM/y^2
 
dy/dt = -GM/y + C1
y = -GM*ln(t) + C1*t + C2

因此
z = y-R
  = -GM*ln(t) + C1*t + C2 - R
C1 和 C2 都是常数