Mathematica中的Newton-Raphson方法

我想在Mathematica中實現Newton-Raphson方法。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計數多次迭代如何在那裏。

然而，這裏是發生了什麼，當我運行此：

enter image description here

似乎有一個與f問題。我現在應該怎麼做？

來源

2013-06-04 Milos

這Mathworld頁 - http://mathworld.wolfram.com/NewtonsMethod.html - 包含一個更Mathematica-這種方法的慣用實現。所以如果內存服務的話，Mathematica文檔也是如此，但我現在沒有這些文檔。 –

定義你的函數是這樣，而不是：

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

> -0.754878

來源

2013-06-04 13:09:17

非常感謝！：） – Milos

Mathematica中的Newton-Raphson方法

回答

相關問題