выч методы не разобраны! / laba22 / INTX_B~1

program lab2;
uses crt;

function simp(a,b:real):real;

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

begin
simp:=(b-a)/6*(f(a)+4*f((b-a)/2)+f(b));
end;

var J1,J2,J,hi,ei,ai,bi,ri,h,R,eb:real;
i:integer;

const a=0.5;
b=2;
e=0.001;

begin
clrscr;
writeln(' ЃҐ§ Ў ­ЄЁа®ў ­Ёп Ї®ЈаҐи­®бвЁ:');
writeln('ЪДДДВДДДДДДДДДДДДДВДДДДДДДДДДДДДВДДДДДДДДДДДДДї');
writeln('і ь і  Ј hi і €­вҐЈа « Ji і Џ®ЈаҐи­. ei і');
writeln('ГДДДЕДДДДДДДДДДДДДЕДДДДДДДДДДДДДЕДДДДДДДДДДДДДґ');
i:=1;
h:=0;
ai:=a;
bi:=b;
R:=0;
repeat
hi:=bi-ai;
ei:=e/(b-a)*hi;
J1:=simp(ai,bi);
J2:=simp(ai,ai+hi/2)+simp(ai+hi/2,bi);
ri:=abs(J2-J1)/15;
if ri<=ei
then
begin
writeln('і',i:2,' ','і',hi:13:10,'і',J1:13:10,'і',ei:13:10,'і');
i:=i+1;
J:=J+J1;
ai:=bi;
bi:=b;
R:=R+ri;
h:=h+hi
end
else
bi:=ai+hi/2;
until ai=b;
writeln('ГДДДЕДДДДДДДДДДДДДЕДДДДДДДДДДДДДЕДДДДДДДДДДДДДґ');
writeln('ібг¬і',h:13:10,'і',J:13:10,'і',R:13:10,'і');
writeln('АДДДБДДДДДДДДДДДДДБДДДДДДДДДДДДДБДДДДДДДДДДДДДЩ');
readkey;
clrscr;
writeln(' C Ў ­ЄЁа®ў ­ЁҐ¬ Ї®ЈаҐи­®бвЁ:');
writeln('ЪДДДВДДДДДДДДДДДДДВДДДДДДДДДДДДДВДДДДДДДДДДДДДї');
writeln('і ь і  Ј hi і €­вҐЈа « Ji і Џ®ЈаҐи­. ei і');
writeln('ГДДДЕДДДДДДДДДДДДДЕДДДДДДДДДДДДДЕДДДДДДДДДДДДДґ');
J:=0;
i:=1;
h:=0;
ai:=a;
bi:=b;
eb:=e;
R:=0;
repeat
hi:=bi-ai;
ei:=eb/(b-a)*hi;
J1:=simp(ai,bi);
J2:=simp(ai,ai+hi/2)+simp(ai+hi/2,bi);
ri:=abs(J2-J1)/15;
if ri<=ei
then
begin
writeln('і',i:2,' ','і',hi:13:10,'і',J1:13:10,'і',ei:13:10,'і');
i:=i+1;
J:=J+J1;
ai:=bi;
bi:=b;
eb:=eb+ei-ri;
R:=R+ri;
h:=h+hi
end
else
bi:=ai+hi/2;
until ai=b;
writeln('ГДДДЕДДДДДДДДДДДДДЕДДДДДДДДДДДДДЕДДДДДДДДДДДДДґ');
writeln('ібг¬і',h:13:10,'і',J:13:10,'і',R:13:10,'і');
writeln('АДДДБДДДДДДДДДДДДДБДДДДДДДДДДДДДБДДДДДДДДДДДДДЩ');

readkey;
end.