/* Roots of ax^2 + bx + c = 0  */ 
#include <iostream> 
#include <cstdlib>
#include <cmath> 

using namespace std;

int main(int argc, char *argv[])
{   double a, b, c; 
    double discriminant, root1_real, root2_real, root_double,
                         root1_imag, root2_imag; 
     

    cout<<"Input the coefficients a, b, c, of the quadratic equation: "; 
    cin>> a >> b >> c;    

    cout<< "Coefficients are "<<a<<", "<<b<<", "<<c<<endl; 

    /* compute two roots for discriminantriminant > 0 */ 
    discriminant = pow(b,2) - 4 * a * c; 

    if (discriminant > 0) // Two real roots
    {  
        root1_real = (-b + sqrt(discriminant)) / (2 * a);  
        root2_real = (-b - sqrt(discriminant)) / (2 * a);  
          cout<<"root1 = "<<root1_real<<endl; 
          cout<<"root2 = "<<root2_real<<endl;    
    } 
    else if (discriminant == 0) //One double root
         { 
           root_double =  -b / (2 * a); 
           cout<< "Double root: "<<root_double<<endl;
         } 
         else //Imaginary roots
         {
          root1_real = -b / (2 * a);  
          root1_imag = +sqrt(-discriminant) / (2 * a);  
          root2_real = -b / (2 * a);  
          root2_imag = -sqrt(-discriminant) / (2 * a);  
          cout<<"Complex root1 = "<<root1_real<<"+I"<<root1_imag<<endl; 
          cout<<"Complex root2 = "<<root2_real<<"+I"<<root2_imag<<endl;  
          } 

        system("PAUSE");
        return 0;  
}