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

[Medvedko0116]

Здравствуйте, написала программу, но почему то ответ не сходится!!Проверьте плиз))
Вот исходная сис.
x1-0.2x2+0.1x3=0.4
-0.1x1+x2-0.1x3=0.8
-0.3x1+0.2x2-x3=0.2

Вот листинг программы
rogram Prac6;
const
n=3;
var
a: array [1..n,1..n] of real;
b: array [1..n] of real;
x1: array [1..n] of real;
x2: array [1..n] of real;
i,j,k,kol:Integer;
s:real;
p:byte;
m,v,c:Real;
begin
{ Write('vvodim el.-ti matrici sistemi');
For i:=1 to n do
For j:=1 to n do
readLn( a[i,j]);}
a[1,1]:=1;
a[1,2]:=0.2;
a[1,3]:=0.1;
a[2,1]:=-0.1;
a[2,2]:=1;
a[2,3]:=-0.1;
a[3,1]:=-0.3;
a[3,2]:=0.2;
a[3,3]:=-1;
{ WriteLn('vvedite svobodnie cheleni');
For i:=1 to n do
ReadLn(b[i]); }
b[1]:=0.4;
b[2]:=0.8;
b[3]:=0.2;

{ For i:=1 to n do
ReadLn(x1[i]);}
x1[1]:=0.41;
x1[2]:=0.81;
WriteLn('Pryamoi xod');
for k:=1 to n-1 do
begin
for i:=k+1 to n do
begin
for j:=1 to n do
begin
c:=-(a[i,j])/a[k,k];
a[i,j]:=(a[k,j]*c)+a[i,j];
end;
b[i]:=(b[k]*c)+b[i];
end;
end;
For i:=1 to n do
begin
for j:=1 to n do
begin
Write(a[i,j]:6:1{x1[i]:3:2});
end;
WriteLn;
end;
WriteLn('obratnyi xod');
x1[n]:=b[n]/a[n,n];
For i:=n-1 downto 1 do
begin
s:=0;
For j:=i+1 to n do
begin
s:=s+a[i,j]*x1[i];
x1[i]:=(b[i]*s)/a[i,i];
end;
end;
For i:=1 to n do
begin
Write(x1[i]:7:2);
end;
WriteLn;
ReadLn;
end.