1+2i(-qpolygamma(i pi/2,e) +qpolygamma(-i pi/2,e))