p228_pm_opposite.m

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

a = [-12 3 3; 3 1 -2; 3 -2 7];
x0 = [1 1 1]';
tol = 0.0001;
p = -13;
[x, lam, it] = pm_opposite(a, x0, tol, p)