编写一个能近似计算地球质量的小程序。

编写一个近似计算地球质量的小程序，我们可以利用万有引力定律，其中地球的质量（\(M_E\)）可以由月球绕地球运动的轨道数据估算。假设已知月球到地球的平均距离（\(r\)），月球公转周期（\(T\)），以及地球的重力加速度（\(g\)）。地球的半径（\(R_E\)）通常是一个已知的固定值。

根据开普勒第三定律，月球的公转周期 \( T \) 和它的轨道半长轴 \( a \)（实际上等于 \( r + R_E \)）之间的关系可以用圆周率（\(\pi\)）表示：
\[ T^2 = \frac{4 \pi^2}{GM_E}a^3 \]

其中 \( G \) 是万有引力常数。简化后，我们可以得到一个近似的方程：
\[ M_E = \frac{4 \pi^2 a^3}{GT^2} \]

下面是用Python编写的一个简单版本：

import math

# 定义常量
G = 6.674 * 10**(-11)  # 万有引力常数 m^3 kg^-1 s^-2
R_E = 6371 * 10**3     # 地球半径 m
r = 384400 * 10**3      # 月球平均距离 m (近似)
T_days = 29.53           # 月球公转周期天数

# 将时间转换为秒
T_seconds = T_days * 24 * 3600

# 计算地球质量
M_E_approx = (4 * math.pi**2 * (r + R_E)**3) / (G * T_seconds**2)

print(f"地球的质量大约为：{M_E_approx:.2f} 吨")

注意，这个计算非常简略，并没有考虑月球自转、地球扁平度等因素的影响，实际地球质量的精确值会通过更复杂的物理模型来确定。