Ref[0][J] := Ref[0][J]+ alpha(Ref[0][J]-Matrix[I,J]);

for K := 0 to 1 do //for each reference
   begin
    Memo1.Lines.Add('Reference: '+IntToStr(K));
    for I := 0 to 9 do  //for each row
       begin
        X := 0;
        for J := 0 to 9 do  //for each column
        X := X+Power(Ref[K,J]-Matrix[I,J], 2);
        X := Sqrt(X);
         EuclDist[K,I]:= X;
        Memo1.Lines.Add( Format('%0:20.18f',[X] ));
       end;
      Memo1.Lines.Add('');
   end;
 //end;
      for I := 0 to 9 do
      begin
         if ((EuclDist[0][I]) < (EuclDist[1][I]))or ((EuclDist[0][I]) = (EuclDist[1][I]))   then
         begin
            target:=0 ;
            for j:= 0 to 9 do
            if target=(matrix[i][j]) then
              begin
               Ref[0][J] := Ref[0][J]+ alpha (Ref[0][J]-Matrix[I,J]);
              end
            else
              begin
              Ref[0][J] := Ref[0][J]- alpha*(Ref[0][J]-Matrix[I,J]);
              end;
            end;
           

for I := 0 to 9 do
      begin
         if ((EuclDist[0][I]) < (EuclDist[1][I]))or ((EuclDist[0][I]) = (EuclDist[1][I]))   then
         begin
            target:=0 ;
            for j:= 0 to 9 do
            if target=(matrix[i][9]) then
              begin

f target=(matrix[i][9]) then