Fixed Kapenga's interface code.  Make sure to change qs = -2 to qs = qsolve_roots(&coef, &root) before testing it with qsolve_roots().