C++二次方程输出：-1.#IND

Question

我的代码工作正常： (1, -2, -8)，当我输入 a=1 b=0 c=1 时，它给我上面的错误，

这是我的代码：

    double x=0,a=0,b=0,c=0,d=0;
    complexType solu1;
    complexType solu2;


    cout << "\n\nEnter values of quadratic a,b,c:";
    cin >> a >> b >> c;

    double solution1 = (-1.0 * b) + (sqrt((b * b) - (4 * a * c)));
    solu1 = solution1 / (2*a);

    cout << setprecision(5) << solu1;

    double solution2 = (-b) - (sqrt((b*b) - (4 * a * c)));
    solu2 = solution2 / (2*a);
    cout << setw(5) << setprecision(5)  << solu2;

我该如何解决这个问题？

Answer 1

您正在尝试找到 x^2 + 1 = 0 的实数解，该解只有虚数解。

如果您想求解所有二次方程，那么您需要检查判别式是否为正，以确定解是实数还是复数：

double d = b*b - 4*a*c;
if (d >= 0) {
    double sol1 = (-b + sqrt(d))/(2*a);
    double sol2 = (-b - sqrt(d))/(2*a);
    std::cout << sol1 << ", " << sol2 << '\n';
} else {
    double real = -b/(2*a);
    double imag = sqrt(-d)/(2*a);
    std::cout << real << " +/- " << imag << "i\n";
}

如果您愿意，您可以使用 std::complex 做一些更简洁的事情。

Answer 2

输入的 sqrt((b*b) - (4 * a * c)) 为 sqrt(-4)。根据 http://www.cplusplus.com/reference/clibrary/cmath/sqrt/ , 如果参数为负数，则会发生域错误，将全局变量 errno 设置为值 EDOM。 在这种情况下，我没有看到它返回的内容的定义。不管怎样，这是错误的。

我发现您的代码中有 complexType 。如果这是 std::complex 的 typedef，那么代码很容易修复。

complexType top = b*b - 4*a*c;
solu1 = (-b + sqrt(top)) / (2*a);
solu2 = (-b - sqrt(top)) / (2*a);

由于 std::sqrt 具有 std::complex 的重载。