выч методы не разобраны! / LAB2-И~1 / МЕТОДС~1 / INT

Program IntAdapt;
Uses CRT, Maths;
Const a=0;
b=3;
e=0.001;
t1=0.5773;
t2=0.7071;

Var a_,b_,r,rest,
int1,int2,int:real;
i:integer;

Function x(t,a,b:real):real;
begin
x:=(a+b)/2+(b-a)/2*t;
end;

Function f(x:real):real;
begin
f:=sqr(x)/(pwr(1-pwr(x,5),3)+2);
end;

Begin
clrscr;
a_:=a;b_:=b;
int:=0;rest:=0;i:=1;
repeat
int1:=(b_-a_)/2*(f(x(t1,a_,b_))+f(x(-t1,a_,b_)));
int2:=(b_-a_)/3*(f(x(t2,a_,b_))+f(x(-t2,a_,b_))+f(x(0,a_,b_)));
r:=abs(int1-int2)/15;
if r<=e*(b_-a_)/(b-a) then
begin
int:=int+int2;
rest:=rest+r;
writeln('a[',i,']=',a_:0:6);
writeln('b[',i,']=',b_:0:6);
writeln('I[',i,']=',int2:0:15);
writeln('R[',i,']=',r:0:15);
writeln('E[',i,']=',e*(b_-a_)/(b-a):0:15);
writeln;
readkey;
a_:=b_;
b_:=b;
i:=i+1;
end
else b_:=(b_+a_)/2;
until b=a_;
writeln('SumI=',int:0:15);
writeln('SumR=',rest:0:15);
readln;
end.