Решение уравнения методом хорд и половинного деления

program zadanie;
 
function F(x: real): real;
begin
  F := exp(x) - 10 * x;
end;
 
function Dichotomy(a, b, e: real; var i: integer): real;
{ Метод дихотомии (деления отрезка пополам) }
var
  c: real;
begin
  i := 0;
  repeat
    c := (a + b) / 2;
    if (F(a) * F(c)) < 0 then
      b := c
    else
      a := c;
    i := i + 1;
  until (b - a) <= e;
  Result := (a + b) / 2;
end;
 
function Chord(a, b, e: real; var i: integer): real; { Метод хорд }
var
  x, x1: real;
begin
  i := 0;
  repeat
    x := (F(b) * a - F(a) * b) / (F(b) - F(a));
    if F(a) * F(x) > 0 then
      a := x
    else
      b := x;
    i := i + 1;
    x1 := (F(b) * a - F(a) * b) / (F(b) - F(a));
  until abs(x1 - x) < e;
  Chord := x;
end;
 
var
  a, b, eps, x: real;
  i: integer;
 
begin
  a := 0;
  b := 1;
  eps := 0.0001;
  x := Dichotomy(a, b, eps, i);
  writeln('X dichotomy=', x:0:5, ' chislo iter=', i);
  x := Chord(a, b, eps, i);
  writeln('X chord=', x:0:5, ' chislo iter=', i);
  readln;
 
end.

program Dihotomia;
 
uses crt;
 
var a, b, c, eps, a1, b1: real;
    k: integer;
    ch: char;
 
    Function f(x:real): real;
 
 begin
  f:=2*x*x*x-3*x*x-12*x; {Вводим выражение для вычисления функции}
 end;
 
 begin
  clrscr;
  {Ввод исходных данных}
  write ('Введите a:'); read(a);
  write ('Введите b:'); read(b);
  write ('Введите eps:'); read(eps); 
  writeln;
  writeln ('a=', a:3:7);
  writeln ('b=', b:3:7);
  writeln ('eps=', eps:3:9); 
  writeln;
 
  a1:=a; b1:=b; {Запоминать исходные данные}
  k:=0; 
 
  writeln('Решение уравнения "2x^3-3x^2-12x=0" методом дихотомии.');
  writeln;
  {Проверка введеных значений}
  if (a1=b1) then
  writeln('Ошибка! Левая часть интервала не должна равняться правой.');
  if (a1>b1) then
  writeln('Ошибка! Левая часть интервала не должна быть больше правой.');
 
 {Цикл}
 while(b-a)/2 >= eps {Пока условие больше или равно точности...}
  do {...выполять}
  begin
   k:=k+1; 
   c:=(a+b)/2; 
    if f(c)=0 then 
    begin    
     writeln ('Уравнение на интервале (',a1:4:3, '; ',b1:4:3,') имеет точный корень x=', c:4:5);
     break; {Прервать вывод, чтобы ответ не повторялся бесконечно}
    end;
    if f(a)*f(c)<0 then b:=c;
    if f(a)*f(c)>0 then a:=c;
  end;
 
 {Проверка и вывод значений}
 begin
  if f(c)<>0 then {Если функция не равно 0, то вывести примерный ответ, иначе корень точный (x=0)}
   begin
    writeln ('Уравнение на интервале (',a1:4:7, '; ',b1:4:7,') имеет корень x=', c:4:5);
    writeln;
    writeln ('Точность ',eps:4:9, ' достигнута за ' ,k, ' цикла(ов)', '.');
   end
 end;
 
readkey;
 
end.