// quadeq.cpp: Solving a quadratic equation with complex values
#include <iostream>
#include <complex>
using namespace std;

int main() {
 complex<double> a, aa, b, c, D, sqrtD, x1, x2;
 cout << "Enter the complex numbers a (nonzero), b, and c:\n";
 cin >> a >> b >> c;
 if (a == 0.0) {
  cout << "a must be nonzero.\n";
  return 1;
 }
 D = b * b - 4.0 * a * c;
 sqrtD = sqrt(D);
 aa = 2.0 * a;
 x1 = (-b + sqrtD)/aa;
 x2 = (-b - sqrtD)/aa;
 cout << "Solution (complex):\n" << x1 << " " << x2 << endl;
 return 0;
}