欧拉公式改进版

import matplotlib.pyplot as plt
import numpy as np
h=0.1
x = 0
xx=np.zeros(100)
yy=np.zeros(100)
y = 1

for i in range (100):
    k = y - 2 * x / y
    x = x + 0.1
    y2=y+0.01*(y-2*x/y)
    k2 = y2-2*x/y2
    y = y + 0.1 * 0.5 * (k + k2)

    xx[i]=x
    yy[i]=y
    print(xx[i],yy[i])
plt.plot(xx,yy, ls='-', lw=2, label='OULA', color='blue')
plt.legend()
plt.show()