int quarticRoots(const double A, const double B, const double C, const double D,
        const double E, double s[4]);