Ax = y
x = y / A
Let B = I - A
x = y / ( I - B)
1 / (I - B) = I + B + B**2 + B**3 + B**4 + ...
你需要调 A such that B 的 eigenvalues 的绝对值都小于 1,则上面的展式是收敛的。
这是well known的算法,可以加上regularization以增加稳定性。