public static double log1p(double x) { return Math.Log(x + 1.0); // TODO: implement log(x+1) for x near to zero }