%% Stab: Точность решения СЛАУ 

clc; clear all; close all;

con2 = []; % cond по 2 норме матрицы
coni = []; % cond по inf норме матрицы
err2 = []; % относительная ошибка СЛАУ по 2 норме
erri = []; % относительная ошибка СЛАУ по inf норме
mn = 4:2:20;       %  для матрицы Гильберта
%mn = 10:50:500;   %  для других тестовых матриц
for n = mn
  n   
  % theta=0.01; A = gallery('dorr',n,theta); 
  % A = gallery('lehmer',n);  
  % A= gallery('minij',n);
  A = hilb(n);
  
  x = rand(n,1); % точное решение СЛАУ
  b = A*x;       % правая часть СЛАУ
  y = A\b;       % решение методом Холецкого 
  
  con2 = [con2; cond(A,2)]; 
  coni = [coni; cond(A,inf)];
  
  err2 = [err2; norm(x-y,2)/norm(x,2)]; 
  erri = [erri; norm(x-y,inf)/norm(x,inf)];
end
disp('Точность решения СЛАУ')
disp('           n      err_2        cond_2      err_inf      cond_inf')
disp([mn' err2 con2 erri coni])