Выход за границы диапазона изменения индекса

Program kl_analiz;
const
n=2;
m=3;
var
k,  i, j, r, t, c: integer;
g,s, min:real;
s1,q :string;
x: array [1..n,1..m] of integer;
y: array [1..n,1..n] of real;
z: array [1..n] of integer;
d: array [1..n-1] of real;
f, g1: text;
 begin
Assign(f,'in.txt');Reset(f);
Assign(g1,'out.txt');Append(g1);
For i:= 1 to n do
 for j:=1 to m do
Read(f, x[i,j]);
 
writeln (g1, ' Исходная матрица: ');
for i:=1 to n do
begin
for j:=1 to m do
write(g1, x[i,j]);
writeln;
end;
 
for i:=1 to n do     //формирование вектора номеров объектов
z[i]:=i;
 
 
for i:=1 to n-1 do   //  построение матрицы сходства
begin
  for k:=i+1 to n do
  begin
  s:=0;
   for j:=1 to m do
   s:=sqr(x[i,j]-x[k,j])+s;
   s:=1/(1+sqrt(s));
   y[i,k]:=s;
   y[k,i]:=y[i,k];
  end;
end;
 
r:=0;          //диагонализация матрицы сходства, построение графа минимальной длины, построение дельта-         последовательности
min:=-1;
for i:=1 to n-1 do
begin
 for j:=i+1 to n do
 if min<y[i,j] then
 begin
 min:=y[i,j];
 r:=j;
 end;
end;
 
               // перестановка местами найденного объекта с первым в матрице сходства и в векторе номеров
for i:=1 to n-1 do
begin
min:=y[1,i+1];
 for j:=1 to i do
  for k:=i+1 to n do
  if  min>y[j,k] then
  begin
  min:=y[j,k];
  r:=k;
  end;
 end;
if r>i+1 then
begin       // меняем столбцы
for k:=1 to n do
begin
g:=y[k,i+1];
y[k,i+1]:=y[k,r];
y[k,r]:=g;
end;
 
begin       // меняем строки
for k:=1 to n do
begin
g:=y[i+1,k];
y[i+1,k]:=y[r,k];
y[r,k]:=g;
end;
 
t:=z[i+1];
z[i+1]:=z[r];
z[r]:=t;
d[i+1]:=min; //дельта-последовательность
end;
 
if r<(i+1) then
d[i+1]:=min;
end;
 
for j:=1 to k*r-1 do    // выделение классов
begin
min:=-1;
 for i:=2 to n do
 if min<d[i] then
 begin
 min:=d[i];
 t:=i;
 end;
d[t]:= -d[t];
end;
 
c:=0;        //вывод результатов
d[n+1]:=-1;
for i:=1 to n do
begin
str(z[i],q);
s1:=s1+q+' ';
if d[i+1]<0 then
begin
c:=c+1;
str(c,q);
s1:=q+'-й  класс:' +s1;
 writeln (g1, s1);
  end;
end;
Close(f);Close(g1);
end.