P_n(x y - sqrt(1 - x^2) sqrt(1 - y^2) cos(α)) = P_n(x) P_n(y) + 2 sum_(k=1)^n ((P_n^k(x) P_n^k(y)) ((-1)^k cos(k α) (n - k)!))/((k + n)!)