выходные числа программы квадратного уравнения

извинения за стену текста.

Из книги, над которой я работаю, мне было дано следующее задание:

«Квадратные уравнения имеют вид: a x ^ 2 + b x + c = 0

Для их решения используется квадратичная формула: (-b ± √ (b ^ 2-4ac)) / 2a

Однако есть проблема: если b ^ 2-4ac меньше нуля, то он потерпит неудачу.
Напишите программу, которая может вычислить х для квадратного уравнения. Создайте функцию, которая выводит корни квадратного уравнения с учетом a, b, c.
Когда программа обнаружит уравнение без реальных корней, распечатайте сообщение. Откуда ты знаешь, что твои результаты правдоподобны? Можете ли вы проверить, что они верны?«

Я написал свою программу. Но все результаты, которые я получаю, странные числа. Например: 7.08141e + 268. Это использует 4, так что sqrt должно быть 2, я считаю.

РЕДАКТИРОВАТЬ: я изменил код и получаю правильные значения для sqrt a, b, c сейчас. Но, x1 и x2 выдают nan, что просто неправильно. Надеюсь, этот код будет немного чище и лучше оформлен.

double quad_x1(double b, double a, double c) {
double x1 = (-b + sqrt(b * b - 4 * a * c)) / (2 *a);
return x1;
}

double quad_x2(double b, double a, double c) {
double x2 = (-b - sqrt(b * b - 4 * a * c)) / (2 * a);
return x2;
}

double square_root_a(double a) {
double sq_a=sqrt(a);
if(sq_a<0){
throw invalid_argument("Input is less than zero.");
}
else {
return sq_a;
}
}

double square_root_b(double b) {
double sq_b=sqrt(b);
if(sq_b<0) {
throw invalid_argument("Input is less than zero.");
}
else {
return sq_b;
}
}

double square_root_c(double c) {
double sq_c=sqrt(c);
if(sq_c<0) {
throw invalid_argument("Input is less than zero.");
}
else {
return sq_c;
}
}

int main ()
{
try{
double a=4;
double b=4;
double c=4;

double sq_a = square_root_a(a);
double sq_b = square_root_b(b);
double sq_c = square_root_c(c);
double x1 = quad_x1( b, a, c);
double x2 = quad_x2( b, a, c);

cout<<"The square roots are: a = "<<sq_a<<" b = "<<sq_b<<" c = "<<sq_c<<'\n';
cout<<"The answer to the quadratic is:\n";
cout<<"x1 = "<<x1<<'\n';
cout<<"x2 = "<<x2<<'\n';
}

catch (invalid_argument) {
cerr<<"Invalid Input - Less than zero"<<'\n';
return 1;
}

return 0;
}