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