loist91的博客

不确定证明过程是否严格，尽力写全面了。

import numpy as npimport matplotlib.pyplot as pltx1 = 0x2 = 0y1 = 1y2 = 1h = 0.1 # x(n+1) - x(n)X1 = np.zeros(100)Y1= np.zeros(100)X2 = np.zeros(100)Y2 = np.zeros(100)def f(x, y): return y - 2 * (x / y)for i in range(0,100): y1 = y1