Currently, in sage-4.7.2: sage: integral_numerical(log(x), 0, 0) (nan, nan)
Mathematically, the integral should certainly be zero: there is a primitive function which is continuous and defined at 0. Symbolically, we can compute the integral correctly: sage: integral(log(x), (x,0,0)) 0 So I would like to add a special-case check for integral_numerical(): if the interval of integration is a point, then always return 0. I realize that this means that also the integral of 1/x from 0 to 0 would be 0, even though 1/x has no continuous primitive at 0. But according to the Lebesgue theory of integration, I think this is not even a problem. -- To post to this group, send an email to sage-devel@googlegroups.com To unsubscribe from this group, send an email to sage-devel+unsubscr...@googlegroups.com For more options, visit this group at http://groups.google.com/group/sage-devel URL: http://www.sagemath.org
