解一元二次方程组时，当解为两个不相等实数根时，运行时出现-1.#IND，如何解决？

Question

int main()
{

double a,b,c,d;

double x1,x2;

cin>>a>>b>>c;

d = b*b-4*a*c;

if( d<0 )
{
    cout<<"error";
}
else if ( d==0 )
{
    x1=x2=-b/(2*a);
    cout<<x1;
}
else
{
    x1=-b/(2*a)+sqrt(4*a*c-b*b)/(2*a);
    x2=-b/(2*a)-sqrt(4*a*c-b*b)/(2*a);
    cout<<x1<<" "<<x2;
}
return 0;

}

Answer 1

sqrt (4*a*c-b*b) 小于0