Quadratic Formula with complex numbers solutions

complex_square_roots.cpp

#include <stdio.h>
#include <complex>
#include <iostream>

int main()
{
    double a1,b1,c1;

    printf("ingrese valor de a,b,c\n");
    scanf("%lf" "%lf" "%lf", &a1, &b1, &c1);
    std::complex<double> a (a1,0);
    std::complex<double> b (b1,0);
    std::complex<double> c (c1,0);
    std::complex<double> r1 (4,0);
    std::complex<double> r2 (2,0);

    std::complex<double> z = sqrt((b*b) - (r1*a*c));
    std::complex<double> x = (-b+z)/(r2*a);
    std::complex<double> y = (-b-z)/(r2*a);

    if(x.imag() == 0 && y.imag() == 0 && x.real() != y.real()) {
        printf("x1= %lf x2=%lf -> el polinomio tiene 2 raices reales distintas",x.real(),y.real());
    }
    else if(x.imag() == 0 && y.imag() == 0 && x.real() == y.real()) {
        printf("%lf=%lf -> el polinomio tiene 2 raices reales iguales",x.real(),y.real());
    }
    else {
        printf("x1= %lfi + %lfj y x2= %lfi + %lfj -> el polinomio tiene 2 raices imaginarias conjugadas",x.real(),x.imag(),y.real(),y.imag());
    }

    return 0;
}