Решить нелинейное уравнение методом Ньютона (Паскаль)

[dmitrylezzz]

Заданный метод должен быть реализован в виде отдельной подпрограммы. Решить нелинейное уравнение методом Ньютона.

program Newton;
uses crt; 
function f(x:real):real; 
begin
 f:=sqr(sqr(x))-5*sqr(x)-x+1;
end;
function f1(x:real):real; 
begin                           
 f1:=4*x*sqr(x)-10*x-1;
end;
var a,b,x,e,en:real;
    i:integer;
begin
 clrscr; 
 writeln ('Решение нелинейного уравнения методом Ньютона');
 writeln ('Уравнение x^4+5x^2-x+1=0');
 write ('Введите левую и правую границы интервала:');
 read (a,b);
 write ('Введите требуемую точность решения:');
 read (e);
 writeln ('Решение:');
 writeln ('Номер шага Значение X');
 en:=abs(a-b);
 x:=b;
 i:=1;
 while (abs(en)>e) do begin 
  x:=x-f(x)/f1(x);
  writeln (i:10,x:20:14);
  en:=abs(x-b);
  b:=x; 
  i:=i+1; 
 end;
end.

Программа работает не верно. Кроме того уравнение не имеет действительных корней. Буду признателен за помощь.



________
Код нужно оформлять по правилам:
тегом [CODE]..[/СODE]
(это кнопочка на панели форматирования с решёточкой #)
Не забывайте об этом!

Модератор.