I was able to obtain some indefinite integrals for the trigonometric function manually,
$\Gamma^{pq}_n = \int \frac {\cos^p x \sin^q x} {(1+A \cos x)^n} dx$
for $(p,q,n)\geq 0$ in recursive form. For example,
$\Gamma^{11}_n=\frac{1}{A}(\Gamma^{01}_{n-1}-\Gamma^{01}_{n}),$ for $n>2$
$\Gamma^{21}_n=\frac{1}{A^2}(\Gamma^{01}_{n}-2\Gamma^{01}_{n-1}+\Gamma^{01}_{n-2}),$ for $n>2$.
How can I generate similar recursive integrals for the trigonometric function
$\Delta^{pq}_n = \int \frac {\cos^p x \sin^q x} {(1+A \cos (x-a))^n} dx$
for $(p,q,n)\geq 0$ using Mathematica?
Simplify[\[Gamma][2, 1, n , A] == -1/ A^2 (\[Gamma][0, 1, n , A] - 2 \[Gamma][0, 1, n - 1 , A] + \[Gamma][0, 1, n - 2, A] ) /. a -> 0]
) $\endgroup$Gam[p_, q_, n_] := (Cos[x]^p Sin[x]^q)/(1 + A Cos[x])^n; Simplify[Gam[2, 1,n] == (Gam[0, 1, n] - 2 Gam[0, 1, n - 1] + Gam[0, 1, n - 2])/A^2]
This is verifying. $\endgroup$