Вычисление коэффициента корреляции между соответствующими данными двух одномерных массивов

const
  N=10;
 
type
  TArray=array[1..N] of real;
 
var
  X,Y:TArray;
  r:real;
  Nominator,Denominator:real;
  temp1,temp2:real;
  MeanX,MeanY:real;
  i:integer;
 
function CalcMean(const Arr:TArray):real;
var
  i:integer;
  res:real;
begin
  res:=0.0;
  for i:=1 to N do res:=res+Arr[i];
  res:=res/N;
  Result:=res;
end;
 
 
begin
  { генерируем выборку в диапазоне от 0 до 1 }
  Randomize;
  for i:=1 to N do
  begin
    X[i]:=Random(100);
    Y[i]:=Random(100);
  end;
  temp1:=X[1];
  temp2:=Y[1];
  for i:=2 to N do
  begin
    if X[i]>temp1 then temp1:=X[i];
    if Y[i]>temp2 then temp2:=Y[i];
  end;
  for i:=1 to N do
  begin
    X[i]:=X[i]/temp1;
    Y[i]:=Y[i]/temp2;
  end;
  { выводим выборки на экран }
  writeln('====== Array X ======');
  for i:=1 to N do write(X[i]:0:4,' ');
  writeln;
  writeln('====== Array Y ======');
  for i:=1 to N do write(Y[i]:0:4,' ');
  writeln;
  { считаем среднее для каждой выборки }
  MeanX:=CalcMean(X);
  MeanY:=CalcMean(Y);
  { вычисляем коэффициент корреляции }
  Nominator:=0.0;
  Denominator:=0.0;
  temp1:=0.0;
  temp2:=0.0;
  for i:=1 to N do
  begin
    Nominator:=Nominator+(X[i]-MeanX)*(Y[i]-MeanY);
    temp1:=temp1+sqr(X[i]-MeanX);
    temp2:=temp2+sqr(Y[i]-MeanY);
  end;
  Denominator:=sqrt(temp1*temp2);
  r:=Nominator/Denominator;
  { выводим на экран }
  writeln('r=',r:0:4);
  readln;
end.