Метод Гаусса Опции

[moto]

Помогите плиз реализовать метод Гаусса  в Turbo pascal как можно более простым способом. Алгоритм решения ясен, вот только возникла проблема в приведении расширенной матрицы к диагональному виду, голову сломала  ??? помогите плиз!!!

[poison]

Программа выглядит так:

Код

Program Gaus;
       const n=4;
       var  A: array [1..n,1..n] of Real;
            B: array [1..n] of Real;
            x: array [1..n] of Real;
            i, j, k:           Integer;
            L, G:              Real;
        Begin
             Writeln;
             for i:=1 to n do
                      x[i]:=0;
             for i:=1 to n do
                 begin
                      Write ('B[',i,']=');
                      Read  (B[i]);
                 end;
             for i:=1 to n do
                 for j:=1 to n do
                     read (A[i,j]);
             for k:=1 to (n-1) do
                 for i:=(k+1) to n do
                     begin
                          L:=A[i,k]/A[k,k];
                          B[i]:=B[i]-L*B[k];
                          for j:=k to n do
                              A[i,j]:=A[i,j]-L*A[k,j];
                     end;
             for i:=1 to n do
                 begin
                      writeln;
                      for j:=1 to n do
                          write (A[i,j]);
                 end;
             x[n]:=B[n]/A[n,n];
             for i:=(n-1) downto 1 do
                 begin
                      for k:=(i+1) to n do
                          G:=G+A[i,k]*x[k];
                      x[i]:=1/A[i,i]*(B[i]-G);
                      G:=0;
                 end;
             writeln;
             writeln;
             for i:=1 to n do
                 begin
                      writeln;
                      write ('X[',i,']=',x[i]);
                 end;
             end.

При выводе трехугольной матрицы случаются ошибки, но ответ всегда правильный.