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

{Программа по теме 4, задание 2, вариант 1.}
{"Интегрирование в физхимии".}
{Интегрирование здесь проводится методом парабол(по формуле Симпсона).
Значение подынтегральной функции при s=0 было вычислено аналитически}
uses wincrt;
var
r,r0,Step,s,h,d,v,u,n,x,y,P,Tl,Th:real;
function{описание функции}
M(s:real):real;
begin
M:=exp(-0.0076*s*s)*sin(2.67*s-0.00004*s*s*s)*sin(s*r0)/s/s/r0;
end;
begin
writeln(' Начальное и конечное значения и шаг по r ?');
write( ' ');
readln(r0,r,Step);
writeln(' Нижний и верхний пределы и шаг интегрирования?');
write( ' ');
readln(Tl,Th,h);
writeln;{r0 и r-соответственно начальное и конечное значения r,
Step-шаг по r; Tl и Th-соответственно нижний и верхний пределы,
h- шаг интегрирования}
writeln(' r P(r)');
writeln;
repeat
n:=(Th-Tl)/h;
d:=1;v:=0;u:=0;
while d<=(n/2)do begin
x:=Tl+(2*d-1)*h;
y:=Tl+2*d*h;
v:=v+M(x);
u:=u+M(y);
d:=d+1;
end;
P:=h/3*(2.67-M(Th)+4*v+2*u);
writeln(r0:20:2,P:20:5);
r0:=r0+Step;
until r0>r;
end.