本文通过一个具体的数值方法示例，展示了如何使用Python中的NumPy和Matplotlib库来解决一个微分方程问题并可视化其结果。代码中采用了一种类似于欧拉方法的数值积分方案来逼近解。
import numpy as np
import matplotlib.pyplot as plt
x = 0
y = 1
step = 0.01
xx = np.zeros(1000)
yy = np.zeros(1000)
for i in range(1000):
k0 = y - 2 * x / y
y = 1.1 * y - 0.2 * x / y
x = x + step
k1 = y - 2 * x / y
y = y + 0.5 * step * (k0 + k1)
yy[i] = y
xx[i] = x
plt.plot(xx,yy)
plt.show()
[外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-NlL20RIF-1632447676282)(output_0_0.png)]
扫一扫
被折叠的 条评论
为什么被折叠?
到【灌水乐园】发言
