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

Uses crt;
var main:real;

Function Func(x:real):real;
begin
 Func:=exp(-1*ln(1-sqr(exp(x))))*exp(x);
end;

FUNCTION NEWCOT(C,D:REAL;N:INTEGER):REAL;
TYPE MS=ARRAY[0..7] OF INTEGER;
     MS2=ARRAY [3..7] OF MS;
     FF=FILE OF MS2;

VAR H1,H2:MS2;
    FC1,FC2:TEXT;
    X1,SUM,H,Y1:REAL;
    I,J:Byte;
    name1,name2:string;

 BEGIN
  writeln('Vvedite imia, soderzhaschego chislitili koefficientov formul Niutona-Kotesa: ');
  readln(name1);
  clrscr;
  writeln('Vvedite imia, soderzhaschego znamenateli koefficientov formul Niutona-Kotesa: ');
  readln(name2);
  clrscr;
  ASSIGN (FC1,name1);
  ASSIGN (FC2,name2);
  RESET(FC1);
  RESET(FC2);
  FOR I := 3 TO 7 DO
  FOR J := 0 TO 7 DO READ(FC1,H1[I,J]);
  FOR I := 3 TO 7 DO
  FOR J := 0 TO 7 DO READ(FC2,H2[I,J]);
  CLOSE(FC1);
  CLOSE(FC2);
  H:=(D-C)/N;
  X1:=C;
  SUM:=0.0;
  FOR I:=0 TO 7 DO
   BEGIN
    SUM:=SUM+FUNC(X1)*H1[N,I]/H2[N,I];
    X1:=X1+H;
   END;
  NEWCOT:=(D-C)*SUM;
END;

begin
clrscr;
main:=NEWCOT(0,1,7);
writeln('Znachenie integrala: ',main:6:2);
readln;
end.