C语言 一元二次方程解的判定及浮点数大小的比较

#include <stdio.h> #include <math.h> #define EXP 0.0000001 /*自定义精度*/ int main() { double a = 0.0, b = 0.0, c = 0.0; double x1 = 0.0, x2 = 0.0; double p = 0.0, q = 0.0; printf("输入方程系数："); scanf("%lf %lf %lf", &a, &b, &c); /*双精度浮点数一般用精度判断是否为0*/ if ((a > -EXP) && (a < EXP)) { printf("这不是一元二次方程\n"); } else { double disc = b*b - 4 * a*c; if ((disc > -EXP) && (disc < EXP)) { printf("方程有两个相等的实根:x=%lf\n",-b/(2*a)); } else if (disc > 0.0) { x1 = (-b + sqrt(disc)) / (2 * a); x2 = (-b - sqrt(disc)) / (2 * a); printf("方程有两个不等的实根:\n x1=%lf x2=%lf\n",x1,x2); } else { p = -b / (2 * a); q = sqrt(-disc) / (2 * a); printf("方程有两个共轭的复根:\n"); printf("x1 = %lf + %lfi\n", p, q); printf("x2 = %lf - %lfi\n", p, q); } } return 0; }

结果如下：

对于float,double这样的浮点数类型，在计算机中不能准确表示，也不能直接与0作比较。通常可以设置一个精度，在正负精度范围之间时，就可以认为该浮点数为0了。

例如：

#include<stdio.h> #define exp 0.000001 //定义的一个精度，可以任意定义，一般取pow(10,-6)就可以 int main() { float a = 0.0; float b = 0.0; scanf("%f %f", &a, &b); printf("a=%f b=%f\n",a,b); if (((a-b) > -exp) && (a-b) < exp) printf("a=b\n"); else if((a-b) > exp) printf("a>b\n"); else if((a-b) < exp) printf("a<b\n"); return 0; } 结果如下：

当要比较的两个数a,b相等或等于0时，必须用这种定义精度的方式。当比较的两个数a>b时，(a-b) > 0.0也是可以的；a<b同a>b类似。