wolfram-mathematica - Mathematica 中的 Newton-Raphson 方法

Question

我想在 Mathematica 中实现 Newton-Raphson 方法。

这是我的代码：

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。我现在该怎么办？

Answer 1

像这样定义你的函数：

f[x_] := x^3 - x^2 + 1
MetodaTangente[-1, .000000000001]

> -0.754878

更多信息： http ://reference.wolfram.com/mathematica/tutorial/DefiningFunctions.html