найти площадь криволинейной трапеции

[icar1992]

люди выручайте, надо найти площадь криволинейной трапеции методом хорд, формула прямоугольников.
каждый раз выдает ошибку 207. посмотрите что не так, пожалуйста

{$F+}
program Integ;
uses crt;
const eps=0.000001;
      n=10;
type fun=function(x:real):real;
var x1,x2,x3,a,b,I1,I2,I3,I:real;w:integer;
function f1(x:real):real;
 begin
  f1:=3*(0.5/(x+1)+1);
 end;
function f2(x:real):real;
 begin
  f2:=2.5*x-9.5;
 end;
function f3(x:real):real;
 begin
  f3:=5/x;
 end;
procedure root (f:fun;  a, b, eps1: real; var x: real);
          var c:real;
begin
clrscr;
 repeat
   c:=(f(b)*a-f(a)*b)/(f(b)-f(a));
   if f(a)*f(c) > 0 then a := c
   else b := c;
 until abs((f(b)*a-f(a)*b)/(f(b)-f(a))-c)<eps1;
 x:=c;
end;
function integral(f:fun;a,b,eps:real):real;
  var h,s,s1,x1,x2,p:real;i,m:integer;
   begin
    m:=n;
    s1:=0;
    p:=1/3;
    repeat
     m:=m*2;
     h:=(b-a)/m;
     s1:=s;
     s:=0;
     i:=0;
     repeat
      i:=i+1;
      x1:=a+i*h;
      x2:=a+(i-1)*h;
      s:=s+f((x1+x2)/2)*h;
     until i=m;
    until abs(p*(s-s1))<eps;
    integral:=s;
   end;
begin
clrscr;
writeln('Enter a,b');
readln(a,b);
root(f1,f3,a,b,eps,x1);
writeln(w);
writeln('X1=',x1:3:6);
root(f2,f3,a,b,eps,x2);
writeln('X2=',x2:3:6);
root(f1,f2,a,b,eps,x3);
writeln('X3=',x3:3:6);
I1:=integral(f3,x1,x2,eps);
writeln('I1=',I1:3:6);
I2:=integral(f1,x1,x3,eps);
writeln('I2=',I2:3:6);
I3:=integral(f2,x2,x3,eps);
writeln('I3=',I3:3:6);
I:=I2-I1-I3;
writeln('I=',I:3:6);
readkey
end.

___________
Код нужно оформлять по правилам:
тегом [CODE]..[/СODE] (это кнопочка с решёточкой #)
Не забывайте об этом!
Модератор.