求二元一次方程的解

#include<stdio.h>
#include<math.h>
int main(){
    double a,b,c,disc,x1,x2,realpart,imagpart;
    while(scanf("%lf%lf%lf",&a,&b,&c)!=EOF){
        if(fabs(a)==0) {
        x1=-c/b,x2=-c/b;
        printf("%lf%lf",x1,x2);
    }
    else{
        disc=b*b-4*c*a;
        if(fabs(disc)==0) {
            x1=-b/(2*a),x2=-b/(2*a);
            printf("%lf%lf",x1,x2);
        }
        else{
            if(disc>0){
                x1=(-b+sqrt(disc))/(2*a);
                x2=(-b-sqrt(disc))/(2*a);
                printf("%lf   %lf",x1,x2);
            }
            else{
                realpart=-b/(2*a);
                imagpart=sqrt(-disc)/(2*a);
                printf("%8.4f+%8.4fi\n",realpart,imagpart);
                printf("%8.4f-%8.4fi\n",realpart,imagpart);
                
            }
        }    
    }

    }
    return 0;
}