生活在南非克鲁格国家公园的黑斑羚种群 x(t)可以用如下方程来建模。 dx/dt=(r-bx sin at)x 其中 r，b 和 a 是常数，输入它们的值和 x 的初值，计算两年时间内每个月的黑斑羚种群数量， 并绘制变化曲线。

根据题目，我们可以得到微分方程：

dx/dt = (r - bx sin(at))x

其中 r，b 和 a 是已知常数，而 x 则是未知函数，表示黑斑羚种群数量。我们需要求解这个微分方程来得到黑斑羚种群数量随时间的变化。

由于题目要求我们计算两年时间内每个月的黑斑羚种群数量，因此我们可以将时间 t 以月为单位来计算，也就是说，每个月相当于时间增加了 1。

为了求解微分方程，我们可以使用欧拉法（Euler's method）进行数值解法。欧拉法是一种简单的数值积分方法，其基本思想是通过根据微分方程的斜率来估计未来的值。具体过程如下：

1. 根据题目所给的参数和初值，设定初始条件 x(0) = x0。

2. 对于每个月 i，计算出黑斑羚种群数量 x(i)。根据欧拉法，我们有：

x(i+1) = x(i) + dx/dt * dt

其中 dx/dt 可以直接根据微分方程计算，即：

dx/dt = (r - bx sin(at))x

而 dt = 1，因为每个月时间增加了 1。

3. 重复步骤 2，直到计算出两年时间内每个月的黑斑羚种群数量。

4. 将计算出的黑斑羚种群数量绘制成变化曲线。

下面是 Python 代码实现：

import matplotlib.pyplot as plt

# 定义微分方程的参数
r = 0.2
b = 0.01
a = 0.1

# 定义初始条件
x0 = 1000

# 定义时间范围（以月为单位）
t0 = 0
t1 = 24

# 定义时间步长
dt = 1

# 计算黑斑羚种群数量随时间的变化
t = [t0]
x = [x0]

for i in range(t0, t1, dt):
    dxdt = (r - b * x[-1] * sin(a * i)) * x[-1]
    x.append(x[-1] + dxdt * dt)
    t.append(i + dt)

# 绘制变化曲线
plt.plot(t, x)
plt.xlabel('Time (months)')
plt.ylabel('Population')
plt.title('Population of Impalas in Kruger National Park')
plt.show()

在运行完上述代码之后，我们将得到一条变化曲线，它展示了两年时间内每个月黑斑羚种群数量的变化。你可以根据自己的需要修改或调整代码，以便更好地理解和掌握这个问题。