Pascal нахождение детерминанты методом гауса Программа не работает(

[Студент13]

procedure INPUT;
begin
assign(f,'DET.txt');
reset(f);
writeln;
write('введите порядок матрицы: n=');
readln(n);
while not eof(f) do
begin
for i:=1 to n do
begin
for j:=1 to n do
read(f,a[i,j]);
readln(f);
end;
end;
close(f);
end;
procedure TREAT;
begin
for k:=1 to n do
begin
if a[k,k]=0 then
begin
for m:=k to n do
begin
buf:=a[k,m];
a[k,m]:=a[k+1,m];
a[k+1,m]:=buf;
end;
end;
for i:=k to n do
for j:=k to n do
if i=k then
c[k,j]:=a[k,j]/a[k,k]
else c[i,j]:=a[i,j]-a[i,k]*c[k,j];
if k<>n then
for i1:=k+1 to n do for j1:=k+1 to n do a[i1,j1]:=c[i1,j1];
end;
end;
procedure OUTPUT;
begin
D:=1;
for i:=1 to n do D:=D*a[i,i];
end;
begin
clrscr;
INPUT;
writeln;
TREAT;
OUTPUT;
writeln('det=',D:7:5);
readln;
end.



отдельно создаём текстовый файл... Туда записываем матрицу.... помогите пожалуйста разобраться почему неработает

[Son Of Pain]

Первая ошибка, которую вижу - убери while not eof(f) do

Дальше не читал.