1058：求一元二次方程

#include <cstdio>
#include <cmath>
#define eps1 1e-10
#define eps2 1e-6
using namespace std;
int main(){
    double a,b,c,x1,x2,g; 
    scanf("%lf%lf%lf",&a,&b,&c);
    g=b*b-4*a*c;
    if(g<0&&fabs(g)>eps1)printf("No answer!");
    else if(fabs(g)<eps1){
        x1=-b/(2*a);
        if(fabs(x1)<eps2) 
            printf("x1=x2=%.5lf",0);
        else
            printf("x1=x2=%.5lf",x1);
    }
    else{
        x1=(-b+sqrt(g))/(2*a);
        x2=(-b-sqrt(g))/(2*a);
        if(fabs(x1)<eps2)x1=fabs(x1);
        if(fabs(x2)<eps2)x2=fabs(x2);
        if(x1<x2)printf("x1=%.5lf;x2=%.5lf",x1,x2);
        else printf("x1=%.5lf;x2=%.5lf",x2,x1);
    }
    return 0;//hhhh写完了233
}