p238_2.m

来自「数值分析, 同济大学教材<<现代数值数学和计算>>数值分」· M 代码 · 共 6 行

A = [2 1 3 4; 1 -3 1 5; 3 1 6 -2; 4 5 -2 -1];
x0 = [1 1 1 1]';
tol = 1e-5;
[x, lam, it] = pm(A, x0, tol)