BUG: psi(x) should return -inf at x=0, not inf

See http://mathworld.wolfram.com/DigammaFunction.html

As x -> 0, psi(x) -> (-inf). Returning infinity is nonintuitive behaviour.