do!*one!*psi(-1/2);


        -1
Psi of ---- is  - 2*log(2) - euler_gamma + 2
        2

do!*one!*psi(-3);



***** Psi undefined for nonpositive integer argument -3 

