function sim = simulationsrRichardson(a1,b1,step1,a2,b2,step2,file)
% Funksjonen skriver de forskjellige matrisene V til en fil, der V
% inneholder spektralradiusene til iterasjonsmatrisen i Richardson med en
% bestemt verdi av µ for hver V. a1 og b1 er start- og sluttverdi på
% intervallet for N, mens step1 er steglengden. a2 og b2 er start og
% sluttverdi på intervallet for µ, mens step2 er steglengden på dette
% intervallet. file er en streng som er navnet på filen informasjonen
% lagres i
e = a2:step2:b2;
f = length(e);
fid = fopen(file,'w');
c = a1:step1:b1;
d = length(c);
V = zeros(d,f);
for i = 1:f
    V(:,i) = spectralradiusRichardson(a1,b1,step1,e(i));
end % for
    dlmwrite(file,V,'precision','%.40f');
sim = fclose(fid);
end % function