Функции. Написать функцию решения квадратного уравнения

#include <iostream>
#include <cmath>
 
//индексы ax^2 + bx + c = 0, в  x1, x2 
//будут записываться корни, 
//bool для определения есть ли корни
bool func(double A, double B, double C, double* x1, double* x2 )
{
double descriminant = sqrt(B*B - 4*A*C, 2);//sqrt работает либо с double  либо с float
if(descriminant < 0) return false;
 
*x1 = (-B - sqrt(descriminant))/2*A;
*x2 = (-B + sqrt(descriminant))/2*A;
return true;
}
 
int main()
{
double A = 2.0, B = 3.0, C = 4.0;//после равно стоит ноль
double X1 = 0, X2 = 0;
if(!func(A, B, C, &X1, &X2))
std::cout<<"no roots!";
else 
std::cout<<"X1 = "<<X1<<"\nX2 = "<<X2;
 
return 0;
 }

#include <ctime>
#include <iomanip>
#include <iostream>
#include <cmath>
 
const int EPSILON = 0.00001;
 
class Complex {
 public:
  Complex() : real(0), imaginary(0) {};
  Complex(double real) : real(real), imaginary(0) {}
  Complex(const double &real, const double &imaginary) : real(real), 
    imaginary(imaginary) {}
  const double &getReal() const { return real; }
  const double &getImaginary() const { return imaginary; }
 private:
  double real, imaginary;
};
 
Complex operator+(const Complex &a, const Complex &b) {
  return Complex(a.getReal() + b.getReal(), a.getImaginary() + b.getImaginary());
}
 
Complex operator-(const Complex &a, const Complex &b) {
  return Complex(a.getReal() - b.getReal(), a.getImaginary() - b.getImaginary());
}
 
Complex operator*(const Complex &a, const Complex &b) {
  return Complex(
    a.getReal() * b.getReal() - a.getImaginary() * b.getImaginary(),
    a.getReal() * b.getImaginary() + a.getImaginary() * b.getReal());
}
 
Complex operator/(const Complex &a, const Complex &b) {
  return Complex(
    (a.getReal() * b.getReal() + a.getImaginary() * b.getImaginary()) /
    (b.getReal() * b.getReal() + b.getImaginary() * b.getImaginary()),
    (a.getImaginary() * b.getReal() - a.getReal() * b.getImaginary()) /
    (b.getReal() * b.getReal() + b.getImaginary() * b.getImaginary()));
}
 
std::ostream &operator<<(std::ostream &stream, const Complex &number) {
  stream << std::setprecision(4) << std::fixed << number.getReal();
  if (fabs(number.getImaginary()) > EPSILON) {
    stream << ((number.getImaginary() > 0)? "+" : "") << number.getImaginary() << "i"; 
  }
  return stream;
};
 
class QuadraticEquation {
 public:
  QuadraticEquation(const double &a, const double &b, const double &c)
    : x1(), x2(), discriminant(0), numberOfRoots(0), a(a), b(b), c(c) {
    solve();
  }
  double getA() const { return a; }
  double getB() const { return b; }
  double getC() const { return c; }
  double getDiscriminant() const { return discriminant; }
  Complex getX1() const { return x1; }
  Complex getX2() const { return x2; }
  double getNumberOfRoots() const { return numberOfRoots; }
 private:
  void solve() {
    if (a == 0 || fabs(a) < EPSILON) {
      if (b == 0 || fabs(b) < EPSILON) {
        numberOfRoots = 0;
      } else {
        numberOfRoots = 1;
        x1 = x2 = (-c) / b;
      }
      return;
    }
    discriminant = b * b - 4 * a * c;
    if (discriminant == 0 || fabs(discriminant) < EPSILON) {
      numberOfRoots = 1;
      x1 = x2 = (b * b) / (2 * a);
    } else if (discriminant < 0) {
      numberOfRoots = 2;
      x1 = Complex((-b) / (2 * a), sqrt(-discriminant) / (2 * a)); 
      x2 = Complex((-b) / (2 * a), -sqrt(-discriminant) / (2 * a)); 
    } else {
      numberOfRoots = 2;
      x1 = (-b + sqrt(discriminant)) / (2 * a);
      x2 = (-b - sqrt(discriminant)) / (2 * a);
    }
  }
  Complex x1, x2;
  double discriminant, numberOfRoots;
  double a, b, c;
};
 
std::ostream &operator<<(std::ostream &stream, const QuadraticEquation &equation) {
  std::cout << equation.getA() << "*x^2" << ((equation.getB() >= 0) ? "+" : "") <<
    equation.getB() << "*x" <<((equation.getC() >= 0) ? "+" : "") <<
    equation.getC() << " = 0";
}
 
int main(int argc, char **argv) {
  QuadraticEquation a(0, 0, 0);
  
  std::cout << a << std::endl <<
    "number of roots: " << a.getNumberOfRoots() << std::endl <<
    "x1 = " << a.getX1() << std::endl <<
    "x2 = " << a.getX2() << std::endl;
 
  std::cin.get();
  return 0;
}