Решение интеграла методом трапеций

Program zadanie2;
   Uses CRT;
   function f(x:real):real;
      begin
           f:=exp(x)*1/2*(1+cos(2*x));
      end;
   function Ff(x:real):real;
      begin
         Ff:=(exp(x)/2)/(1+((2*sin(2*x)+cos(2*x))/5))-0.6;
      end;
   function N(a:real;b:real):real;
      var
         x:real;
      begin
           N:=Ff(b)-Ff(a);
           x:=a;
           a:=0;
           b:=pi;
      end;
   function trapec(a:real;b:real;num:integer):real;
      var
         x,s,h:real;
         i:integer;
      begin
           s:=(f(a)+f(b))/2;
           x:=a;
           h:=(b-a)/num;
           for i:=1 to num-1 do
           begin
              x:=x+(b-a)/num;
              s:=s+f(x);
           end;
           trapec:=s*((b-a)/num);
      end;
   const
        M=0.036;
        E=0.001;
 
   var
      a,b,h,i,z:real;
      num:integer;
      BEGIN
           clrscr;
           Writeln('Vychislenie opredelennogo integrala metodom trapecii');
           Writeln;
           Writeln('x[i]         znachenie x[i]        znachenie pervoobraznoi');
           i:=0;
           z:=1;
           While i<=pi do
              begin
              Writeln(' ',z:2:0,'               ',i:3:4,'                ',ABS(Ff(i)):3:4);
              Writeln;
              i:=i+pi/6;
              z:=z+1;
              end;
           a:=0;
           b:=pi;
           num:=60;
           h:=(b-a)/num;
           Writeln('    Predeli Kolichestvo predelov   shag   znachenie Otkloneniya');
           Writeln('Integrirovaniya  integrala ot tochnogo znacheniya');
           Writeln;
           i:=0;
           num:=60;
           Writeln('a=0','   ','b=pi','   ','num=',num:2,'   ','h=',h:3:2,'    ',trapec(a,b,num):7:5,'        ',
           N(a,b)-trapec(a,b,num));
           num:=154;
           h:=(b-a)/num;
           Writeln;
           Writeln('a=0','   ','b=pi','    ','num=',num:2,'    ','h=',h
           :3:2,'   ',trapec(a,b,num):7:5,'        ',
           N(a,b)-trapec(a,b,num));
           Writeln;
           Writeln('Rezultat vychisleniya po phormule Njutona-Lejbnica:N= ',N(a,b):7:5);
           readln;
           END.
      END.