typedef double Number;
class CubicRealPolynomial
{
public:
static Number computeDiscriminant(Number a, Number b, Number c, Number d);
static std::vector<Number> computeRealRoots(Number a, Number b, Number c, Number d);
private:
static std::vector<Number> computeRealRootsInner(Number a, Number b, Number c, Number d);
};
--------------
Number CubicRealPolynomial::computeDiscriminant(Number a, Number b, Number c, Number d)
{
Number a2 = a * a;
......
阅读全文