%% condA Числа обусловленности матрицы

clc; clear all; close all;

con2 = []; % cond по 2 норме матрицы
con1 = []; % cond по 1 норме матрицы
coni = []; % cond по inf норме матрицы
con1e = []; % оценка cond по 1 норме матрицы
con1r = []; % оценка cond по 1 норме матрицы (LAPACK)
mn = 4:2:20;      %  для матрицы Гильберта
%mn = 10:50:500;   %  для других тестовых матриц
for n = mn
  n 
  % theta=1; A = gallery('dorr',n,theta); 
  % A = gallery('lehmer',n);
  % A= gallery('minij',n);
  A = hilb(n);
  con2 = [con2; cond(A,2)]; 
  con1 = [con1; cond(A,1)]; 
  coni = [coni; cond(A,inf)];
  con1e = [con1e; condest(A)];
  con1r = [con1r; 1/rcond(A)];
end
disp('Числа обусловленности матрицы')
disp('           n     cond_1       cond_1r    cond_1(est)    cond_2      cond_inf')
disp([mn' con1 con1e con1r con2 coni])