解一元二次方程ax^2+bx+c=0

#include<stdio.h>//解一元二次方程ax^2+bx+c=0
#include<math.h>
int main()
{
float a,b,c,disc,x1,x2,p,q;
//输入部分 
printf("请输入a=");
scanf("%f",&a);
fflush(stdin);
printf("请输入b=");
scanf("%f",&b);
fflush(stdin);
printf("请输入c=");
scanf("%f",&c);
fflush(stdin);
printf("\na=%f,b=%f,c=%f\n",a,b,c);
disc=b*b-4*a*c;
printf("disc=%f\n",disc);

//开始运算部分
if(a==0){x1=x2=-c/b;printf("\n结果是a=0,x1=x2=%f\n\n",x1);} 
else if(disc==0){x1=x2=-b/2*a;printf("\n结果是disc=0，x1=x2=%f\n\n",x1);}
else if(disc>0)
{
p=-b/(2*a);
printf("p=%f\n",p);
q=sqrt(disc)/(2*a);
printf("q=%f\n",q);
x1=p+q;
x2=p-q;
printf("\n结果是disc>0，x1=%f,x2=%f\n",x1,x2);
}
else
{
p=-b/(2*a);
printf("p=%f\n",p);
q=sqrt(-disc)/(2*a);
printf("q=%f\n",q);
x1=p+q;
x2=p-q;
printf("\n结果是disc<0，x1=%f,x2=%f\n",x1,x2);
}

return 0;
}