function sim = spectralradiusHilbert(a,b,step)
% Funksjonen regner ut spektralradiusen til alle iterasjonsmatrisene for
% alle metodene i oppgave 3. Denne informasjonen skrives til en matrisen V,
% som deretter skrives til en fil, som lagres. a er startverdien, mens b er
% sluttverdien og step er steglengden på intervallet for N
format long
c = a:step:b;
d = length(c);
V = zeros(d,4);
fid = fopen('Hilbert.asc','w');
for i = 1:d
    A = sparse(Hilbert(c(i)));
    I = eye(c(i),c(i));
    D = diag(diag(A));
    Jacobi = I-(D\A);
    V(i,1) = max(abs(eig(Jacobi)));
    L = tril(A);
    GaussSeidel = I-(L\A);
    V(i,2) = max(abs(eig(GaussSeidel)));
    Richardson = I\(I-A);
    V(i,3) = max(abs(eig(Richardson)));
    V(i,4) = max(eig(abs(Hilbert(c(i)))));
end % for
dlmwrite('Hilbert.asc',V,'precision','%.40f');
sim = fclose(fid);
end % function