public static double erfc(double x) {
		double erfc = 1.0D;
		if (x != 0.0) {
			if (x == 1.0D / 0.0D) {
				erfc = 0.0D;
			} else {
				if (x >= 0) {
					erfc = 1.0D - Stat.incompleteGamma(0.5, x * x);
				} else {
					erfc = 1.0D + Stat.incompleteGamma(0.5, x * x);
				}
			}
		}
		return erfc;
	}