Решение СЛАУ методом Гауса на Паскале

[Klik_1602]

Помогите решить задачу на Паскале, вот текст:
]Выполните «обратный ход» в решении системы линейных алгебраических уравнений методом Гаусса. Вводите только элементы матрицы системы, стоящие выше главной диагонали, и вектор свободных членов.

Если не трудно, можете объяснить или кинуть ссылку, где объяснено, что такое обратный и прямой ход?

заранее спасибо

[k0tone]

Обратный ход, заключается в решении треугольной системы. Из последнего уравнения находим x[m]. По найденному x[m] из (m-1) уравнения находим x[m-1]. Затем по x[m-1] и x[m] из (m-2) уравнения находим x[m-2]. Процесс продолжаем, пока не найдем x[1] из первого уравнения.

Источник:
http://www.mathelp.net/MA5.htm

[Klik_1602]

подскажите плиз, где ошибка в коде?
program L7_23;
uses
crt;
var
f:text;
n,i,j:integer;
h:real;
x:array[1..10] of real;
b:array[1..10] of real;
a:array[1..10,1..10] of real;
procedure vvodmatr(n,i,j:integer; var a:array[1..10,1..10] of real);
begin
assign(f, 'L7_23.txt');
reset(f);
for i:=1 to n do
begin
for j:=1 to n do
read (f, a[i,j]);
readln(f);
end;
close(f);
end;
procedure vivodmatr(n,i,j:integer; var a:array[1..10,1..10] of real);
begin
for i:=1 to n do
begin
for j:=1 to n do
write (a[i,j]:1);
writeln;
end;
end;
procedure diagonal(n,i,j:integer; var a:array[1..10,1..10] of real);
begin
for j:=1 to n do
begin
for i:=1 to n do
if (j<i) then
a[i,j]:=0
else
if (i=j) then
a[i,j]:=1;
end;
writeln;
writeln;
writeln;
for i:=1 to n do
begin
for j:=1 to n do
write (a[i,j]:1);
writeln;
end;
end;
procedure stroka(n,i:integer; var b:array[1..10] of real);
begin
reset (f);
for i:=1 to n do
read (f, b[i]);
readln(f);
close(f);
writeln;
writeln;
writeln;
for i:=1 to n do
write (b[i]:1);
writeln;
end;
procedure korni(n,i,j:integer; h:real;
x:array[1..10] of real;
b:array[1..10] of real;
a:array[1..10,1..10] of real);
begin
for i:=n 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;
writeln;
writeln;
writeln('Vivod korneilя');
for i:=1 to n do writeln('x(',i,')=',x[i]:10:5);
end;
begin
clrscr;
write('Vvedite chislo uravnenij: ');
read(n);
vvodmatr(n,i,j,a);
vivodmatr(n,i,j,a);
diagonal(n,i,j,a);
stroka(n,i,b);
korni(n,i,j,h,x,b,a);
end.

[Klik_1602]

если не трудно, подскажите как исправить ошибку