e361.m

来自「matlab算法集 matlab算法集」· M 代码 · 共 19 行

%-----------------------------------------------------------------------
% Example 3.6.1: QR Method
%-----------------------------------------------------------------------
   clc                              
   clear                            
   n = 4;
   m = 1000;
   tol = 1.e-6;
   A = [-3.9 -0.6 -3.9 0.4; 1 0 0 0; 0 1 0 0; 0 0 1 0];
   c = zeros(n,1);

   fprintf ('Example 3.6.1: QR Method\n');
   show ('A',A)
   [X,C] = eig (A);
   c = diag(C);
   show ('eigenvalues',c);
%-----------------------------------------------------------------------