Написать программу решения системы алгебраических уравнений

double det=a1*b2-a2*b1;
double y=(a1*c2-a2*c1)/det;
double x=(c1*b2-c2*b1)/det;

#define PRECISION    1e-12 // точность (число, меньшее по модулю, считается за 0)
#define ZERO(x)      (fabs(x) < PRECISION) // примерно равно нулю
#define EQUAL(a, b)  ZERO((a) - (b)) // примерно равно
 
// функция решения СЛАУ 2 порядка.
int foo(double a1, double b1, double c1, 
        double a2, double b2, double c2,
        double &x, double &y)
{
  if (ZERO(a1 * a2) && ZERO(b1 * b2)) return -1; // ∞ кол-во корней
  
  double D = a1 * b2 - a2 * b1;
  if (ZERO(D)) return 0; // нет корней
 
  x = (c1 * b2 - c2 * b1) / D, y = (a1 * c2 - a2 * c1) / D;
  return 1; // есть единственное решение
}

#define PRECISION    1e-12 // точность (число, меньшее по модулю, считается за 0)
#define ZERO(x)      (fabs(x) < PRECISION) // примерно равно нулю
#define EQUAL(a, b)  ZERO((a) - (b)) // примерно равно
 
// функция решения СЛАУ 2 порядка. Вернёт:
//   0, если нет корней, 
//  -1, если бесконечное множество корней, 
//  +1, если есть единственное решение.
int foo(double a1, double b1, double c1, 
        double a2, double b2, double c2,
        double &x, double &y)
{
  double D  = a1 * b2 - a2 * b1, 
         dx = c1 * b2 - c2 * b1, 
         dy = a1 * c2 - a2 * c1;
         
  if (ZERO(D)) return (ZERO(dx) && ZERO(dy)) ? -1 : 0;
 
  x = dx / D, y = dy / D;
  return 1; 
}