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

Program Integral;
Uses CRT;
Const a=0;
b=1;
e=0.001;

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

Function func(x:real):real;
begin
func:=x*x/((1-x*x*x*x*x)*(1-x*x*x*x*x)*(1-x*x*x*x*x)+2);
end;

Begin
clrscr;
a_:=a;b_:=b;
int:=0;rest:=0;i:=1;
repeat
int1:=(b_-a_)/2*(func(a_)+func(b_));
int2:=(b_-a_)/4*(func(a_)+func(b_)+2*func((b_+a_)/2));
r:=abs(int1-int2)/3;
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('Int[',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('Znachenie integrala =',int:0:15);
writeln('Summa pogreshnostey =',rest:0:15);
readln;
end.