Нахождение корней квадратного уравнения

#include <cmath>
#include <iostream>
 
int main()
{
    int a, b, c;
 
    std::cout << "A program for solving quadratic equations of the form ax^2 + bx + c = 0 (relative to x).\n"
              << "Enter coefficients \'a\' \'b\' \'c\': ";
 
    std::cin >> a >> b >> c;                    // считываем коэффициенты введенные пользователем
 
    while (a == 0) {                            // коэффициент при первом слагаемом уравнения не может быть равным нулю
        std::cout << "Coefficient \'a\' cannot be zero!\nEnter coefficient \'a\': ";
        std::cin >> a;
    }
 
    double d = pow(b, 2) - 4 * a * c;           // вычисляем дискриминант по формуле: d = b^2 − 4ac
 
    if (d > 0) {                                // дискриминант больше нуля: уравнение имеет два корня
        double x1 = (-b - sqrt(d)) / (2 * a);   // вычисляем первый корень по формуле: x1 = (-b - √d) / (2a)
        double x2 = (-b + sqrt(d)) / (2 * a);   // вычисляем второй корень по формуле: x2 = (-b + √d) / (2a)
 
        std::cout << x1 << " " << x2 << std::endl;
    }
    else if (d == 0) {                          // дискриминант равен нулю: уравнение имеет один корень
        double x = (-b + sqrt(d)) / (2 * a);    // вычисляем корень по формуле: x1 = (-b ± √d) / (2a)
 
        std::cout << x << std::endl;            // повторный вывод корня, сделан из за специфики проверки задания на сайте
    }
    else {                                      // дискриминант меньше нуля: уравнение не имеет корней
        std::cout << "No real roots" << std::endl;
    }
 
    return 0;
}