3.1236 \(\int \cos ^5(c+d x) \sin ^n(c+d x) (a+b \sin (c+d x)) \, dx\)

Optimal. Leaf size=123 \[ \frac {a \sin ^{n+1}(c+d x)}{d (n+1)}-\frac {2 a \sin ^{n+3}(c+d x)}{d (n+3)}+\frac {a \sin ^{n+5}(c+d x)}{d (n+5)}+\frac {b \sin ^{n+2}(c+d x)}{d (n+2)}-\frac {2 b \sin ^{n+4}(c+d x)}{d (n+4)}+\frac {b \sin ^{n+6}(c+d x)}{d (n+6)} \]

[Out]

a*sin(d*x+c)^(1+n)/d/(1+n)+b*sin(d*x+c)^(2+n)/d/(2+n)-2*a*sin(d*x+c)^(3+n)/d/(3+n)-2*b*sin(d*x+c)^(4+n)/d/(4+n
)+a*sin(d*x+c)^(5+n)/d/(5+n)+b*sin(d*x+c)^(6+n)/d/(6+n)

________________________________________________________________________________________

Rubi [A]  time = 0.14, antiderivative size = 123, normalized size of antiderivative = 1.00, number of steps used = 3, number of rules used = 2, integrand size = 27, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.074, Rules used = {2837, 766} \[ \frac {a \sin ^{n+1}(c+d x)}{d (n+1)}-\frac {2 a \sin ^{n+3}(c+d x)}{d (n+3)}+\frac {a \sin ^{n+5}(c+d x)}{d (n+5)}+\frac {b \sin ^{n+2}(c+d x)}{d (n+2)}-\frac {2 b \sin ^{n+4}(c+d x)}{d (n+4)}+\frac {b \sin ^{n+6}(c+d x)}{d (n+6)} \]

Antiderivative was successfully verified.

[In]

Int[Cos[c + d*x]^5*Sin[c + d*x]^n*(a + b*Sin[c + d*x]),x]

[Out]

(a*Sin[c + d*x]^(1 + n))/(d*(1 + n)) + (b*Sin[c + d*x]^(2 + n))/(d*(2 + n)) - (2*a*Sin[c + d*x]^(3 + n))/(d*(3
+ n)) - (2*b*Sin[c + d*x]^(4 + n))/(d*(4 + n)) + (a*Sin[c + d*x]^(5 + n))/(d*(5 + n)) + (b*Sin[c + d*x]^(6 +
n))/(d*(6 + n))

Rule 766

Int[((e_.)*(x_))^(m_.)*((f_.) + (g_.)*(x_))*((a_) + (c_.)*(x_)^2)^(p_.), x_Symbol] :> Int[ExpandIntegrand[(e*x
)^m*(f + g*x)*(a + c*x^2)^p, x], x] /; FreeQ[{a, c, e, f, g, m}, x] && IGtQ[p, 0]

Rule 2837

Int[cos[(e_.) + (f_.)*(x_)]^(p_)*((a_) + (b_.)*sin[(e_.) + (f_.)*(x_)])^(m_.)*((c_.) + (d_.)*sin[(e_.) + (f_.)
*(x_)])^(n_.), x_Symbol] :> Dist[1/(b^p*f), Subst[Int[(a + x)^m*(c + (d*x)/b)^n*(b^2 - x^2)^((p - 1)/2), x], x
, b*Sin[e + f*x]], x] /; FreeQ[{a, b, c, d, e, f, m, n}, x] && IntegerQ[(p - 1)/2] && NeQ[a^2 - b^2, 0]

Rubi steps

\begin {align*} \int \cos ^5(c+d x) \sin ^n(c+d x) (a+b \sin (c+d x)) \, dx &=\frac {\operatorname {Subst}\left (\int \left (\frac {x}{b}\right )^n (a+x) \left (b^2-x^2\right )^2 \, dx,x,b \sin (c+d x)\right )}{b^5 d}\\ &=\frac {\operatorname {Subst}\left (\int \left (a b^4 \left (\frac {x}{b}\right )^n+b^5 \left (\frac {x}{b}\right )^{1+n}-2 a b^4 \left (\frac {x}{b}\right )^{2+n}-2 b^5 \left (\frac {x}{b}\right )^{3+n}+a b^4 \left (\frac {x}{b}\right )^{4+n}+b^5 \left (\frac {x}{b}\right )^{5+n}\right ) \, dx,x,b \sin (c+d x)\right )}{b^5 d}\\ &=\frac {a \sin ^{1+n}(c+d x)}{d (1+n)}+\frac {b \sin ^{2+n}(c+d x)}{d (2+n)}-\frac {2 a \sin ^{3+n}(c+d x)}{d (3+n)}-\frac {2 b \sin ^{4+n}(c+d x)}{d (4+n)}+\frac {a \sin ^{5+n}(c+d x)}{d (5+n)}+\frac {b \sin ^{6+n}(c+d x)}{d (6+n)}\\ \end {align*}

________________________________________________________________________________________

Mathematica [A]  time = 0.21, size = 97, normalized size = 0.79 \[ \frac {\sin ^{n+1}(c+d x) \left (\frac {a \sin ^4(c+d x)}{n+5}-\frac {2 a \sin ^2(c+d x)}{n+3}+\frac {a}{n+1}+\frac {b \sin ^5(c+d x)}{n+6}-\frac {2 b \sin ^3(c+d x)}{n+4}+\frac {b \sin (c+d x)}{n+2}\right )}{d} \]

Antiderivative was successfully verified.

[In]

Integrate[Cos[c + d*x]^5*Sin[c + d*x]^n*(a + b*Sin[c + d*x]),x]

[Out]

(Sin[c + d*x]^(1 + n)*(a/(1 + n) + (b*Sin[c + d*x])/(2 + n) - (2*a*Sin[c + d*x]^2)/(3 + n) - (2*b*Sin[c + d*x]
^3)/(4 + n) + (a*Sin[c + d*x]^4)/(5 + n) + (b*Sin[c + d*x]^5)/(6 + n)))/d

________________________________________________________________________________________

