Методом половинного деления найти корни производной функции

var a,b,x,e,c: real;
function f(x:real):real;
  begin
  f:=7.5*(x*x*x*x)-1.2*(x*x*x)-5.1*sqr(x)+3*x-4.36;
  end;
begin
a:=-4;b:=0;e:=0.0001;
c:=(a+b)/2;
while abs(b-a)>e do 
  begin
  if f(a)*f(c)<0 then
  b:=c else a:=c;
  
  c:=(a+b)/2;
  end;
    
x:=(a+b)/2;
 
writeln('при x=',x:1:4,' с точностью 0,0001 ', 'f(x)=',f(x):1:4);
end.

const
  eps = 1e-5;
  LeftBorder  = -100;
  RightBorder = +100;
  StepSize    = 0.1;
  
function f(x : real) : real;
begin
  Result := 7.5 * x - 1.2;
  Result := Result * x - 5.1;
  Result := Result * x + 3;
  Result := Result * x - 4.36;
end;
 
procedure ansver(x : real);
begin
  WriteLn('f(', x:0:7, ') = ', f(x):0:7);
end;
 
procedure dihotomia(a, b : real);
var
  x : real;
  s : Integer;
begin
  s := sign(f(a));
  repeat
    x := (a + b) / 2;
    if sign(f(x)) = 0 then
      begin ansver(x); Exit; end
    else
      if sign(f(x)) = s then
        a := x
      else
        b := x;
  until abs(f(x)) < eps;
  ansver(x);
end;
 
var
  l, r : real;
  predf, curf : real;
  
begin
  // Будем искать от LeftBorder или ещё меньше, если f(LeftBorder) = 0
  l := LeftBorder;
  predf := f(l);
  while abs(predf) < eps do
    begin
      l -= StepSize;
      predf := f(l);
    end;
  
  // Разобъём функцию на знакопеременные интервалы
  r := LeftBorder + StepSize;
  repeat
    curf := f(r);
    while abs(curf) < eps do
      r += StepSize;
      
    if (sign(predf) <> sign(curf)) and (sign(predf) <> 0) and (sign(curf) <> 0) then
      begin
        // Решение на интервале [l, r]
        dihotomia(l, r);
        
        predf := curf;
        l := r;
      end;
    r += StepSize;
  until r > RightBorder;
  WriteLn('Поиск решений окончен.');
end.