N[c-polygamma(0,c)/polygamma(1,c), 32], where c=b-polygamma(0,b)/polygamma(1,b), b= a-polygamma(0,a)/polygamma(1,a), a=3/2- polygamma(0,3/2)/polygamma(1,3/2)