Add "DiGamma" ψ (psi) function for integer arguments

Well, probably this should bump to 0.8 with a whole new special function, but it's all good since ψ (psi) here is only implemented to the raw bare minimum needed for a mutual information estimate and would clearly suck for large arguments.