Проверьте метод Гаусса

[RUSSTUDENT]

Написать программу решения по методу Гаусса системы линейных уравнений
4х1+0,24х2-0,08х3=8
0,09х1+3х2-0,15х3=9
0,04х1-0,08х2+4х3=20

Правильно ли прога написана?

Прога:

uses crt;
const n=3;
a:array[1..n,1..n] of real=((4,0.24,-0.08),
(0.09,3,-0.15),
(0.04,-0.08,4));
b:array[1..n] of real=(8,9,20);
var x:array[1..3] of real;
i,j,k:integer;z,r,g:real;
begin
clrscr;
{приведение матрицы к треугольному виду}
for k:=1 to n do
begin
for j:=k+1 to n do
begin
r:=a[j,k]/a[k,k];
for i:=k to n do
a[j,i]:=a[j,i]-r*a[k,i];
b[j]:=b[j]-r*b[k];
end;
end;
{вычисление корней}
for k:=n downto 1 do
begin
r:=0;
for j:=k+1 to n do
begin
g:=a[k,j]*x[j];
r:=r+g;
end;
x[k]:=(b[k]-r)/a[k,k];
end;
writeln('Корни системы:');
for i:=1 to n do
write('x[',i,']=',x[i]:0:2,' ');
readln;
end.

[whatever]

Ну все правильно, только массивы из const в type перенес, так как по ходу программы они изменяются => должны быть объявлены не как константы, а как переменные.

А так ответы дает верные:

[whatever]

Во избежание вопросов в ЛС:
Вторая прога на картинке - моя лабораторная по численым методам.

Формат входных данных:
matrix.txt

4 0.24 -0.08
0.09 3 -0.15
0.04 -0.08 4

vector.txt

8
9
20

Как правило считает верно, но за 100% работоспособность не ручаюсь. В общем пользуйтесь на свой страх и риск.