Скачиваний:
27
Добавлен:
09.12.2013
Размер:
1.91 Кб
Скачать
Program laba2;
Uses CRT;
Const a=0.5;
b=2;
eps=0.001;

Var x1,x2,pogr,pogrS,e,
int1,int2,int3,int,eps1:real;
i:integer;

function f(x:real):real;
begin
f:=x*x*x*(cos(3*ln(x)));
end;

Begin
clrscr;
Writeln('‚лзЁб«Ґ­ЁҐ Ё­вҐЈа «  ЎҐ§ Ў ­ЄЁа®ў ­Ёп');
WriteLn('h(i)':17,'Q(i)':21,'e(i)':17,'R(i)':17);
WriteLn('-----------------------------------------------------------------------------');
x1:=a;x2:=b;e:=eps;
int:=0;pogrS:=0;i:=1;
repeat
int1:=(x2-x1)/2*(f(x1)+f(x2));
int2:=(x2-x1)/4*(f(x1)+f(x2)+2*f((x1+x2)/2));
pogr:=abs(int2-int1)/3;
eps1:=e*(x2-x1)/(b-a);
if pogr<=eps1 then
begin
int:=int+int2;
pogrS:=pogrS+pogr;
WriteLn(' ', (x2-x1):5:5, ' | ' ,int:7:7, ' | ' , e*(x2-x1)/(b-a):5:9,' | ',pogr);
WriteLn;
x1:=x2;
x2:=b;
i:=i+1;
end
else x2:=(x2+x1)/2;
until b=x1;
writeln('Љ®«ЁзҐбвў® и Ј®ў=',i-1);
writeln('‘㬬  Ё­вҐЈа « = ',int:0:12);
writeln('‘㬬  Ї®ЈаҐи­®б⥩= ',pogrS:0:12,' < 0.001');
readln;
Writeln('‚лзЁб«Ґ­ЁҐ Ё­вҐЈа «  c Ў ­ЄЁа®ў ­ЁҐ¬');
WriteLn('h(i)':17,'Q(i)':21,'e(i)':17,'R(i)':17);
WriteLn('-----------------------------------------------------------------------------');
x1:=a;x2:=b;e:=eps;
int:=0;pogrS:=0;i:=1;
repeat
int1:=(x2-x1)/2*(f(x1)+f(x2));
int2:=(x2-x1)/4*(f(x1)+f(x2)+2*f((x1+x2)/2));
pogr:=abs(int1-int2)/3;
eps1:=e*(x2-x1)/(b-a);
if pogr<=eps1 then
begin
int:=int+int2;
pogrS:=pogrS+pogr;
WriteLn(' ', (x2-x1):5:5, ' | ' ,int:7:7, ' | ' , e*(x2-x1)/(b-a):5:9,' | ',pogr);
WriteLn;
e:=e+e*(x2-x1)/(b-a)-pogr;
x1:=x2;
x2:=b;
i:=i+1;
end
else x2:=(x2+x1)/2;
until b=x1;
writeln('Љ®«ЁзҐбвў® и Ј®ў=',i-1);
writeln('‘㬬  Ё­вҐЈа « = ',int:0:12);
writeln('‘㬬  Ї®ЈаҐи­®б⥩= ',pogrS:0:12,' < 0.001');
readln;
end.
Соседние файлы в папке ЛАБ_№2