Метод Зейделя: Ошибка времени выполнения, в экземпляре объекта не задана ссылка на объект

const Eps=0.001;
 
Type
  TMatrix = array of array of real;
  TVector = array of real;
  
Procedure FillA (var A: TMatrix; var n: integer);
var
   i,j: integer;
begin
     for i:=1 to n do
         for j:=1 to n do
             readln (A[i,j]);
end;
 
procedure PrintA (A: Tmatrix;var n: integer);
var
   i,j: integer;
begin
     for i:=0 to n do
     begin
          for j:=0 to n do
          Write (A[i,j]:4,' ');
     writeln;
     end;
     writeln;
end;
 
Procedure FillB(var b: TVector; var n: integer);
var
  i: integer;
begin
    for i:=0 to n do
    readln(b[i]);
end;
 
Procedure PrintB (b: Tvector; var n: integer);
var
  i: integer;
begin
    for i:=0 to n do
    Write(b[i]:2,'');
end;
 
Procedure PreobrazAvD1 (var A,D1: TMatrix; var n: integer);
var
  i,j: integer;
begin
    for i:=1 to n do 
    begin
       for j:=1 to n do
       begin
        if i=j then D1[i,j]:=0 else
          D1[i,j]:=(-A[i,j])/(A[i,i]);
        write(D1[i,j]:0:2,'  ');
       end;
    writeln;
    end;
end;
 
Procedure max1_sum_stroka (var D1: TMatrix; var sum, max1: real;
  var imax1,n: integer);
var
  i,j: integer;
begin
  max1:=0;
  for j:=0 to n do
  begin
    sum:=0;
    for i:=0 to n do
       sum:=sum+abs(D1[j,i]);
       writeln('Sum [',j,'] = ',sum:0:2);
       if sum>max1 then
       begin
         max1:=sum;
         imax1:=j;
       end;
  end;
end;
 
Procedure PreobrazbvC(var b,C:TVector; A:TMatrix; var n:integer);
var
   i: integer;
begin
  for i:=0 to n do
    begin
    C[i]:=b[i]/A[i,i];
    writeln(C[i]:0:2);
    end;
end;
     
Procedure PreobrazD1vD2 (var D1,D2,A: TMatrix; var n: integer);
var
  i,j: integer;
begin
    for i:=1 to n do 
    begin
       for j:=1 to n do
       begin
        if i<j then D2[i,j]:=0 else
          D2[i,j]:=(-A[i,j])/(A[i,i]);
        write(D2[i,j]:0:2,'  ');
       end;
    writeln;
    end;
end;
 
 
Procedure max2_sum_stroka (var D2: TMatrix; var sum, max2: real;
  var imax2,n: integer);
var
  i,j: integer;
begin
  max2:=0;
  for j:=0 to n do
  begin
    sum:=0;
    for i:=0 to n do
       sum:=sum+abs(D2[j,i]);
       writeln('Sum [',j,'] = ',sum:0:2);
       if sum>max2 then
       begin
         max2:=sum;
         imax2:=j;
       end;
  end;
end;
 
  
var
  i,j: integer;
  A: TMatrix;
  b: TVector;
  C: TVector;
  x0,x1,r:TVector;
  D1,D2: TMatrix;
  sum, max1,max2,E,maxR: real;
  imax1, imax2: integer;
  n:integer;
  
  
begin
  Writeln ('Enter the dimension of the matrix');
  Readln (n);
 
  Writeln ('Enter the matrix A');
  FillA(A,n);
  PrintA(A,n);
  
  Writeln ('Enter the vector b');
  FillB(b,n);
  PrintB(b,n);
  
  Writeln ('Transform matrix A into the matrix D1');
  PreobrazAvD1 (A,D1,n);
 
  max1_sum_stroka (D1,sum, max1,imax1,n);
  writeln('Norm = ',max1:0:2,'<1, str = ',imax1);
 
  Writeln ('Transform vector b in vector C');
  PreobrazbvC(b,C,A,n);
 
  Writeln ('Transform matrix D1 to D2');
  PreobrazD1vD2 (D1,D2,A,n);
  
  max2_sum_stroka (D2,sum, max2,imax2,n);
  writeln('Norm = ',max2:0:2,'<1, str = ',imax2);
 
  writeln;
 
  E:=((1-max1)/max2)*Eps;
  writeln('E = ',E:0:4);
 
//---------x0-------------
  for i:=0 to n do
  x0[i]:=C[i];
  Writeln;
//---------x1-------------
  repeat
  for i:=1 to n do
  begin
    x1[i]:=c[i];
    for j:=1 to n do
    x1[i]:=D2[i,j]*x0[j]+x1[i];
    write('x1[',i,']= ',x1[i]:0:4,'  ');
    end;
    writeln;
for i:=1 to n do
    begin
    maxR:=0;
    r[i]:=abs(x1[i]-x0[i]);
    if r[i]>maxR then
    maxR:=r[i];
    end;
if maxR<E then
writeln('Otvet: x1=',x1[1]:0:4,' x2=',x1[2]:0:4,' x3=',x1[3]:0:4) else
    begin
    for i:=1 to n do
    x0[i]:=x1[i];
    end;
until maxR<E;
 
end.

Procedure FillA (var A: TMatrix; var n: integer);
var
   i,j: integer;
begin
     SetLength(A, n, n);
     for i:=0 to n-1 do // Обрати внимание на индексы !!!
         for j:=0 to n-1 do
             readln (A[i,j]);
end;