Метод хорд

[-=Andriushka=-]

подскажите пожалуйста...

Program MethodHorda;
Uses Crt;
Var
    a,b,t,x,eps,x1,pog_prac,pog_teor: Real;
    n,i:Integer;
    Q:array[1..100] of real;
Function func(x: Real): Real;
Begin
    func:=x*x*x - 18*x -83;
End;
Function f(x:Real):Real;//производная функции
Begin
     f:=x*x-18;
End;
Begin
    ClrScr;
    Write('Введите левую границу отрезка: ');
    ReadLn(a);
    Write('Введите правую границу отрезка: ');
    ReadLn(b);
    Write('Введите необходимое приближение (точность): ');
    ReadLn(eps);
    n:=0;
    Repeat
        inc(n);
        x:=a-func(a)*(b-a)/(func(b)-func(a));
        If func(a)*func(x)<=0 Then b:=x
                   Else a:=x;
    Until abs(func(x))<=eps;
    Writeln('Корень уравнения x = ', x:10:8);
    x1:=5.705;
    pog_prac:=abs(x1-x);
    writeln('Практическая погрешность = ',pog_prac:10:8);
    writeln('Число итераций = ',n);
    Repeat Until KeyPressed
End.

есть такой код...подскажите как вклинить можно подсчитывание теоретической погрешности а именно (|func(X[n+1])|/f(X[n]))<eps...а именно значение функции в следующей точке деленное на производную функции в предыдущей точке,и затем сравнивалось с eps,и если меньше eps то пошел счет в след точках...