Документ взят из кэша поисковой машины. Адрес оригинального документа : http://www.students.chemport.ru/materials/pascal/p9_5.pas Дата изменения: Thu Jan 15 18:04:43 2009 Дата индексирования: Mon Oct 1 21:19:35 2012 Кодировка: Windows-1251

uses wincrt;
var
a,b,c,e,g,h,k,r,s1,s2,t,w:real;
x,y:array [0..15] of real;
f:text;
n,i:integer;
begin
assign(f,'p9_5.txt');
reset(f);
read(f,n);
for i:=1 to n do
begin
read(f,x[i],y[i]);
t:=t+x[i]/n;
w:=w+y[i]/n;
end;
s1:=0;s2:=0;h:=0;k:=0;c:=0;e:=0;g:=0;
for i:=1 to n do begin
s1:=s1+(x[i]-t)*(y[i]-w);
s2:=s2+(x[i]-t)*(x[i]-t);
h:=h+x[i];{сумма иксов}
k:=k+y[i];{сумма игреков}
c:=c+x[i]*y[i];{сумма x*y}
e:=e+x[i]*x[i];{сумма квадратов иксов}
g:=g+y[i]*y[i];{сумма квадратов игреков}
end;
a:=s1/s2;
b:=w-a*t;
r:=abs(((n*c)-(h*k))/sqrt((n*e-h*h)*(n*g-k*k)));{коэффициент корреляции}
writeln;
writeln(' a=',a:5:3,' b=',b:6:3,' r=',r:5:3);
writeln;
writeln(' N x y a*x+b a*x+b-y');
writeln;
for i:=1 to n do
begin
writeln(' ',i:3,x[i]:12:3,y[i]:12:3,a*x[i]+b:12:3,a*x[i]+b-y[i]:12:3);
end;
end.