Вещественное переполнение (Вычисление корней уравнения методом Гаусса)

function Zeid(s: integer; a: Tmas1; b: Tmas2): Tmas2; //Метод Зейделя
const
        ee = 0.001; //Погрешность метода
        m = 1000;  //Максимальное число итераций
var
        i,j,k: integer;
        tmp: real;
        bool: boolean;
        l: integer;
        x: Tmas2;
        dd: real;
begin
        //Устранение нулевых диагональных элементов
        repeat
                bool:= true;
                for k:= 0 to s-1 do
                begin
                        if a[k,k]=0 then
                        begin
                                bool:= false;
                                for i:= 0 to s-1 do
                                begin
                                        if a[i,i]<>0 then
                                                j:= i;
                                end;
 
 
                                for i:= 0 to s-1 do
                                begin
                                        tmp:= a[j,i];
                                        a[j,i]:= a[k,i];
                                        a[k,i]:= tmp;
                                end;
                                tmp:= b[j];
                                b[j]:= b[k];
                                b[k]:= tmp;
                        end;
                end;
        until bool;
 
        for i:= 0 to s-1 do
                x[i]:= 0;
 
        l:= 0;
        k:= 1;
        repeat
                dd:= 0;
                for i:= 0 to s-1 do
                begin
                        tmp:= 0;
                        for j:= 0 to i-1 do
                        begin
                                tmp:= tmp + a[i,j]*x[j];   //Ссылается сюда
                        end;
                        for j:= i+1 to s-1 do
                        begin
                                tmp:= tmp + a[i,j]*x[j];
                        end;
                        tmp:= (b[i]-tmp)/a[i,i];
                        if Abs(tmp - x[i]) > dd then
                        begin
                                dd:= Abs(tmp - x[i]);
                        end;
                        x[i]:= tmp;
                end;
                if dd < ee then
                begin
                        l:= 1;
                end
                else
                begin
                        if k=m then
                        begin
                                l:= 2;
                        end
                        else
                        begin
                                k:= k+1;
                        end;
                end;
        until l<>0;
 
        if l=1 then
        begin
                Zeid:= x;
        end
        else
        begin
                for i:= 0 to s-1 do
                begin
                        x[i]:= 0;
                end;
                Zeid:= x;
        end;
end;