нахождение вектора невязки

[Baburek]

помогите пожалуйста. Решаю СЛАУ методом Гаусса. Как вывести вектор невязки? (Вектор невязки: r=Ay-f)

...
var a:Tmatrix; //массивы коэффициентов системы
 x,b,r:Tmass; //массивы корней уравнений и свободных членов
 k,j,row:integer;

for k:=0 to n-1 do
    r[k]:=b[k];
for k:=0 to n-1 do
    for j:=0 to n-1 do
    begin
        r[k]:=r[k]-a[k,j]*x[k];
        tab2.Cells[2,k]:=floattostrF(r[k],ffFixed,5,3);
        end

помогите пожалуйста разобраться. Что нужно исправить?