Нелинейное уравнение методом Ньютона (исправить) - Turbo Pascal

Формулировка задачи:

В программе необходимо метод оформить в виде подпрограммы. Вывод на экран должен быть полным, с указание уравнения и количества итераций. Помогите сделать,пожалуйста.

Program Newton;
Function f(x:real):real;
begin 
  f:=x*x*x+x*x+x+1; 
end;
Function f1(x:real):real;
begin 
  f1:=3*x*x+2*x+1; 
end;
var
  eps, x: real; 
BEGIN
write ('Начальное приближение:'); readln(x);
write ('Точность:'); readln(eps);
Repeat
  x:=x- f(x)/f1(x);
Until abs(f(x))<=eps;
writeln ('x=',x:3:3,' f(x)=',f(x):4:4);
readln(x);
END.

Решение задачи: «Нелинейное уравнение методом Ньютона (исправить)»

textual
                            Листинг программы 
                        

                            repeat
write ('Начальное приближение в интервале [-2;0] x:');
readln(x);
until(x>=-2)and(x<=0);
write ('Точность:'); readln(eps);
y:=Newton(x,eps,i);
writeln ('x=',y:6:3,' f(x)=',f(y):7:4,' количество итераций=',i);
readln
                        