Verhulst模型的Python代码

Verhulst模型是一种常见的人口增长模型，可以用以下公式表示：

dN/dt = rN(1 - N/K)

其中，N是种群数量，t是时间，r是种群的增长率，K是种群的环境容量。

以下是Verhulst模型的Python代码：

import matplotlib.pyplot as plt
import numpy as np

# 设置参数
r = 0.5
K = 100
N0 = 10
dt = 0.1
t = np.arange(0, 50, dt)

# 定义函数
def verhulst(N, r, K):
    return r * N * (1 - N / K)

# 初始化种群数量变量
N = np.zeros(len(t))
N[0] = N0

# 循环计算种群数量
for i in range(1, len(t)):
    N[i] = N[i-1] + dt * verhulst(N[i-1], r, K)

# 绘制图像
plt.plot(t, N)
plt.xlabel('Time')
plt.ylabel('Population')
plt.title('Verhulst Model')
plt.show()

运行该代码可以得到一个人口数量随时间变化的图像。你可以通过修改参数来观察不同条件下的种群数量变化情况。