Решение СЛАУ Гауссом, используя данные из файла

uses crt;
 
const
 e = 0.000001;
 
var f:text;
    a,a1: array[1..10, 1..10] of real;
    b,b1: array[1..10] of real;
    n, i, j, k: integer;
    z, r, g: real;
    nf,s:string;
begin
clrscr;
write('imya fayla: ');
readln(nf);
nf:=nf+'.txt';
assign(f,nf);
reset(f);
if IOResult <> 0 then
 begin
  writeLn('file',nf,' ne nayden');
  readln;
  exit;
 end;
read(f,n);
for i:=1 to n do
for j:=1 to n do
read(f,a[i,j]);
for i:=1 to n do
read(f,b[i]);
close(f);
writeln('rashirenaya matrica: ');
for i:=1 to n do
 begin
  for j:=1 to n do
  write(a[i,j]:6:2);
  writeln(b[i]:8:2);
 end;
readln
end.
 
 
  a1[i,j]:=a[i,j];
  for k := 1 to n do
  begin
    for j := k + 1 to n do
    begin
      r := a[j, k] / a[k, k];
      for i := k to n do
      begin
        a[j, i] := a[j, i] - r * a[k, i];
      end;
      b[j] := b[j] - r * b[k];
    end;
  end;
  for k := n downto 1 do
  begin
    r := 0;
    for j := k + 1 to n do
    begin
      g := a[k, j] * x[j];
      r := r + g;
    end;
    x[k] := (b[k] - r) / a[k, k];
  end;
  writeln('Корни системы:');
  for i := 1 to n do
    write('x[', i, ']=', x[i]:0:2, '   ');
  writeln;
  writeln('Проверка в соответствии с матрицей B: ');
 
  for i:=1 to n do
   begin 
    for j:=1 to n do  
      b1[i]:=b1[i]+a1[i,j]*x[j];     
     writeln(b1[i]:3:2,' ');
    end; 
  readln;
end.