Напишите программу полного исследования совокупности корней БИКВАДРАТНОГО уравнения

#include <cstdio>
#include <cmath>
#include <cassert>
// Квадратное уравнение
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;
}
 
// биквадратное уравнение, возвращает кол-во корней.
// массив roots должен иметь минимум 4 элемента под найденные корни.
int qq_eq(double a, double b, double c, 
          double* roots, double precision)
{
  double y[2];
  int total = 0, rc = sq_eq(a, b, c, y[0], y[1], precision);
  
  for (int idx = 0; idx != rc; idx++) 
    if (y[idx] >= 0.0) 
      roots[total] = sqrt(y[idx]), roots[total + 1] = -sqrt(y[idx]), total += 2;
  
  return total; // кол-во корней
}    
 
#define PREC 1e-9 // точность 
 
int main()
{
  double roots[4] = {0}; // массив под корни
  const double a = 1.0, b = 4.0, c = -21.0; // коэффициенты 
 
  switch (qq_eq(a, b, c, roots, PREC))
  {
    case 0: printf("# Нет корней.\n"); break;
    case 2: printf("# 2 корня: x1 = %lf, x2 = %lf\n", 
                   roots[0], roots[1]); break;
    case 4: printf("# 4 корня: x1 = %lf, x2 = %lf', x3 = %lf, x4 = %lf\n", 
                   roots[0], roots[1], roots[2], roots[3]); break;
    default: assert(false);
  }
  
  return 0;
}