высчитывание площади

[nimbus2000]

Помогите пожалуйста разобраться!Нужно сделать программу ,которая считает интегралы.Какие способы на програмной языке бывают?
Есть пробный вариант программы ,но она не компилится.

program Project2;

{$APPTYPE CONSOLE}

uses
SysUtils;

Type
TFunc=function(x:real):real;
procedure root(a,b,eps:real;f:TFunc;var res:real);
begin
repeat
res:=(a+b)/2;
if (f(a)*f(b)<0) then b:=res else a:=res;
until (abs(b-a)<eps*abs(res)) or (abs(f(res))<eps);
end;
procedure integer(a,b,eps:real;f:TFunc;var res:real);
var
i,n:byte;
s,h,x:real;
begin
h:=1;
s:=(f(a)+f(b))*(b-a)/2;
repeat
s:=res;
n:=n*2;
h:=(b-a)/n;
res:=(f(a)+f(b))/2;
x:=a;
for i:=1 to (n-1) do
begin
x:=x+h;
res:=res+f(x);
end;
res:=res*h;
until abs((s-res)/res)<eps;
end;
var
x1,x2,eps:real;
f1,f2:TFunc;
Function f1(x:real):real;
begin
f1:=cos(x);
end;
Function f2(x:real):real;
begin
if (not(x=0)) then
f2:=1/x;
end;
begin
Function f3(x:real):real;
begin
f3:=f1(x)-f2(x);
end;
begin
readln(x1,x2,eps);
root(a,b,eps,f3,x1);
readln(a,b,eps);
root(a,b,eps,f3,x2);
integer(x1,x2,eps,f3,res);
writeln(res);
readln
end;
end.
Поодскажите ,где что неправильно а?