% test skripta za tridijagonalizaciju simetricne matrice
% Numericka matematika (Marko Hajba)
n = 200; %red matrice A

%neka je matrica A matrica slučajnih brojeva reda n
A = rand(n, n);

% simetriziramo A
A = 0.5 * ( A + A' );

%pozivamo skriptu za tridijagonalizaciju
[H, U, sim] = TridijagForm(A);

if sim == false break
  end
  
%graficki prikaz nenul elemenata matrice (zelimo tridijagonalnu)
spy(H)

%izračunajmo koliko je dobiveno rješenje iz algoritma daleko od matrice A
norma = norm(A -  U * H * U')/norm(A)

%broj nenul elementata je n + 2* (n-1)
%za dimenziju 200 je to 598