fricas [B]  time = 0.84, size = 282, normalized size = 2.29 \[ -\frac {{\left ({\left (b n^{5} + 15 \, b n^{4} + 85 \, b n^{3} + 225 \, b n^{2} + 274 \, b n + 120 \, b\right )} \cos \left (d x + c\right )^{6} - {\left (b n^{5} + 11 \, b n^{4} + 41 \, b n^{3} + 61 \, b n^{2} + 30 \, b n\right )} \cos \left (d x + c\right )^{4} - 8 \, b n^{3} - 72 \, b n^{2} - 4 \, {\left (b n^{4} + 9 \, b n^{3} + 23 \, b n^{2} + 15 \, b n\right )} \cos \left (d x + c\right )^{2} - 184 \, b n - {\left ({\left (a n^{5} + 16 \, a n^{4} + 95 \, a n^{3} + 260 \, a n^{2} + 324 \, a n + 144 \, a\right )} \cos \left (d x + c\right )^{4} + 8 \, a n^{3} + 96 \, a n^{2} + 4 \, {\left (a n^{4} + 13 \, a n^{3} + 56 \, a n^{2} + 92 \, a n + 48 \, a\right )} \cos \left (d x + c\right )^{2} + 352 \, a n + 384 \, a\right )} \sin \left (d x + c\right ) - 120 \, b\right )} \sin \left (d x + c\right )^{n}}{d n^{6} + 21 \, d n^{5} + 175 \, d n^{4} + 735 \, d n^{3} + 1624 \, d n^{2} + 1764 \, d n + 720 \, d} \]

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate(cos(d*x+c)^5*sin(d*x+c)^n*(a+b*sin(d*x+c)),x, algorithm="fricas")

[Out]

-((b*n^5 + 15*b*n^4 + 85*b*n^3 + 225*b*n^2 + 274*b*n + 120*b)*cos(d*x + c)^6 - (b*n^5 + 11*b*n^4 + 41*b*n^3 +
61*b*n^2 + 30*b*n)*cos(d*x + c)^4 - 8*b*n^3 - 72*b*n^2 - 4*(b*n^4 + 9*b*n^3 + 23*b*n^2 + 15*b*n)*cos(d*x + c)^
2 - 184*b*n - ((a*n^5 + 16*a*n^4 + 95*a*n^3 + 260*a*n^2 + 324*a*n + 144*a)*cos(d*x + c)^4 + 8*a*n^3 + 96*a*n^2
+ 4*(a*n^4 + 13*a*n^3 + 56*a*n^2 + 92*a*n + 48*a)*cos(d*x + c)^2 + 352*a*n + 384*a)*sin(d*x + c) - 120*b)*sin
(d*x + c)^n/(d*n^6 + 21*d*n^5 + 175*d*n^4 + 735*d*n^3 + 1624*d*n^2 + 1764*d*n + 720*d)

________________________________________________________________________________________

giac [F(-1)]  time = 0.00, size = 0, normalized size = 0.00 \[ \text {Timed out} \]

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate(cos(d*x+c)^5*sin(d*x+c)^n*(a+b*sin(d*x+c)),x, algorithm="giac")

[Out]

Timed out

________________________________________________________________________________________

maple [F]  time = 11.38, size = 0, normalized size = 0.00 \[ \int \left (\cos ^{5}\left (d x +c \right )\right ) \left (\sin ^{n}\left (d x +c \right )\right ) \left (a +b \sin \left (d x +c \right )\right )\, dx \]

Verification of antiderivative is not currently implemented for this CAS.

[In]

int(cos(d*x+c)^5*sin(d*x+c)^n*(a+b*sin(d*x+c)),x)

[Out]

int(cos(d*x+c)^5*sin(d*x+c)^n*(a+b*sin(d*x+c)),x)

________________________________________________________________________________________

maxima [A]  time = 0.33, size = 109, normalized size = 0.89 \[ \frac {\frac {b \sin \left (d x + c\right )^{n + 6}}{n + 6} + \frac {a \sin \left (d x + c\right )^{n + 5}}{n + 5} - \frac {2 \, b \sin \left (d x + c\right )^{n + 4}}{n + 4} - \frac {2 \, a \sin \left (d x + c\right )^{n + 3}}{n + 3} + \frac {b \sin \left (d x + c\right )^{n + 2}}{n + 2} + \frac {a \sin \left (d x + c\right )^{n + 1}}{n + 1}}{d} \]

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate(cos(d*x+c)^5*sin(d*x+c)^n*(a+b*sin(d*x+c)),x, algorithm="maxima")

[Out]

(b*sin(d*x + c)^(n + 6)/(n + 6) + a*sin(d*x + c)^(n + 5)/(n + 5) - 2*b*sin(d*x + c)^(n + 4)/(n + 4) - 2*a*sin(
d*x + c)^(n + 3)/(n + 3) + b*sin(d*x + c)^(n + 2)/(n + 2) + a*sin(d*x + c)^(n + 1)/(n + 1))/d

________________________________________________________________________________________

