метод гаусса в паскале

[zenid]

помогите исправить программу., она не правильно считает и не работает, выдает ошибку деление на ноль или что то в этом роде.
код программы:

program Gayss;
const n=3;
var i,j,k:integer;
c,s:real;
a:array [1..n,1..n] of real;
b:array [1..n] of real;
x:array [1..n] of real;
begin
for i:=1 to n do
for j:=1 to n do read (a[i,j]);
writeln ('vvedite koaf cvobodnih chlenov');
for i:=1 to n do read (b[i]);
k:=1;
repeat i:=k;
repeat c:=a[i+1,k]/a[k,k];
j:=1;
repeat a[i+1,j]:=a[i+1,j]-c*a[k,j];
j:=j+1;
until j>=n;
b[i+1]:=b[i+1]-c*b[k];
i:=i+1
until i>n-1;
k:=k+1;until k>=(n-1);
x[n]:=b[n]/a[n,n];
i:=n-1;
repeat s:=0;
j:=i+1;
repeat s:=s+x[j]*a[i,j];
j:=j+1 until j>n;
x[i]:=(b[i]-s)/a[i,i];
i:=i-1;
until i<0;
for i:=1 to n do writeln ('x[',i,']=',x[i]:9:7);
readln;
readln; 
end.

[zenid]

ребят, помогите пожалуйста

[McFrey]

Вроде так, но вручную не проверял

program Gayss;

const
  n = 3;

var
  i, j, k: integer;
  h: real;
  a: array [1..n, 1..n] of real;
  b: array [1..n] of real;
  x: array [1..n] of real;

begin
  for i := 1 to n do
    for j := 1 to n do read(a[i, j]);
  writeln('vvedite koaf cvobodnih chlenov');
  for i := 1 to n do read(b[i]);
  
  for i := 1 to n do
  begin
    for j := i + 1 to n do
    begin
      a[j, i] := -a[j, i] / a[i, i];
      for k := i + 1 to n do a[j, k] := a[j, k] + a[j, i] * a[i, k];
      b[j] := b[j] + a[j, i] * b[i];
    end;
  end;
  x[n] := b[n] / a[n, n];
  for i := n - 1 downto 1 do
  begin
    h := b[i];
    for j := I + 1 to n do
      h := h - x[j] * a[i, j];
    x[i] := h / a[i, i];
  end;
  writeln('Решение системы:');
  for i := 1 to n do
    writeln('X(', i, ')=', x[i]:3:5);
  readln;
end.

[zenid]

спасибо большое