Question

编程求解一元二次方程ax²+bx+c=0的根，其中三个系数可以是任意实数

Answer 1

// 任意实数一元二次方程的解，精度比用直接用一元二次方程通解公式高。
// 当方程无实数解，或a为0.0时返回NULL
double * Quadratic(double a, double  b, double  c, double  * x)
{
	if(a == 0.0) return NULL; // 不是二次方程
	double aa = 1.0 / a;
	double A  = 0.5 * b * aa;
	double B  = 1.0 / A;
	B = B * B * c * aa;
	if(B > 1.0) return NULL; // 无实数根
	double C = -1.0 - sqrt(1.0 - B);
	x[0] = A * C;
	x[1] = A * B / C;
	return x;
}