public static double expm1(double x) { return Math.Exp(x) - 1.0; // TODO: implement exp(x)-1 for x near to zero }