Вычисление корней уравнений и определенных интегралов

Program zadacha1;
uses crt;
Const
    eps=0.001;
    inf=1.0e9;
    eps1=eps*0.01;
    eps2=eps*0.001;
type
    TF=function(x:real):real;
var
    res:real;
    r1,r2,r3:real;
 
 
{$F+}
{Function 1}
 
function f1(x:real):real;
begin
     f1:=ln(x);
end;
 
function ft1(x:real):real;
begin
     ft1:=sqr(x);
end;
 
function ft2(x:real):real;
begin
     ft2:=8/x;
end;
 
function f1p(x:real):real;
begin
 If x<>0 then
     f1p:=1/x
 else f1p:=Inf;
end;
 
function f1pp(x:real):real;
begin
     f1pp:=-1/(x*x);
end;
 
{Function 2}
 
function f2(x:real):real;
begin
     f2:=-2*x+14;
end;
 
function f2p(x:real):real;
begin
     f2p:=-2;
end;
 
function f2pp(x:real):real;
begin
     f2pp:=0;
end;
 
{Function 3}
 
function f3(x:real):real;
begin
 If x<>2   then
     f3:=(1/(2-x))+6
 else f3:=Inf;
end;
 
function f3p(x:real):real;
begin
 If x<>-2 then
     f3p:=1/(x*x-4*x+4)
 else f3p:=Inf;
end;
 
function f3pp(x:real):real;
begin
 If x<>2 then
     f3pp:=1/(x*x*x-6*x*x+12*x-8)
 else f3pp:=Inf;
end;
 
 
{Procedure ROOT}
 
procedure root(f,g,f1,g1:TF;a,b,d,eps:real;var x:real);
var
   c:real;
 
   function H(x:real):real;
   begin
        H:=f(x)-g(x);
           end;
 
    function H1(x:real):real;
   begin
        H1:=f(x)-g(x);
           end;
 
Begin
   repeat
   c:=d-H(d)/H1(d);
   write('c=',c,' ');
   write('a=',a,'b=',b,' ',(H((a+b)/2)-(H(a)+H(b))/2*H(a)),' ');
   writeln;
   If ((H((a+b)/2) > (H(a)+H(b))/2) and (H(a) > 0)) or ((H((a+b)/2) < (H(a)+H(b))/2) and (H(a) < 0)) then
      d:=b
   else
      d:=a;
      Write('|a-b|=',abs(b-a),' '); writeln;
   until abs(b-a) < eps;
   x:=a;
   writeln('ok!')
   end;
 
{function INTEGRAL}
 
function integral(f:TF;a,b,eps:real):real; {Simpson}
const
     p=0.06666;
var
     n,i:integer;
     h,sum,i1,i2:real;
 
Begin
     integral:=0;
     n:=16;
     h:=(b-a)/n;
     sum:=(f(a)+f(b))/2;
     for i:=1 to n-1 do
         sum:=sum+f(a+i*h);
     I2:=h*sum;
 
     repeat
     I1:=I2;
     n:=2*n;
     h:=(b-a)/n;
     sum:=(f(a)+f(b))/2;
     for i:=1 to n-1 do
         sum:=sum+f(a+i*h);
     I2:=h*sum;
     integral:=I2;
 
     until p*abs(I1-I2) < eps;
end;
 
     begin
     WriteLn('Test');
     WriteLn;
     WriteLn('f1=x^2');
     WriteLn('f2==8/x');
     WriteLn;
 
     root(ft1,ft2,eps1,x1);
 
     WriteLn('(f1,f2):x=',x1);
     WriteLn;
     WriteLn('integral y=x^2 ot 0 do 1 raven ',integral(ft1,0,1,eps2):5:3);
     WriteLn('__________');
     WriteLn;
 
     WriteLn('f1=ln(x)');
     WriteLn('f2=-2x+14');
     WriteLn('f3=1/(2-x)+6');
     WriteLn;
 
     I1:=integral(f1,x1,x3,eps2);
     I2:=integral(f2,x2,x3,eps2);
     I3:=intefral(f3,x1,x2,eps2);
 
     WriteLn('I1=',I1:6:4);
     WriteLn('I2=',I2:6:4);
     WriteLn('I3=',I3:6:4);
     WriteLn;
 
     WriteLn('S=',I1-I2-I3:5:3);
END.