改进欧拉算法解题

最新推荐文章于 2022-06-21 11:33:19 发布

qq_45971115

最新推荐文章于 2022-06-21 11:33:19 发布

阅读量223

点赞数

分类专栏：笔记文章标签：算法

本文链接：https://blog.csdn.net/qq_45971115/article/details/120448403

版权

笔记专栏收录该内容

11 篇文章 0 订阅

订阅专栏

代码如下

import matplotlib.pyplot as plt
y = 1
x = 0
xx = [0] * 101
yy = [0] * 101
t = 0.1
for i in range (1,101):
    y = 1.1 * y - 0.2 * x / y
    x = x + 0.1
    k1 = (y - yy[i-1])/(x - xx[i-1])
    k2 = ((y+k1*t)-yy[i-1])/((x+t)-xx[i-1])
    xx[i] = x
    yy[i] = y + 0.1 * 0.5 * ( k1 + k2 )
    print(yy[i],xx[i])
plt.plot(xx,yy)
plt.show()

运行结果如下图