要使用编程绘制引力图形，你可以遵循以下步骤：
确定地球模型
选择一个地球模型，该模型应包含不同半径处的密度信息。
例如，可以使用分层密度模型，其中每一层都有不同的密度值。
计算引力
使用万有引力公式计算每个质点受到的引力。公式为：
\[
F = G \frac{m_1 m_2}{r^2}
\]
其中，\（F\） 是引力大小，\（G\） 是万有引力常数，\（m_1\） 和 \（m_2\） 是两个质点的质量，\（r\） 是它们之间的距离。
绘制图形
根据密度分层，先按最小密度绘制大圆图。
逐层减去上一层的密度，并绘制新的图形，直到地核。
将所有层的图形叠加起来，得到最终的引力图形。
使用编程语言实现
选择一种编程语言，如Python，来实现上述计算和绘图过程。
使用Python的库，如`numpy`和`matplotlib`，来进行数值计算和图形绘制。
```python
import numpy as np
import matplotlib.pyplot as plt
地球模型参数
radii = np.array（[2.75, 3.4, 5.35, 10.7, 12.3, 12.5]）
densities = np.array（[0.99482, 0.846178, 0.544812, 0.262282, 0.194318, 0.1]）
计算引力
G = 6.67430e-11 万有引力常数，单位：m^3 kg^-1 s^-2
def calculate_gravity（r, m1, m2）:
return G * m1 * m2 / r2
绘制引力图形
fig, ax = plt.subplots（）
for i in range（len（radii） - 1）:
r1, r2 = radii[i], radii[i + 1]
density_diff = densities[i + 1] - densities[i]
for m2 in radii:
F = calculate_gravity（r2, 1, m2）
F_divided_by_r2 = F / r2
ax.plot（[r1, r2], [F_divided_by_r2, 0], label=f'r = {r2} km'）
ax.set_xlabel（'Radius （km）'）
ax.set_ylabel（'Gravitational Force （N/kg）'）
ax.legend（）
plt.show（）
```
这个代码示例计算并绘制了地球不同半径处的引力值。你可以根据需要调整地球模型参数和绘图设置。