Решить нелинейное уравнение методом простой итерации

uses crt;
function F(x:real):real;
begin
F:=x*x*x+cos(8*x);
end;
procedure Iter(a,b,e:real);
var i,n:longint;
    x1,x2,f1,f2:real;
begin
writeln('Точность=',e:0:8);
if F(a)=0 then
 begin
  writeln('x=a=',a:0:4);
  readln;
  exit;
 end
else if F(b)=0 then
 begin
  writeln('x=b=',b:0:4);
  readln;
  exit;
 end
else if F(A)*F(B)>0 then WriteLn('Границы интервала заданы некорректно!')
else
 begin
  n:=round((b-a)/e);
  x1:=a; f1:=F(x1);
  for i:=1 to n do
   begin
    x2:=x1+e;
    f2:=F(x2);
    if f1*f2 <= 0 then writeln('x = ', (x2 + x1)/2:0:8);
    x1:=x2;
    f1:=f2;
   end;
 end;
end;
var a,b:real;
begin
clrscr;
a:=0.5;
b:=0.6;
Iter(a,b,0.00001);
Iter(a,b,0.0000001);
Iter(a,b,0.00000001);
readln
end.

{Решение уравнения x^3 + Cos(8*x) = 0 методом простых итераций}
Uses Crt;
Var x1, x2, f1, f2, k, eps, h :real;
    xX : array[1..20] of real;
    i, j, n : integer;
 
Function Func(arg : real): real;
begin
  Func:=arg*arg*arg + Cos(8*arg);
end;
 
Procedure PrimeIter;     {Метод простых итераций}
begin
  i:=0;
  repeat
    inc(i);
    x2:=x1 - k*Func(x1);
    If abs(x2-x1)<eps then Break;
    x1:=x2;
    If i>100000 then Break;      {Защита от зацикливания}
  until false;
end;
 
Begin
  h:=0.01;      {Шаг отделения корней}
  i:=0;
  n:=0;
  {Отделение корней}
  x1:=-20;
  f1:=Func(x1);
  repeat
    f2:=f1;
    x1:=x1+h;
    f1:=Func(x1);
    If i>100000 then Break;      {Защита от зацикливания}
    inc(i);
    If f1*f2<0 then
      begin
        Inc(n); xX[n]:=x1-h;     {Здесь корень}
      end;
  until false;
 
  {Уточнение корней с заданной точностью}
  k:=0.1;            {Коэффициент простых итераций}
  eps:=0.000001;     {Заданная точность}
  For j:=1 to n do
    begin
      x1:=xX[j];
      PrimeIter;
      f1:=Func(x1);
      WriteLn('Число итераций: ',i:6,' | X = ',x1:15:11,' | F(x) = ',f1:15:11);
    end;
  ReadLn;
end.