p238_1_1.m

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

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

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