Интегрирование Ньютон-Котес

[Gonzo]

Коэффициеты считаются заранее и после этого берутся из файлов (они в атаче). Подскажите кто что знает про адаптивный алгоритм, применительно к этой задаче.

 USES CRT;
var c,d,int,zn:real;
    N_:byte;

{funkcia nahozhdenia znachenia podintegralnoi funkcii}
Function Func(x:real):real;
begin
 Func:=exp((-1)*ln(1+sqr(exp(x))))*exp(x);
end;

{osnovnaia rabochaia funkcia}
FUNCTION NEWCOT (A,B:REAL;N:INTEGER):REAL;
TYPE MS = ARRAY [0..7] OF INTEGER;
VAR F1,F2:TEXT;
    X1,SUM,H,Y1:REAL;
    H1,H2:MS;
    I,J:BYTE;
BEGIN
 ASSIGN(F1,'KOEFF1.txt');
 RESET(F1);
 FOR J:=0 TO n DO READLN(F1,H1[J]);
 CLOSE(F1);

 ASSIGN(F2,'KOEFF2.txt');
 RESET(F2);
 FOR J:=0 TO n DO READLN(F2,H2[J]);
 CLOSE(F2);

 H:=(B-A)/N;
 X1:=A;
 SUM:=0.0;
 FOR I:=0 TO n DO
  BEGIN
   SUM:=SUM+FUNC(X1)*H1[I]/H2[I];
   X1:=X1+H;
  END;
 NEWCOT:=(B-A)*SUM;
END;

{osnovnaia programma}
begin

 c:=0;
 d:=1;
 N_:=7;

 int:=NEWCOT(c,d,N_);

 clrscr;
 writeln('INTEGRIROVANIE PO FORMULE NEWTONA-COTESA');
 writeln;
 writeln('C=',c:5:2);
 writeln('D=',d:5:2);
 writeln('N= ',N_);
 writeln;
 writeln('INTEGRAL = ',int:8:5);
 write('Press any key to continue...');
 readkey;
 clrscr;
 write('Vvedite prointegrirovannloe znachenie funkcii: ');
 readln(zn);
 zn:=abs(zn-int);
 writeln('Raznica: ',zn:8:5);
 write('Press any key to EXIT...');
 readkey;
end.