Метод Гаусса с получением средне арифметического

program Gaus;
uses crt;
var k,n,i,j,k1,kn:integer;
    r:real;
    x,b,d,f:array[1..10]of real;
    a,c:array[1..10,1..10]of real;
label 1;
begin
clrscr;
write('vvedie por.sist.n=');
readln(n);
for i:=1 to n do
begin
write('b[',i,']=');
readln(b[i]);
d[i]:=b[i];
for j:=1 to n do
begin
write('a[',i,j,']=');
readln(a[i,j]);
c[i,j]:=a[i,j];
end;
end;
for k:=1 to n-1 do
begin
if a[k,k]=0 then begin
for k1:=k+1 to n do
begin
if a[k1,k]<>0 then
begin
kn:=k1;
for j:=1 to n do
begin
r:=a[k,j];
a[k,j]:=a[kn,j];
a[kn,j]:=r;
end;
r:=b[k];
b[k]:=b[kn];
b[kn]:=r;
goto 1
end;
end;
writeln('ved.element=0');
readln;
halt;
end;
1:b[k]:=b[k]/a[k,k];
for i:=k+1 to n do
b[i]:=b[i]-a[i,k]*b[k];
for j:=k+1 to n do
begin
a[k,j]:=a[k,j]/a[k,k];
for i:=k+1 to n do
a[i,j]:=a[i,j]-a[i,k]*a[k,j];
end;
end;
x[n]:=b[n]/a[n,n];
for i:=n-1 downto 1 do
begin
for j:=i+1 to n do
b[i]:=b[i]-x[j]*a[i,j];
x[i]:=b[i];
end;
for i:=1 to n do
writeln('x[',i,']=',x[i]:7:5);
for i:=1 to n do
begin
f[i]:=-d[i];
for j:=1 to n do
f[i]:=f[i]+c[i,j]*x[j];
writeln ('F[',i,']=',f[i]:7:5);
end;
readln;
end.