[C++] Metoda Laguerre'a - wyznaczanie zer wielomianu

[Maath]

Cześć,
Piszę program służący wyznaczaniu wszystkich rzeczywistych zer wielomianu, korzystam z algorytmu laguerre'a zamieszczonego tu:

Kod: Zaznacz cały

http://th-www.if.uj.edu.pl/zfs/gora/metnum11/wyklad08.pdf

,
a dokładnie ze wzoru (11) ze strony 14. Lecz przy niektórych wielomianach zdąża się że albo wartość pod pierwiastkiem jest ujemna albo albo dziele 0/0 co jest komplętną bzdurą...
Zamieszczam kod szukania zera:

[steal]

Nie wgłębiałem się w treść wykładu, ale czy wzór (11) nie jest dla podany dla liczb zespolonych? Wtedy podany przez Ciebie kod jest wybrakowany, bowiem nie obsługuje działań na tych liczbach.
Nie podałeś postaci wielomianu, który powoduje błąd.

[Maath]

Chociażby \(\displaystyle{ (x+1)^{3}}\), dla tego wielomianu powinno znajdować pierwiastek a tak nie jest..

[steal]

A obniżasz stopień wielomianu po znalezieniu miejsca zerowego?
Ciężko jest znaleźć błąd jeżeli widzi się tylko kawałek kodu źródłowego; z resztą wydaje mi się, że mało komu na forum będzie się chciało go analizować, poświęcając swój wolny czas. Niestety szukanie błędów to najmozolniejsza i najnudniejsza część programowania.

[Maath]

Obniżam stopień.. Zauważyłem że powodzenie algorytmu zalezy od punktu startowego..