[C++] Pierwiastki równania 3 rzędu

Czy jest jakaś biblioteka, która mi policzy pierwiastki z równania 3 rzędu (ax^3+bx^2+cx+d=0), czy muszę sobie sam zaimplementować jakąś metodę numeryczną (np. splajny)?

zanotowane.pl

doc.pisz.pl

pdf.pisz.pl

zsf.htw.pl