double a, b; ... if (a == b) /* 错! */
#include <math.h> if (fabs(a - b) <= epsilon * fabs(a))
参考资料: [Knuth, Sec. 4.2.2 pp. 217-8]。