4楼: Originally posted by lurencyj at 2012-04-20 16:20:50:
上面发的帖子格式没对，呵呵，故重复一下：


“赋值A*，代入A=f（A*），获得A，若A=A*，则结束，不等则重新赋值A*=A*-x，反复循环，直到获得最终结果”

就把A看作一个数，不是一个矩阵。

ErrorA = 1E-3 ...

ErrorA = 1E-3;   % 迭代的容忍误差
A = A0;                % A0是初始数值，一定程度上随便选择的一个数
A_pre = A0 + 10*ErrorA;             % 定义上次迭代的结果的变量。

while   abs(A - A_pre) > ErrorA
        A_pre = A;            %保存上次计算现场
        A = f(A + x);            %计算新的A数值
end