выч методы не разобраны! / LAB2-И~1 / 2ЛАБА~1 / МОЕ / 1

program thebuch;
uses crt;
const e=0.001;
a=0;
b=1;
t=0.7947;
var xa,xb,e1,s1,s2,eq:real;
q1:array[0..100] of real;
q2:array[0..100] of real;
ep1:array[0..100] of real;
ep2:array[0..100] of real;
h1:array[0..100] of real;
h2:array[0..100] of real;
i,i1,i2:integer;

function f(x:real):real;
begin
f:=(x*x)*sin(4*ln(x))/ln(x)
end;

function epsi(a,b:real):real;
var f4,x9:real;
begin
x9:=(b+a)/2;
f4:=((2*x9*sin(4*ln(x9))/ln(x9))+(4*x9*cos(4*ln(x9))/ln(x9))+(x9*sin(4*ln(x9))/ln(x9)*ln(x9)));
epsi:=(1/2721600)*exp(5*ln(b-a))*f4;
end;

function Theb(a1,b1:extended):extended;
var x1,x2,x3:real;
begin
x1:=(b1+a1)/2-(b1-a1)/2*t;
x2:=(b1+a1)/2+(b1-a1)/2*t;
x3:=(b1+a1)/2;
Theb:=(0.5*(f(x1)+f(x2)+f(x3)));
end;

procedure Out;
begin
ClrScr;
Writeln('‚лзЁб«Ґ­ЁҐ Ё­вҐЈа «  ЎҐ§ Ў ­ЄЁа®ў ­Ёп');
WriteLn('j':5, '|', 'h(i)':17,'Q(i)':21,'e(i)':5);
WriteLn('-----------------------------------------------------------------------------');
WriteLn('-----------------------------------------------------------------------------');
for i:=0 to i1-1 do begin
WriteLn(' ':5, ' ', h1[i]:5:5, ' | ', q1[i]:7:7, ' | ' , ep1[i]:5:11);
end;
WriteLn;
WriteLn('‡ ¤ ­  Ї®ЈаҐи­®бвм : ',e:5:5);
WriteLn('‡­ зҐ­ЁҐ Ё­вҐЈа «  : ',S1:7:7);readln;
Writeln('‚лзЁб«Ґ­ЁҐ Ё­вҐЈа «  c Ў ­ЄЁа®ў ­ЁҐ¬');
WriteLn('j':5, '|', 'h(i)':17,'Q(i)':21,'e(i)':5);
WriteLn('-------------------------------------------------------------------------------');
WriteLn('-------------------------------------------------------------------------------');
for i:=0 to i2-1 do begin
WriteLn(' ':5, ' ', h2[i]:5:5, ' | ', q2[i]:5:5, ' | ' , ep2[i]:5:11);
end;
WriteLn('‡ ¤ ­  Ї®ЈаҐи­®бвм : ',e:5:5);
WriteLn('‡­ зҐ­ЁҐ Ё­вҐЈа «  : ',S2:5:5);readln


end;


BEGIN
xa:=a;xb:=b;i:=0;
repeat
if abs(e*(b-a)/(xb-xa))>abs(epsi(xa,xb)) then
begin q1[i]:=Theb(xa,xb);ep1[i]:=abs(epsi(xa,xb));h1[i]:=xb-xa;xa:=xb;xb:=b;i:=i+1;i1:=i;end
else xb:=(xa+xb)/2;
if xa<>xb then eq:=epsi(xa,xb);
until xa=b;
{Bank}
xa:=a;xb:=b;i:=0;e1:=e;
repeat
if abs(e1*(b-a)/(xb-xa))>abs(epsi(xa,xb))
then begin q2[i]:=Theb(xa,xb);ep2[i]:=abs(epsi(xa,xb));h2[i]:=xb-xa;e1:=e1-epsi(xa,xb);xa:=xb;xb:=b;i:=i+1;i2:=i;end
else xb:=(xa+xb)/2;
if xa<>xb then eq:=epsi(xa,xb);
until xa=b;
s1:=0;
for i:=0 to i1 do s1:=q1[i]+s1;
{writeln(s1);}s2:=0;
for i:=0 to i2 do s2:=q2[i]+s1;
{writeln(s2); }
out;
END.