qpolygamma(3,1,e^pi)-lim[sum[ln(x)^4 sigma(3,n)/x^n,{n,1,infty}],x->exp(pi)]