mupad [B]  time = 16.34, size = 550, normalized size = 4.47 \[ \frac {b\,{\sin \left (c+d\,x\right )}^n\,\left (n^5+23\,n^4+237\,n^3+1129\,n^2+2234\,n+1320\right )}{16\,d\,\left (n^6+21\,n^5+175\,n^4+735\,n^3+1624\,n^2+1764\,n+720\right )}-\frac {b\,{\sin \left (c+d\,x\right )}^n\,\cos \left (6\,c+6\,d\,x\right )\,\left (n^5+15\,n^4+85\,n^3+225\,n^2+274\,n+120\right )}{32\,d\,\left (n^6+21\,n^5+175\,n^4+735\,n^3+1624\,n^2+1764\,n+720\right )}-\frac {b\,{\sin \left (c+d\,x\right )}^n\,\cos \left (4\,c+4\,d\,x\right )\,\left (n^5+23\,n^4+173\,n^3+553\,n^2+762\,n+360\right )}{16\,d\,\left (n^6+21\,n^5+175\,n^4+735\,n^3+1624\,n^2+1764\,n+720\right )}-\frac {a\,\sin \left (c+d\,x\right )\,{\sin \left (c+d\,x\right )}^n\,\left (n^5\,1{}\mathrm {i}+n^4\,24{}\mathrm {i}+n^3\,263{}\mathrm {i}+n^2\,1476{}\mathrm {i}+n\,3876{}\mathrm {i}+3600{}\mathrm {i}\right )\,1{}\mathrm {i}}{8\,d\,\left (n^6+21\,n^5+175\,n^4+735\,n^3+1624\,n^2+1764\,n+720\right )}-\frac {b\,{\sin \left (c+d\,x\right )}^n\,\cos \left (2\,c+2\,d\,x\right )\,\left (-n^5-15\,n^4+43\,n^3+927\,n^2+2670\,n+1800\right )}{32\,d\,\left (n^6+21\,n^5+175\,n^4+735\,n^3+1624\,n^2+1764\,n+720\right )}-\frac {a\,{\sin \left (c+d\,x\right )}^n\,\sin \left (5\,c+5\,d\,x\right )\,\left (n^5\,1{}\mathrm {i}+n^4\,16{}\mathrm {i}+n^3\,95{}\mathrm {i}+n^2\,260{}\mathrm {i}+n\,324{}\mathrm {i}+144{}\mathrm {i}\right )\,1{}\mathrm {i}}{16\,d\,\left (n^6+21\,n^5+175\,n^4+735\,n^3+1624\,n^2+1764\,n+720\right )}-\frac {a\,{\sin \left (c+d\,x\right )}^n\,\sin \left (3\,c+3\,d\,x\right )\,\left (n^5\,3{}\mathrm {i}+n^4\,64{}\mathrm {i}+n^3\,493{}\mathrm {i}+n^2\,1676{}\mathrm {i}+n\,2444{}\mathrm {i}+1200{}\mathrm {i}\right )\,1{}\mathrm {i}}{16\,d\,\left (n^6+21\,n^5+175\,n^4+735\,n^3+1624\,n^2+1764\,n+720\right )} \]

Verification of antiderivative is not currently implemented for this CAS.

[In]

int(cos(c + d*x)^5*sin(c + d*x)^n*(a + b*sin(c + d*x)),x)

[Out]

(b*sin(c + d*x)^n*(2234*n + 1129*n^2 + 237*n^3 + 23*n^4 + n^5 + 1320))/(16*d*(1764*n + 1624*n^2 + 735*n^3 + 17
5*n^4 + 21*n^5 + n^6 + 720)) - (b*sin(c + d*x)^n*cos(6*c + 6*d*x)*(274*n + 225*n^2 + 85*n^3 + 15*n^4 + n^5 + 1
20))/(32*d*(1764*n + 1624*n^2 + 735*n^3 + 175*n^4 + 21*n^5 + n^6 + 720)) - (b*sin(c + d*x)^n*cos(4*c + 4*d*x)*
(762*n + 553*n^2 + 173*n^3 + 23*n^4 + n^5 + 360))/(16*d*(1764*n + 1624*n^2 + 735*n^3 + 175*n^4 + 21*n^5 + n^6
+ 720)) - (a*sin(c + d*x)*sin(c + d*x)^n*(n*3876i + n^2*1476i + n^3*263i + n^4*24i + n^5*1i + 3600i)*1i)/(8*d*
(1764*n + 1624*n^2 + 735*n^3 + 175*n^4 + 21*n^5 + n^6 + 720)) - (b*sin(c + d*x)^n*cos(2*c + 2*d*x)*(2670*n + 9
27*n^2 + 43*n^3 - 15*n^4 - n^5 + 1800))/(32*d*(1764*n + 1624*n^2 + 735*n^3 + 175*n^4 + 21*n^5 + n^6 + 720)) -
(a*sin(c + d*x)^n*sin(5*c + 5*d*x)*(n*324i + n^2*260i + n^3*95i + n^4*16i + n^5*1i + 144i)*1i)/(16*d*(1764*n +
1624*n^2 + 735*n^3 + 175*n^4 + 21*n^5 + n^6 + 720)) - (a*sin(c + d*x)^n*sin(3*c + 3*d*x)*(n*2444i + n^2*1676i
+ n^3*493i + n^4*64i + n^5*3i + 1200i)*1i)/(16*d*(1764*n + 1624*n^2 + 735*n^3 + 175*n^4 + 21*n^5 + n^6 + 720)
)

________________________________________________________________________________________

sympy [A]  time = 100.87, size = 8675, normalized size = 70.53 \[ \text {result too large to display} \]

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate(cos(d*x+c)**5*sin(d*x+c)**n*(a+b*sin(d*x+c)),x)

[Out]

