Метод Ньютона-Рафсона в Mathematica

Я хочу реализовать метод Ньютона-Рафсона в Mathematica.

Вот мой код:

f[x] = x^3 - x^2 + 1

MetodaTangente[x0_, eps_] := Block[{p0, p1, dp, k},
   p0 = N[x0];
   p1 = p0;
   dp = 1;
   k = 0;
   While[dp > eps,
    p0 = p1;
    p1 = p0 - f[p0]/f'[p0];
    dp = Abs[p1 - p0];
    k = k + 1;
    ];
   Print[p1];
   ];

k подсчитывает, сколько было итераций.

Однако вот что происходит, когда я запускаю это:

Кажется, проблема с f. Что мне теперь делать?