Napisao sam funkciju za testiranje funkcije roots. Radi na nacin da jos se prosledi polinom i koreni tog polinoma a funkcija od datih korena rekonstruise zadani polinom.
Znaci rekonstrukcija je

. Dalje se porede koeficijenti novonastalog polinoma i originalnog i racuna se greska kao

. I ovo je izlaz funkcije (greske).
Malo sam testirao funkciju i sve radi savrseno do nekog 70-og stepena polinoma. Za 94-i stepen polinoma sa koeficijentima 1 (

) ja dobijam gresku 90 i nesto! Sto vise povecavam stepen polinoma dobijam ogromnu gresku. Za polinom 500-og stepena greska je na 52-i stepen!
Da li moj algoritam za proveru ne valja ili pa njihova funkcija nije bas kako treba?
Funkciju koju sam napisao je u attach-u.
Eureka!
MILAN DJUKIC
D J U K A