改进欧拉公式

代码实现

from matplotlib import pyplot as plt


def y(x, y):
    y = 1.1 * y - (0.2 * x) / y
    return float("%.2f" % y)


def collection(x0, y0, h):
    x1 = x0 + h
    yp = y0 + h * y(x0, y0)
    yc = y0 + h * y(x1, yp)
    y1 = (1 / 2) * (yc + yp)
    return float("%.2f" % x1), float("%.2f" % y1)


def xy_list(x, y):
    x_list.append(x)
    y_list.append(y)


def draw(x_list, y_list):
    plt.figure(figsize=(10, 10), dpi=100)
    plt.scatter(x_list, y_list)
    plt.show()


if __name__ == '__main__':
    x_list = []
    y_list = []
    result = collection(0, 1, h=0.1)
    print('x1:', result[0], 'y1:', result[1])
    xy_list(result[0], result[1])
    for i in range(100):
        result = collection(result[0], result[1], 0.1)
        xy_list(result[0], result[1])
        print('x{}:'.format(i + 2), result[0], 'y{}:'.format(i + 2), result[1])
    draw(x_list, y_list)

运行截图