выч методы не разобраны! / 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.