Метод наименьших квадратов в паскале

program mnk;
uses crt;
var {Mx:Array[0..8] of integer;}Mx,My,Mz:array[0..7] of real;
i,j:integer;
A,B,k:real;
 
x1,x2,y1,y2,y3,r:real;
q1,q2,q3,q4,q5:real;
 
function ff(x:real;a:real;b:real):real;
var fun:real;
begin
fun:=a*(1/x)+b;
ff:=fun;
end;
 
begin
clrscr;
 
  Mx[0]:=3.80;   My[0]:=-21.41;
  Mx[1]:=0.25;   My[1]:=-9.9;
  Mx[2]:=0.48;   My[3]:=1.14;
  Mx[3]:=5.78;   My[4]:=-12.04;
  Mx[4]:=4.91;   My[5]:=-19.23;
  Mx[5]:=1.56;   My[6]:=-0.3;
  Mx[6]:=0.92;   My[2]:=11.26;
  Mx[7]:=5.73;   My[7]:=-19.56;
 
  for i:=0 to 7 do
  begin
  writeln('x[',i+1,']=',Mx[i]:2:2,' ','y[',i+1,']=',My[i]:2:2);
  end;
 
  for i:=0 to 7 do
  begin
  x1:=x1+Mx[i];
  x2:=x2+(Mx[i]*Mx[i]);
  y1:=y1+1/(My[i]*Mx[i]);
  y2:=y2+1/(My[i]);
  k:=k+((Mx[i])*(1/My[i]));
  end;
  A:=(8*k-k)/(8*x2-x2);
  B:=(y2-A*x1)/8;
  for i:=0 to 7 do
  begin Mz[i]:=ff(Mx[i],a,b);writeln(mz[i]);end;
  for i:=0 to 7 do
  y3:=y3+(Mz[i]-My[i])*(Mz[i]-My[i]);
  for i:=0 to 7 do
  begin
  q1:=q1+(Mx[i]-x1/8);
  q2:=q2+ff(Mx[i],a,b);
  q3:=q3+sqr((Mx[i]-q1/8));
  end;
  for i:=0 to 7 do
  begin
  q4:=q4+sqr(ff(Mx[i],a,b)-q2/8);
  q5:=q5+(ff(Mx[i],a,b)-q2/8);
  end;
  r:=q1*q5/sqrt(q3*q4);
 
  writeln('a=',a:2:5,' ','b=',b:2:5);
  writeln(Коэффицент корреляции равен=',' ', y3:2:2);
  writeln('Сумма квадратов отклонений =',' ', r:0:2);
 
  for i:=0 to 7 do
  writeln('x[',i+1,']=',Mx[i]:2:2,' ','z[',i+1,']=',mz[i]:2:2);
  readln;