Универсальная программа для решения уравнений методом половинного деления

uses crt;
function F(x:real):real;
begin
F:=cos(2/x)-2*sin(1/x)+1/x;
end;
var a,b,c,x,e:real;
begin
clrscr;
a:=1;
b:=2;
e:=0.001;
repeat
c:=(a+b)/2;
if F(a)*F(c)<=0 then b:=c
else a:=c;
until b-a<e;
x:=(a+b)/2;
writeln('x=',x:0:3);
readln
end.

unit unit1;
interface
type Func=function(x:real):real;
procedure Bisec(F:Func;a,b,e:real);
implementation
procedure Bisec(F:Func;a,b,e:real);
var c,x:real;
begin
if F(a)*F(b)<=0 then
 begin
  repeat
   c:=(a+b)/2;
   if F(a)*F(c)<=0 then b:=c
   else a:=c;
  until abs(b-a)<e;
  x:=(a+b)/2;
  writeln('x=',x:0:3);
 end
else writeln('Интервал выбран неверно');
end;
end.

uses crt,unit1;
{$F+}
function Y(x:real):real;
begin
Y:=x+sqrt(x)+exp(ln(x)/3)-2.5;
end;
function Z(x:real):real;
begin
Z:=cos(2/x)-2*sin(1/x)+1/x;
end;
{$F-}
begin
clrscr;
writeln('Функция y=x+x^(1/2)+x^(1/3)-2.5');
writeln('Интервал [0.4;1] точность 0.001');
Bisec(Y,0.4,1,0.001);
writeln('Функция y=cos(2/x)-2*sin(1/x)+1/x');
writeln('Интервал [1;2] точность 0.001');
Bisec(Z,1,2,0.001);
readln
end.