c++求一元二次方程的解

#include<iostream>
#include<math.h>
using namespace std;
int main(){
	//求解方程ax^2+bx+c=0的解
	double a,b,c,disc,p,q,x1,x2;
	cout<<"输入a，b，c的值"<<endl;
	cin>>a;
	cin>>b;
	cin>>c;
	disc=b*b-4*a*c;
	
	p=-b/(2.0*a);
	q=sqrt(disc)/(2.0*a);
	x1=p+q;
	x2=p-q;
	if(disc<0)
	{
		cout<<"您输入的数据有误"<<endl; 
	} 
	else if(disc==0){
		cout<<"方程有两个相同的解"<<endl;
		cout<<"x1=x2="<<x1<<endl;
		 
	}
	else if(disc>0){
		cout<<"方程有两个不同的解"<<endl;
		cout<<"x1="<<x1<<" "<<"x2="<<x2<<endl; 
	}
	
	system("pause");
	return 0;
	
}