改进欧拉公式

最新推荐文章于 2024-11-03 20:27:08 发布

柠檬上神

最新推荐文章于 2024-11-03 20:27:08 发布

阅读量477

点赞数 1

文章标签： python 欧拉公式

本文链接：https://blog.csdn.net/qq_58146842/article/details/120448270

版权

# 改进欧拉公式
import numpy as np
import matplotlib.pyplot as plt

y = 1  # 初值
x = 0
a = 1000 # 循环次数 即所求点数
x1 = np.zeros(a)
y1 = np.zeros(a)
h = 0.01 # 步长
for i in range(a):
    y0 = 1.1 * y - 0.2 * x / y
    k1 = y - 2 * x / y
    k2 = y0 - 2 * (x + h) / y0
    k = (k1 + k2) / 2
    y = y + h * k
    x = x + h
    x1[i] = x
    y1[i] = y
    print(x, y)
plt.plot(x1, y1)
plt.show()