用到了一元二次方程,一元三次方程的求解.
class QuarticRealPolynomial
{
public:
static Number computeDiscriminant(Number a, Number b, Number c, Number d, Number e);
static std::vector<Number> computeRealRoots(Number a, Number b, Number c, Number d, Number e);
private:
static std::vector<Number> original(Number a3, Number a2, Number a1, Number a0);
static std::vector<Number> neumark(Number a3, Number a2, Number a1, Number a0);
};
-----------------------------......
阅读全文