Wolfram Alphaに入力すると
トリガンマ関数 ψ_1(x)=d2/dx2 lnΓ(x)
を使った解が出てくる

課金すれば導出過程を表示できる