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

uses wincrt;
var
a,b,c,e,Ea,g,h,k,k0,r,s1,s2,t,w:real;
x,y:array [0..9] of real;
f,q:text;
n,i:integer;
begin
assign(f,'p10_5.txt');
reset(f);
read(f,n);
for i:=1 to n do
begin
read(f,x[i],y[i]);
x[i]:=1/x[i];y[i]:=ln(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:8:1,' b=',b:5:1,' r=',r:7:4);
writeln;
writeln(' N 1/T ln k ln k выч Погрешность');
writeln;
assign(q,'p10_5_1.txt');rewrite(q);
for i:=1 to n do
begin
writeln(i:3,x[i]:12:6,y[i]:12:4,a*x[i]+b:12:4,y[i]-a*x[i]-b:12:4);
writeln(q,x[i]:12:6,y[i]:12:4,a*x[i]+b:12:4);
end;
close(q);
k0:=exp(b);Ea:=-a*8.31441;
writeln;
writeln(' k0 =',k0:9,' Ea =',Ea:10:2);
end.