改进的欧拉公式

import matplotlib.pyplot as plt


x = [0 for _ in range(100)]
y = [0 for _ in range(100)]
x[0] = 0
y[0] = 1


def k(i):
    return y[i] - 2.0 * x[i] / y[i]


def he(t):
    for i in range(t):
        y[i + 1] = 1.1 * y[i] - 0.2 * x[i] / y[i]

        x[i + 1] = x[i] + 0.1

        y[i + 1] = y[i] + 0.1 * (k(i) + k(i + 1)) / 2.0


he(99)
for i in range(99):
    plt.plot(x, y)
plt.show()