ψ(y)=12y+π2coth(πy),