Newton's method, hitta noll punkten

Ursprungligen inskrivet av Denniz PoP
Kan inte själv något om Perl, men kan detta vara något: http://www.ddj.com/web-development/184416187 ?

Ursprungligen inskrivet av vb
Jag antar att det är http://en.wikipedia.org/wiki/Newton's_method du är ute efter att implementera? Isf är det inte alls speciellt svårt om du har en lätt-deriverad funktion.

def do_newtons(x0, a):
    xn = x0 - f(x0)/f_der(x0)
    if abs(xn - x0) < a:
        return xn
    else:
        return do_newtons(xn, a)

(i python då jag inte kan perl, men du kan säkert läsa ut vad den gör ändå)

I mitt exempel har jag antagit att funktionen f som du vill hitta nollvärde för finns definierad och att du även tagit fram en derivata av den och kallat derivatan för f_der. x0 som skickas in till do_newtons() är en första gissning som du eller ditt program kan göra (beroende på hur din funktion ser ut), och a är ett gräns för när metoden har fastnat i ett minimum och inte gör fler framsteg. Ju lägre desto fler iterationer och längre tid kommer det att ta, men slutresultatet blir exaktare.