Функция находящая корни квадратного уравнения

#include <iostream>
#include <math.h>
#include <iomanip>
using namespace std;
 
 
int main()
{
    double a=1.5, b=1.5, c=1.5, d, x1, x2;
    cin >> a >> b >> c;
    d = (b * b) - (4 * a * c);
    if(d == 0)
    {
        x1=(-b+sqrt(d))/(2*a);
        x2=(-b-sqrt(d))/(2*a);
        cout << fixed << setprecision(6) << x1 << endl;
        cout << fixed << setprecision(6) << x2 << endl;
    }
    if(d > 0)
    {
        x1=(-b+sqrt(d))/(2*a);
        x2=(-b-sqrt(d))/(2*a);
        cout << fixed << setprecision(6) << x1 << endl;
        cout << fixed << setprecision(6) << x2 << endl;
    }
    if(d < 0)
    {
        cout<<" "<<endl;
    }    return 0;
}

#include <iostream>
#include <cmath>
#include <iomanip>
 
void getRoots(double a, double b, double c, int &n, double &x1, double &x2) {
    double d = b*b - 4*a*c;
    n = (d > 0 ? 2 : !(d<0) );
    x1 = (sqrt(d) - b) / ( 2*a );
    x2 = ( d>0 ? (-sqrt(d)-b) / (2*a) : x1);
}
 
using namespace std;
 
int main()
{
    int n = 0;
    double a, b, c, x1, x2;
    cin >> a >> b >> c;
    getRoots(a, b, c, n, x1, x2);
    cout << "n = " << fixed << setprecision(6) << n << ", x1 = " << x1 << ", x2 = " << x2 << endl;
}

#include <iostream>
#include <cmath>
#include <iomanip>
using namespace std;
void numswap(double *a, double *b) {
  double c=*a;
  *a=*b;
  *b =c;
}
int quadratic_solve(double a, double b, double c, double *x1, double *x2) {
  double d = pow(b,2) - 4*a*c;
  if (d<0) return 0;
  if (fabs(d) <= 1e-7) {
    *x1=-b/(2*a);
    return 1;
  }
  *x1 = (-b + sqrt(d)) / (2*a);
  *x2 = (-b - sqrt(d)) / (2*a);
  if (*x1 > *x2) numswap(x1,x2);
  return 2;
}
 
int main() {
  double a, b, c, x1, x2;
  cin >> a >> b >> c;
  int x =quadratic_solve(a, b, c, &x1, &x2);
  cout <<fixed << setprecision(6);
  if (x==1) cout<<x1;
  if (x==2) cout<<x1<<" "<<x2;
}