Комплексные числа в с++

#include <iostream>
#include <math.h>
using namespace std;
int main()
{
    float x1,x2, i;
    double a, b, c;
    double real;
    float imag;
    cout<<"The quadratic decision\nEnter value of factors\n";
    cout << "Enter 'a': \n"; 
    cin>>a;
    cout << "Enter 'b': \n";
    cin>>b;
    cout << "Enter 'c': \n"; 
    cin>>c;
    if(a==0)
    {
        cout<<"The equation isn't square"; 
    }
    else
    {
        double d=b*b-4*a*c;
        if(d<0)
        {
        cout<<"Complex roots of the equation:\n"; 
        i='i';
        real=(-b)/(2*a);
        imag=i*sqrt(fabs(d));
        x1=real+imag;
        x2=real-imag;
        cout<<"x1 = "<<x1<<"\n"<<"x2 = "<<x2;
        }
        if(d==0)
        {
            cout<<"The equation has 2 identical roots of the equation:\n"; 
            x1=-b/(2*a);
            cout<<"x = "<<x1;
        }
        if(d>0)
        {
            cout<<"Equation roots:\n"; 
            d=sqrt(d);
            x1=(-b+d)/(2*a);
            x2=(-b-d)/(2*a);
            cout<<"x1 = "<<x1<<"\n"<<"x2 = "<<x2;
        }
    }
    return 1;
}