Программа решения квадратного уравнения

#include <iostream>
 
using namespace std;
 
int main()
{
    double a, b, c;
 
    cin >> a >> b >> c;
 
    double D = b * b - 4 * a * c;
 
    if (D < 0)
        cout << "D < 0. No decisions!" << endl;
    else if (D == 0)
         {
             double x1 = -b / (2 * a);
             cout << "x1 = " << x1 << endl;
         }
         else
         {
             double x1 = (-b + sqrt(D)) / (2 * a),
                    x2 = (-b - sqrt(D)) / (2 * a);
             cout << "x2 = " << x1 << endl
                  << "x2 = " << x2 << endl;
         }
 
    system("pause");
 
    return 0;
}

#include <stdio.h>
#include <math.h>
 
int sq_eq(double a, double b, double c, 
          double &x1, double &x2, double precision)
{
  if (!a)
  {
    if (b)
    { 
      x1 = -c / b, x2 = x1;
      return 1;
    }  
    return 0;
  }  
  
  double d = b * b - 4 * a * c;
  if (fabs(d) < precision) d = 0; // <--- точность учитывается здеся
  if (d < 0) return 0;
  
  x1 = (-b - sqrt(d)) / (2 * a);
  x2 = (-b + sqrt(d)) / (2 * a);
  return d ? 2 : 1;
}
 
int main()
{
  double a = 1.0, b = -0.4, c = 0.04;
  double x1, x2;
  switch (sq_eq(a, b, c, x1, x2, 5e-17))
  {
    case 2:  printf("# 2 корня: x1 = %g,  x2 = %g\n", x1, x2); break;
    case 1:  printf("# 1 корень: x1 = x2 = %g\n", x1); break;
    default: printf("# Нет корней\n");
  }
 
  return 0;
}