function sim = spectralradiusSORHilbert(a1,b1,step1,a2,b2,step2)
% Denne funksjonen regner ut spektralradiusen til iterasjonsmatrisene til
% SOR i oppgave 3, der Hilbert er en funksjon som genererer Hilbertmatriser
% med gitte størrelser.Deretter skrives all informasjon til en bestemt fil.
% a1 er startverdien på intervallet for N og b1 er sluttverdien. step1 er
% steglengden for N. a2 og b2 er hneholdvis start- og sluttverdi for ω,
% mens step 2 er steglengden for ω

c = a1:step1:b1;
d = length(c);
e = a2:step2:b2;
f = length(e);
V = zeros(f,d);
fid = fopen('HilbertSOR.asc','w');
for w = 1:f
    for v = 1:d
        A = sparse(Hilbert(c(v)));
        B = tril(A);
        for i = 1:(c(v))
            B(i,i) = B(i,i)/e(w);
        end % for
        I = eye(c(v),c(v));
        C = B\I;
        V(w,v) = max(abs(eig(I-C*A)));
    end % for
end % for
dlmwrite('HilbertSOR.asc',V,'precision','%.40f');
sim = fclose(fid);
end % function