Не могу разобраться в коде. Решить нелинейное уравнение методом половинного деления.

[Леонардо]

Ребят вот код задачи на Delphi "Решить нелинейное уравнение методом половинного деления"

по поводу задачи, она написана для решение Тестового примера: f (x) = x^3 − sin(x) + cos(x)

unit Unit1;
 
interface
 
uses
  Windows, Messages, SysUtils, Variants, Classes, Graphics, Controls, Forms,
  Dialogs, StdCtrls, TeEngine, Series, ExtCtrls, TeeProcs, Chart;
 
type
  TForm1 = class(TForm)
    EquationSolveButton: TButton;
    Label1: TLabel;
    Label2: TLabel;
    Label3: TLabel;
    Edit1: TEdit;
    Edit2: TEdit;
    Edit3: TEdit;
    Label4: TLabel;
    Label5: TLabel;
    procedure EquationSolveButtonClick(Sender: TObject);
  private
    { Private declarations }
  public
    { Public declarations }
  end;
 
var
  Form1: TForm1;
 
implementation
 
{$R *.dfm}
//описание заданной функции
function equation(x:extended):extended;
begin
result:=x*x*x-sin(x)+cos(x);
end;
//------------------------------------------------------------------
{****Решение произвольного нелинейного уравнения методом
половинного деления****}
 
procedure TForm1.EquationSolveButtonClick(Sender: TObject);
var
  x,error:array of double;
  x0,x1,eps:real;
  i,p:integer;
begin
 
  label4.Caption:='';
  label5.Caption:='';
 
 
  x0:=StrToFloat(edit1.Text);
  x1:=StrToFloat(edit2.Text);
  eps:=StrToFloat(edit3.Text);
 
 setlength(x,10000);
 setlength(error,10000);
 
  x[0]:=x0;
  i:=0;
 
  begin
 
    repeat
 
      i:=i+1;
      error[i]:=abs(equation(x[i-1]));
    until (error[i]<eps) or (i>9990);
 
    p:=i;
 
    label4.Caption:='x='+FloatToStrF(x[p],ffFixed,5,5);
    label5.Caption:='Число итераций равно'+IntToStr((p));
//================================
        x[1]:=x1;
 
      x[2]:=(x[0]+x[1])/2;
 
        i:=2;
 
        repeat
 
          if equation(x[i])*equation(x[i-1])<=0 then
            x[i+1]:=(x[i]+x[i-1])/2
          else
            if equation(x[i])*equation(x[i-2])<=0 then
              x[i+1]:=(x[i]+x[i-2])/2
            else
              if equation(x[i])*equation(x[i-3])<=0 then
                x[i+1]:=(x[i]+x[i-3])/2
              else
                x[i+1]:=(x[i]+x[i-4])/2;
 
          error[i]:=abs(x[i+1]-x[i]);
          i:=i+1;
 
        until (error[i-1]<2*eps) or (i>9990);
 
        p:=i;
 
        label4.Caption:='x='+FloatToStrF(x[p],ffFixed,5,5);
 
        label5.Caption:='Число итераций равно  '+IntToStr(p);
 
end
end;
end.

Объясните мне пожалуйста весь код, что-как считает, кто за что отвечает ...