代码
from matplotlib import pyplot as plt
import numpy as np
h = 0.01
yy = []
xx = []
yy.append(1)
xx.append(0)
def y1(m,n):
return m - 2*n/m
for i in range(100):
k1 = y1(yy[i],xx[i])
y = yy[i] + hy1(yy[i],xx[i])
xx.append(xx[i]+h)
k2 = y1(y,xx[i+1])
yy.append(yy[i]+h1.0/2.0*(k2+k1))
print(“y=”,yy[i],“x=%.1f”%xx[i])
plt.plot(xx,yy)
plt.show
截图