ex8.m

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

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

A = [4 -2 7 3 -1 8; 
  -2 5 1 1 4 7; 
  7 1 7 2 3 5; 
  3 1 2 6 5 1; 
  -1 4 3 5 3 2; 
  8 7 5 1 2 4;]; 
v = [1;0;1;0;0;1;]; 
[x lam it] = pm(A, v, tol)

tol = 1e-4; 
A = [2 -1 0 0 0; 
  -1 2 -1 0 0; 
  0 -1 2 -1 0; 
  0 0 -1 2 -1; 
  0 0 0 -1 2;]; 
v = [1;1;1;1;1;]; 
[x lam it] = pm(A, v, tol)

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

tol = 1e-4; 
A = [-1 2 1; 
    2 -4 1; 
    1 1 -6;]; 
v = [1;1;1;]; 
[x lam it] = pm(A, v, tol)
[x lam it] = pm_opposite(A, v, tol, 2)