Добавил:
Источник:
Uman
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз:
Предмет:
Файл:01 семестр / К экзамену-зачёту / Все задачи с экзамена в формате PAS / Alex / INTEG
.PASprogram int;
var a,b,ik,is,H,IT:real;
i,p,N:integer ;
procedure kvadr(n:INTEGER;a,H,it:real;var ik:real;var p:integer);
VAR K:INTEGER;X:ARRAY[1..3435] OF REAL;
IH,IHD,E:REAL;
begin
k:=1;p:=0;
repeat
;n:=n*k;
for i:=1 to n do
begin
x[i]:=a+(i-1/2)*h;ih:=ih+h*(exp(1/x[I])/(x[I]*x[I]));
end;WRITELN(P);
for i:=1 to 2*n do
x[i]:=a+(i-1/2)*h;
ihd:=ihd+(h/2)*(exp(1/x[I])/(x[I]*x[I]));
ik:=ihd+(1/3)*(ihd-ih);WRITELN(IK);READLN;e:=abs(it-ik);p:=p+1;h:=h/2;K:=K*2;
until e<0.0001;
end;
procedure sim(n:integer;a,B,it:real;var is:real;
var p:intEGER);
var i,K:integer;x,f:array[1..3435] of real;ih,E,IHD,H:real
; BEGIN
h:=(b-a)/(2*n);
repeat k:=1;n:=n*k;
for i:=1 to 2*n+1 do
begin
x[i]:=a+h*i;f[i]:=(exp(1/x[I])/(x[I]*x[I]));end;
ih:=h/3*(f[1]+f[2*n+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 4*n+1 do
begin
x[i]:=a+h*i;f[i]:=(exp(1/x[I])/(x[I]*x[I]));
ihd:=h/3*(f[1]+f[4*n+1]);END;
i:=0;repeat i:=i+2;ihd:=ihd+(4/3)*h*f[i];until i=2*n;
i:=1;repeat i:=i+2;ihd:=ihd+(2/3)*f[i]*h;until i=2*n-1 ;
is:=ihd+(ihd-ih)/15
; e:=abs(it-is);P:=P+1;K:=K*2;
until e<0.0001;
end;
begin{osn progr}
writeln('vvedite n,a,b');P:=0;
readln(n,a,b);h:=(b-a)/n;it:=exp(1/a)-exp(1/b);WRITELN('IT--',IT);
KVADR(N,A,H,IT,IK,P);
WRITELN('METOD KVADRATOV----- ',IK,'PROGONOK---- ', P);
P:=0 ;SIM(N,A,B,IT,IS,P) ;
WRITELN(' METOD SIMPSONA--- ',IS,'PROGONOK--- ',P);
READLN;READLN
END.
Соседние файлы в папке Alex
