import matplotlib.pyplot as plt
x = [0 for _ in range(100)]
y = [0 for _ in range(100)]
x[0] = 0
y[0] = 1
def k(i):
return y[i] - 2.0 * x[i] / y[i]
def he(t):
for i in range(t):
y[i + 1] = 1.1 * y[i] - 0.2 * x[i] / y[i]
x[i + 1] = x[i] + 0.1
y[i + 1] = y[i] + 0.1 * (k(i) + k(i + 1)) / 2.0
he(99)
for i in range(99):
plt.plot(x, y)
plt.show()
改进的欧拉公式
最新推荐文章于 2023-04-23 15:10:27 发布