Нелинейное уравнение методом Ньютона(реализовать метод)

[dmitrylezzz]

Товарищи,помогите реализовать метод в виде подпрограммы. Буду очень благодарен.

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.