Piecewise((x*(a + b*sin(c))*sin(c)**n*cos(c)**5, Eq(d, 0)), (-8*a/(15*d*sin(c + d*x)) + 4*a*cos(c + d*x)**2/(1
5*d*sin(c + d*x)**3) - a*cos(c + d*x)**4/(5*d*sin(c + d*x)**5) + b*log(sin(c + d*x))/d + b*cos(c + d*x)**2/(2*
d*sin(c + d*x)**2) - b*cos(c + d*x)**4/(4*d*sin(c + d*x)**4), Eq(n, -6)), (a*log(sin(c + d*x))/d + a*cos(c + d
*x)**2/(2*d*sin(c + d*x)**2) - a*cos(c + d*x)**4/(4*d*sin(c + d*x)**4) + 8*b*sin(c + d*x)/(3*d) + 4*b*cos(c +
d*x)**2/(3*d*sin(c + d*x)) - b*cos(c + d*x)**4/(3*d*sin(c + d*x)**3), Eq(n, -5)), (-a*tan(c/2 + d*x/2)**10/(24
*d*tan(c/2 + d*x/2)**7 + 48*d*tan(c/2 + d*x/2)**5 + 24*d*tan(c/2 + d*x/2)**3) + 19*a*tan(c/2 + d*x/2)**8/(24*d
*tan(c/2 + d*x/2)**7 + 48*d*tan(c/2 + d*x/2)**5 + 24*d*tan(c/2 + d*x/2)**3) + 110*a*tan(c/2 + d*x/2)**6/(24*d*
tan(c/2 + d*x/2)**7 + 48*d*tan(c/2 + d*x/2)**5 + 24*d*tan(c/2 + d*x/2)**3) + 110*a*tan(c/2 + d*x/2)**4/(24*d*t
an(c/2 + d*x/2)**7 + 48*d*tan(c/2 + d*x/2)**5 + 24*d*tan(c/2 + d*x/2)**3) + 19*a*tan(c/2 + d*x/2)**2/(24*d*tan
(c/2 + d*x/2)**7 + 48*d*tan(c/2 + d*x/2)**5 + 24*d*tan(c/2 + d*x/2)**3) - a/(24*d*tan(c/2 + d*x/2)**7 + 48*d*t
an(c/2 + d*x/2)**5 + 24*d*tan(c/2 + d*x/2)**3) + 48*b*log(tan(c/2 + d*x/2)**2 + 1)*tan(c/2 + d*x/2)**7/(24*d*t
an(c/2 + d*x/2)**7 + 48*d*tan(c/2 + d*x/2)**5 + 24*d*tan(c/2 + d*x/2)**3) + 96*b*log(tan(c/2 + d*x/2)**2 + 1)*
tan(c/2 + d*x/2)**5/(24*d*tan(c/2 + d*x/2)**7 + 48*d*tan(c/2 + d*x/2)**5 + 24*d*tan(c/2 + d*x/2)**3) + 48*b*lo
g(tan(c/2 + d*x/2)**2 + 1)*tan(c/2 + d*x/2)**3/(24*d*tan(c/2 + d*x/2)**7 + 48*d*tan(c/2 + d*x/2)**5 + 24*d*tan
(c/2 + d*x/2)**3) - 48*b*log(tan(c/2 + d*x/2))*tan(c/2 + d*x/2)**7/(24*d*tan(c/2 + d*x/2)**7 + 48*d*tan(c/2 +
d*x/2)**5 + 24*d*tan(c/2 + d*x/2)**3) - 96*b*log(tan(c/2 + d*x/2))*tan(c/2 + d*x/2)**5/(24*d*tan(c/2 + d*x/2)*
*7 + 48*d*tan(c/2 + d*x/2)**5 + 24*d*tan(c/2 + d*x/2)**3) - 48*b*log(tan(c/2 + d*x/2))*tan(c/2 + d*x/2)**3/(24
*d*tan(c/2 + d*x/2)**7 + 48*d*tan(c/2 + d*x/2)**5 + 24*d*tan(c/2 + d*x/2)**3) - 3*b*tan(c/2 + d*x/2)**9/(24*d*
tan(c/2 + d*x/2)**7 + 48*d*tan(c/2 + d*x/2)**5 + 24*d*tan(c/2 + d*x/2)**3) + 54*b*tan(c/2 + d*x/2)**5/(24*d*ta
n(c/2 + d*x/2)**7 + 48*d*tan(c/2 + d*x/2)**5 + 24*d*tan(c/2 + d*x/2)**3) - 3*b*tan(c/2 + d*x/2)/(24*d*tan(c/2
+ d*x/2)**7 + 48*d*tan(c/2 + d*x/2)**5 + 24*d*tan(c/2 + d*x/2)**3), Eq(n, -4)), (48*a*log(tan(c/2 + d*x/2)**2
+ 1)*tan(c/2 + d*x/2)**8/(24*d*tan(c/2 + d*x/2)**8 + 72*d*tan(c/2 + d*x/2)**6 + 72*d*tan(c/2 + d*x/2)**4 + 24*
d*tan(c/2 + d*x/2)**2) + 144*a*log(tan(c/2 + d*x/2)**2 + 1)*tan(c/2 + d*x/2)**6/(24*d*tan(c/2 + d*x/2)**8 + 72
*d*tan(c/2 + d*x/2)**6 + 72*d*tan(c/2 + d*x/2)**4 + 24*d*tan(c/2 + d*x/2)**2) + 144*a*log(tan(c/2 + d*x/2)**2
+ 1)*tan(c/2 + d*x/2)**4/(24*d*tan(c/2 + d*x/2)**8 + 72*d*tan(c/2 + d*x/2)**6 + 72*d*tan(c/2 + d*x/2)**4 + 24*
d*tan(c/2 + d*x/2)**2) + 48*a*log(tan(c/2 + d*x/2)**2 + 1)*tan(c/2 + d*x/2)**2/(24*d*tan(c/2 + d*x/2)**8 + 72*
d*tan(c/2 + d*x/2)**6 + 72*d*tan(c/2 + d*x/2)**4 + 24*d*tan(c/2 + d*x/2)**2) - 48*a*log(tan(c/2 + d*x/2))*tan(
c/2 + d*x/2)**8/(24*d*tan(c/2 + d*x/2)**8 + 72*d*tan(c/2 + d*x/2)**6 + 72*d*tan(c/2 + d*x/2)**4 + 24*d*tan(c/2
+ d*x/2)**2) - 144*a*log(tan(c/2 + d*x/2))*tan(c/2 + d*x/2)**6/(24*d*tan(c/2 + d*x/2)**8 + 72*d*tan(c/2 + d*x
/2)**6 + 72*d*tan(c/2 + d*x/2)**4 + 24*d*tan(c/2 + d*x/2)**2) - 144*a*log(tan(c/2 + d*x/2))*tan(c/2 + d*x/2)**
4/(24*d*tan(c/2 + d*x/2)**8 + 72*d*tan(c/2 + d*x/2)**6 + 72*d*tan(c/2 + d*x/2)**4 + 24*d*tan(c/2 + d*x/2)**2)
- 48*a*log(tan(c/2 + d*x/2))*tan(c/2 + d*x/2)**2/(24*d*tan(c/2 + d*x/2)**8 + 72*d*tan(c/2 + d*x/2)**6 + 72*d*t
an(c/2 + d*x/2)**4 + 24*d*tan(c/2 + d*x/2)**2) - 3*a*tan(c/2 + d*x/2)**10/(24*d*tan(c/2 + d*x/2)**8 + 72*d*tan
(c/2 + d*x/2)**6 + 72*d*tan(c/2 + d*x/2)**4 + 24*d*tan(c/2 + d*x/2)**2) + 63*a*tan(c/2 + d*x/2)**6/(24*d*tan(c
/2 + d*x/2)**8 + 72*d*tan(c/2 + d*x/2)**6 + 72*d*tan(c/2 + d*x/2)**4 + 24*d*tan(c/2 + d*x/2)**2) + 63*a*tan(c/
2 + d*x/2)**4/(24*d*tan(c/2 + d*x/2)**8 + 72*d*tan(c/2 + d*x/2)**6 + 72*d*tan(c/2 + d*x/2)**4 + 24*d*tan(c/2 +
d*x/2)**2) - 3*a/(24*d*tan(c/2 + d*x/2)**8 + 72*d*tan(c/2 + d*x/2)**6 + 72*d*tan(c/2 + d*x/2)**4 + 24*d*tan(c
/2 + d*x/2)**2) - 12*b*tan(c/2 + d*x/2)**9/(24*d*tan(c/2 + d*x/2)**8 + 72*d*tan(c/2 + d*x/2)**6 + 72*d*tan(c/2
+ d*x/2)**4 + 24*d*tan(c/2 + d*x/2)**2) - 144*b*tan(c/2 + d*x/2)**7/(24*d*tan(c/2 + d*x/2)**8 + 72*d*tan(c/2
+ d*x/2)**6 + 72*d*tan(c/2 + d*x/2)**4 + 24*d*tan(c/2 + d*x/2)**2) - 200*b*tan(c/2 + d*x/2)**5/(24*d*tan(c/2 +
d*x/2)**8 + 72*d*tan(c/2 + d*x/2)**6 + 72*d*tan(c/2 + d*x/2)**4 + 24*d*tan(c/2 + d*x/2)**2) - 144*b*tan(c/2 +
d*x/2)**3/(24*d*tan(c/2 + d*x/2)**8 + 72*d*tan(c/2 + d*x/2)**6 + 72*d*tan(c/2 + d*x/2)**4 + 24*d*tan(c/2 + d*
x/2)**2) - 12*b*tan(c/2 + d*x/2)/(24*d*tan(c/2 + d*x/2)**8 + 72*d*tan(c/2 + d*x/2)**6 + 72*d*tan(c/2 + d*x/2)*
*4 + 24*d*tan(c/2 + d*x/2)**2), Eq(n, -3)), (-3*a*tan(c/2 + d*x/2)**10/(6*d*tan(c/2 + d*x/2)**9 + 24*d*tan(c/2
+ d*x/2)**7 + 36*d*tan(c/2 + d*x/2)**5 + 24*d*tan(c/2 + d*x/2)**3 + 6*d*tan(c/2 + d*x/2)) - 39*a*tan(c/2 + d*
x/2)**8/(6*d*tan(c/2 + d*x/2)**9 + 24*d*tan(c/2 + d*x/2)**7 + 36*d*tan(c/2 + d*x/2)**5 + 24*d*tan(c/2 + d*x/2)
**3 + 6*d*tan(c/2 + d*x/2)) - 86*a*tan(c/2 + d*x/2)**6/(6*d*tan(c/2 + d*x/2)**9 + 24*d*tan(c/2 + d*x/2)**7 + 3
6*d*tan(c/2 + d*x/2)**5 + 24*d*tan(c/2 + d*x/2)**3 + 6*d*tan(c/2 + d*x/2)) - 86*a*tan(c/2 + d*x/2)**4/(6*d*tan
(c/2 + d*x/2)**9 + 24*d*tan(c/2 + d*x/2)**7 + 36*d*tan(c/2 + d*x/2)**5 + 24*d*tan(c/2 + d*x/2)**3 + 6*d*tan(c/
2 + d*x/2)) - 39*a*tan(c/2 + d*x/2)**2/(6*d*tan(c/2 + d*x/2)**9 + 24*d*tan(c/2 + d*x/2)**7 + 36*d*tan(c/2 + d*
x/2)**5 + 24*d*tan(c/2 + d*x/2)**3 + 6*d*tan(c/2 + d*x/2)) - 3*a/(6*d*tan(c/2 + d*x/2)**9 + 24*d*tan(c/2 + d*x
/2)**7 + 36*d*tan(c/2 + d*x/2)**5 + 24*d*tan(c/2 + d*x/2)**3 + 6*d*tan(c/2 + d*x/2)) - 6*b*log(tan(c/2 + d*x/2
)**2 + 1)*tan(c/2 + d*x/2)**9/(6*d*tan(c/2 + d*x/2)**9 + 24*d*tan(c/2 + d*x/2)**7 + 36*d*tan(c/2 + d*x/2)**5 +
24*d*tan(c/2 + d*x/2)**3 + 6*d*tan(c/2 + d*x/2)) - 24*b*log(tan(c/2 + d*x/2)**2 + 1)*tan(c/2 + d*x/2)**7/(6*d
*tan(c/2 + d*x/2)**9 + 24*d*tan(c/2 + d*x/2)**7 + 36*d*tan(c/2 + d*x/2)**5 + 24*d*tan(c/2 + d*x/2)**3 + 6*d*ta
n(c/2 + d*x/2)) - 36*b*log(tan(c/2 + d*x/2)**2 + 1)*tan(c/2 + d*x/2)**5/(6*d*tan(c/2 + d*x/2)**9 + 24*d*tan(c/
2 + d*x/2)**7 + 36*d*tan(c/2 + d*x/2)**5 + 24*d*tan(c/2 + d*x/2)**3 + 6*d*tan(c/2 + d*x/2)) - 24*b*log(tan(c/2
+ d*x/2)**2 + 1)*tan(c/2 + d*x/2)**3/(6*d*tan(c/2 + d*x/2)**9 + 24*d*tan(c/2 + d*x/2)**7 + 36*d*tan(c/2 + d*x
/2)**5 + 24*d*tan(c/2 + d*x/2)**3 + 6*d*tan(c/2 + d*x/2)) - 6*b*log(tan(c/2 + d*x/2)**2 + 1)*tan(c/2 + d*x/2)/
(6*d*tan(c/2 + d*x/2)**9 + 24*d*tan(c/2 + d*x/2)**7 + 36*d*tan(c/2 + d*x/2)**5 + 24*d*tan(c/2 + d*x/2)**3 + 6*
d*tan(c/2 + d*x/2)) + 6*b*log(tan(c/2 + d*x/2))*tan(c/2 + d*x/2)**9/(6*d*tan(c/2 + d*x/2)**9 + 24*d*tan(c/2 +
d*x/2)**7 + 36*d*tan(c/2 + d*x/2)**5 + 24*d*tan(c/2 + d*x/2)**3 + 6*d*tan(c/2 + d*x/2)) + 24*b*log(tan(c/2 + d
*x/2))*tan(c/2 + d*x/2)**7/(6*d*tan(c/2 + d*x/2)**9 + 24*d*tan(c/2 + d*x/2)**7 + 36*d*tan(c/2 + d*x/2)**5 + 24
*d*tan(c/2 + d*x/2)**3 + 6*d*tan(c/2 + d*x/2)) + 36*b*log(tan(c/2 + d*x/2))*tan(c/2 + d*x/2)**5/(6*d*tan(c/2 +
d*x/2)**9 + 24*d*tan(c/2 + d*x/2)**7 + 36*d*tan(c/2 + d*x/2)**5 + 24*d*tan(c/2 + d*x/2)**3 + 6*d*tan(c/2 + d*
x/2)) + 24*b*log(tan(c/2 + d*x/2))*tan(c/2 + d*x/2)**3/(6*d*tan(c/2 + d*x/2)**9 + 24*d*tan(c/2 + d*x/2)**7 + 3
6*d*tan(c/2 + d*x/2)**5 + 24*d*tan(c/2 + d*x/2)**3 + 6*d*tan(c/2 + d*x/2)) + 6*b*log(tan(c/2 + d*x/2))*tan(c/2
+ d*x/2)/(6*d*tan(c/2 + d*x/2)**9 + 24*d*tan(c/2 + d*x/2)**7 + 36*d*tan(c/2 + d*x/2)**5 + 24*d*tan(c/2 + d*x/
2)**3 + 6*d*tan(c/2 + d*x/2)) - 24*b*tan(c/2 + d*x/2)**7/(6*d*tan(c/2 + d*x/2)**9 + 24*d*tan(c/2 + d*x/2)**7 +
36*d*tan(c/2 + d*x/2)**5 + 24*d*tan(c/2 + d*x/2)**3 + 6*d*tan(c/2 + d*x/2)) - 24*b*tan(c/2 + d*x/2)**5/(6*d*t
an(c/2 + d*x/2)**9 + 24*d*tan(c/2 + d*x/2)**7 + 36*d*tan(c/2 + d*x/2)**5 + 24*d*tan(c/2 + d*x/2)**3 + 6*d*tan(
c/2 + d*x/2)) - 24*b*tan(c/2 + d*x/2)**3/(6*d*tan(c/2 + d*x/2)**9 + 24*d*tan(c/2 + d*x/2)**7 + 36*d*tan(c/2 +
d*x/2)**5 + 24*d*tan(c/2 + d*x/2)**3 + 6*d*tan(c/2 + d*x/2)), Eq(n, -2)), (-15*a*log(tan(c/2 + d*x/2)**2 + 1)*
tan(c/2 + d*x/2)**10/(15*d*tan(c/2 + d*x/2)**10 + 75*d*tan(c/2 + d*x/2)**8 + 150*d*tan(c/2 + d*x/2)**6 + 150*d
*tan(c/2 + d*x/2)**4 + 75*d*tan(c/2 + d*x/2)**2 + 15*d) - 75*a*log(tan(c/2 + d*x/2)**2 + 1)*tan(c/2 + d*x/2)**
8/(15*d*tan(c/2 + d*x/2)**10 + 75*d*tan(c/2 + d*x/2)**8 + 150*d*tan(c/2 + d*x/2)**6 + 150*d*tan(c/2 + d*x/2)**
4 + 75*d*tan(c/2 + d*x/2)**2 + 15*d) - 150*a*log(tan(c/2 + d*x/2)**2 + 1)*tan(c/2 + d*x/2)**6/(15*d*tan(c/2 +
d*x/2)**10 + 75*d*tan(c/2 + d*x/2)**8 + 150*d*tan(c/2 + d*x/2)**6 + 150*d*tan(c/2 + d*x/2)**4 + 75*d*tan(c/2 +
d*x/2)**2 + 15*d) - 150*a*log(tan(c/2 + d*x/2)**2 + 1)*tan(c/2 + d*x/2)**4/(15*d*tan(c/2 + d*x/2)**10 + 75*d*
tan(c/2 + d*x/2)**8 + 150*d*tan(c/2 + d*x/2)**6 + 150*d*tan(c/2 + d*x/2)**4 + 75*d*tan(c/2 + d*x/2)**2 + 15*d)
- 75*a*log(tan(c/2 + d*x/2)**2 + 1)*tan(c/2 + d*x/2)**2/(15*d*tan(c/2 + d*x/2)**10 + 75*d*tan(c/2 + d*x/2)**8
+ 150*d*tan(c/2 + d*x/2)**6 + 150*d*tan(c/2 + d*x/2)**4 + 75*d*tan(c/2 + d*x/2)**2 + 15*d) - 15*a*log(tan(c/2
+ d*x/2)**2 + 1)/(15*d*tan(c/2 + d*x/2)**10 + 75*d*tan(c/2 + d*x/2)**8 + 150*d*tan(c/2 + d*x/2)**6 + 150*d*ta
n(c/2 + d*x/2)**4 + 75*d*tan(c/2 + d*x/2)**2 + 15*d) + 15*a*log(tan(c/2 + d*x/2))*tan(c/2 + d*x/2)**10/(15*d*t
an(c/2 + d*x/2)**10 + 75*d*tan(c/2 + d*x/2)**8 + 150*d*tan(c/2 + d*x/2)**6 + 150*d*tan(c/2 + d*x/2)**4 + 75*d*
tan(c/2 + d*x/2)**2 + 15*d) + 75*a*log(tan(c/2 + d*x/2))*tan(c/2 + d*x/2)**8/(15*d*tan(c/2 + d*x/2)**10 + 75*d
*tan(c/2 + d*x/2)**8 + 150*d*tan(c/2 + d*x/2)**6 + 150*d*tan(c/2 + d*x/2)**4 + 75*d*tan(c/2 + d*x/2)**2 + 15*d
) + 150*a*log(tan(c/2 + d*x/2))*tan(c/2 + d*x/2)**6/(15*d*tan(c/2 + d*x/2)**10 + 75*d*tan(c/2 + d*x/2)**8 + 15
0*d*tan(c/2 + d*x/2)**6 + 150*d*tan(c/2 + d*x/2)**4 + 75*d*tan(c/2 + d*x/2)**2 + 15*d) + 150*a*log(tan(c/2 + d
*x/2))*tan(c/2 + d*x/2)**4/(15*d*tan(c/2 + d*x/2)**10 + 75*d*tan(c/2 + d*x/2)**8 + 150*d*tan(c/2 + d*x/2)**6 +
150*d*tan(c/2 + d*x/2)**4 + 75*d*tan(c/2 + d*x/2)**2 + 15*d) + 75*a*log(tan(c/2 + d*x/2))*tan(c/2 + d*x/2)**2
/(15*d*tan(c/2 + d*x/2)**10 + 75*d*tan(c/2 + d*x/2)**8 + 150*d*tan(c/2 + d*x/2)**6 + 150*d*tan(c/2 + d*x/2)**4
+ 75*d*tan(c/2 + d*x/2)**2 + 15*d) + 15*a*log(tan(c/2 + d*x/2))/(15*d*tan(c/2 + d*x/2)**10 + 75*d*tan(c/2 + d
*x/2)**8 + 150*d*tan(c/2 + d*x/2)**6 + 150*d*tan(c/2 + d*x/2)**4 + 75*d*tan(c/2 + d*x/2)**2 + 15*d) - 60*a*tan
(c/2 + d*x/2)**8/(15*d*tan(c/2 + d*x/2)**10 + 75*d*tan(c/2 + d*x/2)**8 + 150*d*tan(c/2 + d*x/2)**6 + 150*d*tan
(c/2 + d*x/2)**4 + 75*d*tan(c/2 + d*x/2)**2 + 15*d) - 120*a*tan(c/2 + d*x/2)**6/(15*d*tan(c/2 + d*x/2)**10 + 7
5*d*tan(c/2 + d*x/2)**8 + 150*d*tan(c/2 + d*x/2)**6 + 150*d*tan(c/2 + d*x/2)**4 + 75*d*tan(c/2 + d*x/2)**2 + 1
5*d) - 120*a*tan(c/2 + d*x/2)**4/(15*d*tan(c/2 + d*x/2)**10 + 75*d*tan(c/2 + d*x/2)**8 + 150*d*tan(c/2 + d*x/2
)**6 + 150*d*tan(c/2 + d*x/2)**4 + 75*d*tan(c/2 + d*x/2)**2 + 15*d) - 60*a*tan(c/2 + d*x/2)**2/(15*d*tan(c/2 +
d*x/2)**10 + 75*d*tan(c/2 + d*x/2)**8 + 150*d*tan(c/2 + d*x/2)**6 + 150*d*tan(c/2 + d*x/2)**4 + 75*d*tan(c/2
+ d*x/2)**2 + 15*d) + 30*b*tan(c/2 + d*x/2)**9/(15*d*tan(c/2 + d*x/2)**10 + 75*d*tan(c/2 + d*x/2)**8 + 150*d*t
an(c/2 + d*x/2)**6 + 150*d*tan(c/2 + d*x/2)**4 + 75*d*tan(c/2 + d*x/2)**2 + 15*d) + 40*b*tan(c/2 + d*x/2)**7/(
15*d*tan(c/2 + d*x/2)**10 + 75*d*tan(c/2 + d*x/2)**8 + 150*d*tan(c/2 + d*x/2)**6 + 150*d*tan(c/2 + d*x/2)**4 +
75*d*tan(c/2 + d*x/2)**2 + 15*d) + 116*b*tan(c/2 + d*x/2)**5/(15*d*tan(c/2 + d*x/2)**10 + 75*d*tan(c/2 + d*x/
2)**8 + 150*d*tan(c/2 + d*x/2)**6 + 150*d*tan(c/2 + d*x/2)**4 + 75*d*tan(c/2 + d*x/2)**2 + 15*d) + 40*b*tan(c/
2 + d*x/2)**3/(15*d*tan(c/2 + d*x/2)**10 + 75*d*tan(c/2 + d*x/2)**8 + 150*d*tan(c/2 + d*x/2)**6 + 150*d*tan(c/
2 + d*x/2)**4 + 75*d*tan(c/2 + d*x/2)**2 + 15*d) + 30*b*tan(c/2 + d*x/2)/(15*d*tan(c/2 + d*x/2)**10 + 75*d*tan
(c/2 + d*x/2)**8 + 150*d*tan(c/2 + d*x/2)**6 + 150*d*tan(c/2 + d*x/2)**4 + 75*d*tan(c/2 + d*x/2)**2 + 15*d), E
q(n, -1)), (a*n**5*sin(c + d*x)*sin(c + d*x)**n*cos(c + d*x)**4/(d*n**6 + 21*d*n**5 + 175*d*n**4 + 735*d*n**3
+ 1624*d*n**2 + 1764*d*n + 720*d) + 4*a*n**4*sin(c + d*x)**3*sin(c + d*x)**n*cos(c + d*x)**2/(d*n**6 + 21*d*n*
*5 + 175*d*n**4 + 735*d*n**3 + 1624*d*n**2 + 1764*d*n + 720*d) + 20*a*n**4*sin(c + d*x)*sin(c + d*x)**n*cos(c
+ d*x)**4/(d*n**6 + 21*d*n**5 + 175*d*n**4 + 735*d*n**3 + 1624*d*n**2 + 1764*d*n + 720*d) + 8*a*n**3*sin(c + d
*x)**5*sin(c + d*x)**n/(d*n**6 + 21*d*n**5 + 175*d*n**4 + 735*d*n**3 + 1624*d*n**2 + 1764*d*n + 720*d) + 68*a*
n**3*sin(c + d*x)**3*sin(c + d*x)**n*cos(c + d*x)**2/(d*n**6 + 21*d*n**5 + 175*d*n**4 + 735*d*n**3 + 1624*d*n*
*2 + 1764*d*n + 720*d) + 155*a*n**3*sin(c + d*x)*sin(c + d*x)**n*cos(c + d*x)**4/(d*n**6 + 21*d*n**5 + 175*d*n
**4 + 735*d*n**3 + 1624*d*n**2 + 1764*d*n + 720*d) + 96*a*n**2*sin(c + d*x)**5*sin(c + d*x)**n/(d*n**6 + 21*d*
n**5 + 175*d*n**4 + 735*d*n**3 + 1624*d*n**2 + 1764*d*n + 720*d) + 416*a*n**2*sin(c + d*x)**3*sin(c + d*x)**n*
cos(c + d*x)**2/(d*n**6 + 21*d*n**5 + 175*d*n**4 + 735*d*n**3 + 1624*d*n**2 + 1764*d*n + 720*d) + 580*a*n**2*s
in(c + d*x)*sin(c + d*x)**n*cos(c + d*x)**4/(d*n**6 + 21*d*n**5 + 175*d*n**4 + 735*d*n**3 + 1624*d*n**2 + 1764
*d*n + 720*d) + 352*a*n*sin(c + d*x)**5*sin(c + d*x)**n/(d*n**6 + 21*d*n**5 + 175*d*n**4 + 735*d*n**3 + 1624*d
*n**2 + 1764*d*n + 720*d) + 1072*a*n*sin(c + d*x)**3*sin(c + d*x)**n*cos(c + d*x)**2/(d*n**6 + 21*d*n**5 + 175
*d*n**4 + 735*d*n**3 + 1624*d*n**2 + 1764*d*n + 720*d) + 1044*a*n*sin(c + d*x)*sin(c + d*x)**n*cos(c + d*x)**4
/(d*n**6 + 21*d*n**5 + 175*d*n**4 + 735*d*n**3 + 1624*d*n**2 + 1764*d*n + 720*d) + 384*a*sin(c + d*x)**5*sin(c
+ d*x)**n/(d*n**6 + 21*d*n**5 + 175*d*n**4 + 735*d*n**3 + 1624*d*n**2 + 1764*d*n + 720*d) + 960*a*sin(c + d*x
)**3*sin(c + d*x)**n*cos(c + d*x)**2/(d*n**6 + 21*d*n**5 + 175*d*n**4 + 735*d*n**3 + 1624*d*n**2 + 1764*d*n +
720*d) + 720*a*sin(c + d*x)*sin(c + d*x)**n*cos(c + d*x)**4/(d*n**6 + 21*d*n**5 + 175*d*n**4 + 735*d*n**3 + 16
24*d*n**2 + 1764*d*n + 720*d) + b*n**5*sin(c + d*x)**2*sin(c + d*x)**n*cos(c + d*x)**4/(d*n**6 + 21*d*n**5 + 1
75*d*n**4 + 735*d*n**3 + 1624*d*n**2 + 1764*d*n + 720*d) + 4*b*n**4*sin(c + d*x)**4*sin(c + d*x)**n*cos(c + d*
x)**2/(d*n**6 + 21*d*n**5 + 175*d*n**4 + 735*d*n**3 + 1624*d*n**2 + 1764*d*n + 720*d) + 19*b*n**4*sin(c + d*x)
**2*sin(c + d*x)**n*cos(c + d*x)**4/(d*n**6 + 21*d*n**5 + 175*d*n**4 + 735*d*n**3 + 1624*d*n**2 + 1764*d*n + 7
20*d) + 8*b*n**3*sin(c + d*x)**6*sin(c + d*x)**n/(d*n**6 + 21*d*n**5 + 175*d*n**4 + 735*d*n**3 + 1624*d*n**2 +
1764*d*n + 720*d) + 60*b*n**3*sin(c + d*x)**4*sin(c + d*x)**n*cos(c + d*x)**2/(d*n**6 + 21*d*n**5 + 175*d*n**
4 + 735*d*n**3 + 1624*d*n**2 + 1764*d*n + 720*d) + 137*b*n**3*sin(c + d*x)**2*sin(c + d*x)**n*cos(c + d*x)**4/
(d*n**6 + 21*d*n**5 + 175*d*n**4 + 735*d*n**3 + 1624*d*n**2 + 1764*d*n + 720*d) + 72*b*n**2*sin(c + d*x)**6*si
n(c + d*x)**n/(d*n**6 + 21*d*n**5 + 175*d*n**4 + 735*d*n**3 + 1624*d*n**2 + 1764*d*n + 720*d) + 308*b*n**2*sin
(c + d*x)**4*sin(c + d*x)**n*cos(c + d*x)**2/(d*n**6 + 21*d*n**5 + 175*d*n**4 + 735*d*n**3 + 1624*d*n**2 + 176
4*d*n + 720*d) + 461*b*n**2*sin(c + d*x)**2*sin(c + d*x)**n*cos(c + d*x)**4/(d*n**6 + 21*d*n**5 + 175*d*n**4 +
735*d*n**3 + 1624*d*n**2 + 1764*d*n + 720*d) + 184*b*n*sin(c + d*x)**6*sin(c + d*x)**n/(d*n**6 + 21*d*n**5 +
175*d*n**4 + 735*d*n**3 + 1624*d*n**2 + 1764*d*n + 720*d) + 612*b*n*sin(c + d*x)**4*sin(c + d*x)**n*cos(c + d*
x)**2/(d*n**6 + 21*d*n**5 + 175*d*n**4 + 735*d*n**3 + 1624*d*n**2 + 1764*d*n + 720*d) + 702*b*n*sin(c + d*x)**
2*sin(c + d*x)**n*cos(c + d*x)**4/(d*n**6 + 21*d*n**5 + 175*d*n**4 + 735*d*n**3 + 1624*d*n**2 + 1764*d*n + 720
*d) + 120*b*sin(c + d*x)**6*sin(c + d*x)**n/(d*n**6 + 21*d*n**5 + 175*d*n**4 + 735*d*n**3 + 1624*d*n**2 + 1764
*d*n + 720*d) + 360*b*sin(c + d*x)**4*sin(c + d*x)**n*cos(c + d*x)**2/(d*n**6 + 21*d*n**5 + 175*d*n**4 + 735*d
*n**3 + 1624*d*n**2 + 1764*d*n + 720*d) + 360*b*sin(c + d*x)**2*sin(c + d*x)**n*cos(c + d*x)**4/(d*n**6 + 21*d
*n**5 + 175*d*n**4 + 735*d*n**3 + 1624*d*n**2 + 1764*d*n + 720*d), True))

________________________________________________________________________________________