Добавил:
Источник:
Uman
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз:
Предмет:
Файл:01 семестр / К экзамену-зачёту / Все задачи с экзамена в формате PAS / Alex / SERGE82
.PASprogram int;
var a,b,ik,is:real;
i,pk,ps:integer
procedure kvadr(n:integer,a,it:real;var ik:real;var p:integer);
begin
k:=1/2;p:=0;
repeat
k:=k*2;n:=n*k;
for i:=1 to n do
begin
x[i]:=a+(i-1/2)*h;ih:=ih+h*(exp(1/x)/(x*x));
end;
for i:=1 to 2*n do
x[i]:=a+(i-1/2)*h;
ihd:=ihd+(h/2)(exp(1/x)/(x*x));
ik:=ihd+1/3(ihd-ih);e:=abs(it-ik);p:=p+1;h:=h/2;
until e<0.0001;
end;
procedure sim(n:integer,a,it:real;var is:real;
var p:int);
var i:integer;x,f:array[1..70]of real,k,n,ih:real
h:=(b-a)/(2*n);
repeat k:=k*2;n:=n*k;
for i:=1 to 2*n+1 do
begin
x[i]:=a+h*i;f[i]:=(exp(1/x)/(x*x));end;
ih:=h/3*(f[1]+f[2n+1])
;i:=0;repeat i:=i+2;ih:=ih+4/3*h*f[i];until i=2*n;
i:=1; repeat i:=i+2;ih:=ih+2/3*f[i];until i:=2*n-1;
h:=h/2; for i:=1 to 4n+1 do
begin
x[i]:=a+h*i;f[i]:=(exp(1/x)/(x*x));
ihd:=h/3*(f[1]+f[4n+1];
i:=0;repeat i:=i+2;ihd:=ihd+(4/3)*h*f[i];until i=2n;
i:=1;repeat i:=i+2;ihd:=ihd+(2/3)*f[i]*h;until i=2n-1 ;
is:=ihd+(ihd-ih)/15
e:=abs(it-is);
until e<0.0001;
end;
begin{osn progr}
writeln('vvedite n,a,b');
readln(n,a,b);h:=(b-a)/n;it:=exp(1/a)-exp(1/b);
Соседние файлы в папке Alex
