LegendreQ[1,x] - LegendreP[1,x] * integrate[1/((1-x^2)LegendreP[1,x]^2)]