\(\int (a+b \cos (c+d x))^4 (A+B \cos (c+d x)+C \cos ^2(c+d x)) \sec ^8(c+d x) \, dx\) [973]

Optimal result
Mathematica [A] (verified)
Rubi [A] (verified)
Maple [A] (verified)
Fricas [A] (verification not implemented)
Sympy [F(-1)]
Maxima [A] (verification not implemented)
Giac [B] (verification not implemented)
Mupad [B] (verification not implemented)
Reduce [B] (verification not implemented)

Optimal result

Integrand size = 41, antiderivative size = 454 \[ \int (a+b \cos (c+d x))^4 \left (A+B \cos (c+d x)+C \cos ^2(c+d x)\right ) \sec ^8(c+d x) \, dx=\frac {\left (5 a^4 B+36 a^2 b^2 B+8 b^4 B+8 a b^3 (3 A+4 C)+4 a^3 b (5 A+6 C)\right ) \text {arctanh}(\sin (c+d x))}{16 d}+\frac {\left (224 a^3 b B+280 a b^3 B+35 b^4 (2 A+3 C)+84 a^2 b^2 (4 A+5 C)+8 a^4 (6 A+7 C)\right ) \tan (c+d x)}{105 d}+\frac {\left (5 a^4 B+36 a^2 b^2 B+8 b^4 B+8 a b^3 (3 A+4 C)+4 a^3 b (5 A+6 C)\right ) \sec (c+d x) \tan (c+d x)}{16 d}+\frac {\left (4 A b^4+112 a^3 b B+91 a b^3 B+4 a^4 (6 A+7 C)+3 a^2 b^2 (50 A+63 C)\right ) \sec ^2(c+d x) \tan (c+d x)}{105 d}+\frac {a \left (24 A b^3+175 a^3 B+336 a b^2 B+a^2 (412 A b+504 b C)\right ) \sec ^3(c+d x) \tan (c+d x)}{840 d}+\frac {\left (4 A b^2+21 a b B+2 a^2 (6 A+7 C)\right ) (a+b \cos (c+d x))^2 \sec ^4(c+d x) \tan (c+d x)}{70 d}+\frac {(4 A b+7 a B) (a+b \cos (c+d x))^3 \sec ^5(c+d x) \tan (c+d x)}{42 d}+\frac {A (a+b \cos (c+d x))^4 \sec ^6(c+d x) \tan (c+d x)}{7 d} \] Output:

1/16*(5*B*a^4+36*B*a^2*b^2+8*B*b^4+8*a*b^3*(3*A+4*C)+4*a^3*b*(5*A+6*C))*ar 
ctanh(sin(d*x+c))/d+1/105*(224*B*a^3*b+280*B*a*b^3+35*b^4*(2*A+3*C)+84*a^2 
*b^2*(4*A+5*C)+8*a^4*(6*A+7*C))*tan(d*x+c)/d+1/16*(5*B*a^4+36*B*a^2*b^2+8* 
B*b^4+8*a*b^3*(3*A+4*C)+4*a^3*b*(5*A+6*C))*sec(d*x+c)*tan(d*x+c)/d+1/105*( 
4*A*b^4+112*B*a^3*b+91*B*a*b^3+4*a^4*(6*A+7*C)+3*a^2*b^2*(50*A+63*C))*sec( 
d*x+c)^2*tan(d*x+c)/d+1/840*a*(24*A*b^3+175*B*a^3+336*B*a*b^2+a^2*(412*A*b 
+504*C*b))*sec(d*x+c)^3*tan(d*x+c)/d+1/70*(4*A*b^2+21*B*a*b+2*a^2*(6*A+7*C 
))*(a+b*cos(d*x+c))^2*sec(d*x+c)^4*tan(d*x+c)/d+1/42*(4*A*b+7*B*a)*(a+b*co 
s(d*x+c))^3*sec(d*x+c)^5*tan(d*x+c)/d+1/7*A*(a+b*cos(d*x+c))^4*sec(d*x+c)^ 
6*tan(d*x+c)/d
                                                                                    
                                                                                    
 

Mathematica [A] (verified)

Time = 3.68 (sec) , antiderivative size = 341, normalized size of antiderivative = 0.75 \[ \int (a+b \cos (c+d x))^4 \left (A+B \cos (c+d x)+C \cos ^2(c+d x)\right ) \sec ^8(c+d x) \, dx=\frac {105 \left (5 a^4 B+36 a^2 b^2 B+8 b^4 B+8 a b^3 (3 A+4 C)+4 a^3 b (5 A+6 C)\right ) \text {arctanh}(\sin (c+d x))+\tan (c+d x) \left (105 \left (5 a^4 B+36 a^2 b^2 B+8 b^4 B+8 a b^3 (3 A+4 C)+4 a^3 b (5 A+6 C)\right ) \sec (c+d x)+70 a \left (24 A b^3+5 a^3 B+36 a b^2 B+4 a^2 b (5 A+6 C)\right ) \sec ^3(c+d x)+280 a^3 (4 A b+a B) \sec ^5(c+d x)+16 \left (105 \left (4 a^3 b B+4 a b^3 B+a^4 (A+C)+6 a^2 b^2 (A+C)+b^4 (A+C)\right )+35 \left (A b^4+8 a^3 b B+4 a b^3 B+6 a^2 b^2 (2 A+C)+a^4 (3 A+2 C)\right ) \tan ^2(c+d x)+21 a^2 \left (6 A b^2+4 a b B+a^2 (3 A+C)\right ) \tan ^4(c+d x)+15 a^4 A \tan ^6(c+d x)\right )\right )}{1680 d} \] Input:

Integrate[(a + b*Cos[c + d*x])^4*(A + B*Cos[c + d*x] + C*Cos[c + d*x]^2)*S 
ec[c + d*x]^8,x]
 

Output:

(105*(5*a^4*B + 36*a^2*b^2*B + 8*b^4*B + 8*a*b^3*(3*A + 4*C) + 4*a^3*b*(5* 
A + 6*C))*ArcTanh[Sin[c + d*x]] + Tan[c + d*x]*(105*(5*a^4*B + 36*a^2*b^2* 
B + 8*b^4*B + 8*a*b^3*(3*A + 4*C) + 4*a^3*b*(5*A + 6*C))*Sec[c + d*x] + 70 
*a*(24*A*b^3 + 5*a^3*B + 36*a*b^2*B + 4*a^2*b*(5*A + 6*C))*Sec[c + d*x]^3 
+ 280*a^3*(4*A*b + a*B)*Sec[c + d*x]^5 + 16*(105*(4*a^3*b*B + 4*a*b^3*B + 
a^4*(A + C) + 6*a^2*b^2*(A + C) + b^4*(A + C)) + 35*(A*b^4 + 8*a^3*b*B + 4 
*a*b^3*B + 6*a^2*b^2*(2*A + C) + a^4*(3*A + 2*C))*Tan[c + d*x]^2 + 21*a^2* 
(6*A*b^2 + 4*a*b*B + a^2*(3*A + C))*Tan[c + d*x]^4 + 15*a^4*A*Tan[c + d*x] 
^6)))/(1680*d)
 

Rubi [A] (verified)

Time = 2.98 (sec) , antiderivative size = 425, normalized size of antiderivative = 0.94, number of steps used = 21, number of rules used = 20, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.488, Rules used = {3042, 3526, 3042, 3526, 3042, 3526, 3042, 3510, 25, 3042, 3500, 27, 3042, 3227, 3042, 4254, 24, 4255, 3042, 4257}

Below are the steps used by Rubi to obtain the solution. The rule number used for the transformation is given above next to the arrow. The rules definitions used are listed below.

\(\displaystyle \int \sec ^8(c+d x) (a+b \cos (c+d x))^4 \left (A+B \cos (c+d x)+C \cos ^2(c+d x)\right ) \, dx\)

\(\Big \downarrow \) 3042

\(\displaystyle \int \frac {\left (a+b \sin \left (c+d x+\frac {\pi }{2}\right )\right )^4 \left (A+B \sin \left (c+d x+\frac {\pi }{2}\right )+C \sin \left (c+d x+\frac {\pi }{2}\right )^2\right )}{\sin \left (c+d x+\frac {\pi }{2}\right )^8}dx\)

\(\Big \downarrow \) 3526

\(\displaystyle \frac {1}{7} \int (a+b \cos (c+d x))^3 \left (b (2 A+7 C) \cos ^2(c+d x)+(6 a A+7 b B+7 a C) \cos (c+d x)+4 A b+7 a B\right ) \sec ^7(c+d x)dx+\frac {A \tan (c+d x) \sec ^6(c+d x) (a+b \cos (c+d x))^4}{7 d}\)

\(\Big \downarrow \) 3042

\(\displaystyle \frac {1}{7} \int \frac {\left (a+b \sin \left (c+d x+\frac {\pi }{2}\right )\right )^3 \left (b (2 A+7 C) \sin \left (c+d x+\frac {\pi }{2}\right )^2+(6 a A+7 b B+7 a C) \sin \left (c+d x+\frac {\pi }{2}\right )+4 A b+7 a B\right )}{\sin \left (c+d x+\frac {\pi }{2}\right )^7}dx+\frac {A \tan (c+d x) \sec ^6(c+d x) (a+b \cos (c+d x))^4}{7 d}\)

\(\Big \downarrow \) 3526

\(\displaystyle \frac {1}{7} \left (\frac {1}{6} \int (a+b \cos (c+d x))^2 \left (2 b (10 A b+21 C b+7 a B) \cos ^2(c+d x)+\left (35 B a^2+68 A b a+84 b C a+42 b^2 B\right ) \cos (c+d x)+3 \left (2 (6 A+7 C) a^2+21 b B a+4 A b^2\right )\right ) \sec ^6(c+d x)dx+\frac {(7 a B+4 A b) \tan (c+d x) \sec ^5(c+d x) (a+b \cos (c+d x))^3}{6 d}\right )+\frac {A \tan (c+d x) \sec ^6(c+d x) (a+b \cos (c+d x))^4}{7 d}\)

\(\Big \downarrow \) 3042

\(\displaystyle \frac {1}{7} \left (\frac {1}{6} \int \frac {\left (a+b \sin \left (c+d x+\frac {\pi }{2}\right )\right )^2 \left (2 b (10 A b+21 C b+7 a B) \sin \left (c+d x+\frac {\pi }{2}\right )^2+\left (35 B a^2+68 A b a+84 b C a+42 b^2 B\right ) \sin \left (c+d x+\frac {\pi }{2}\right )+3 \left (2 (6 A+7 C) a^2+21 b B a+4 A b^2\right )\right )}{\sin \left (c+d x+\frac {\pi }{2}\right )^6}dx+\frac {(7 a B+4 A b) \tan (c+d x) \sec ^5(c+d x) (a+b \cos (c+d x))^3}{6 d}\right )+\frac {A \tan (c+d x) \sec ^6(c+d x) (a+b \cos (c+d x))^4}{7 d}\)

\(\Big \downarrow \) 3526

\(\displaystyle \frac {1}{7} \left (\frac {1}{6} \left (\frac {1}{5} \int (a+b \cos (c+d x)) \left (175 B a^3+(412 A b+504 C b) a^2+336 b^2 B a+24 A b^3+2 b \left (6 (6 A+7 C) a^2+98 b B a+b^2 (62 A+105 C)\right ) \cos ^2(c+d x)+\left (24 (6 A+7 C) a^3+497 b B a^2+2 b^2 (244 A+315 C) a+210 b^3 B\right ) \cos (c+d x)\right ) \sec ^5(c+d x)dx+\frac {3 \tan (c+d x) \sec ^4(c+d x) \left (2 a^2 (6 A+7 C)+21 a b B+4 A b^2\right ) (a+b \cos (c+d x))^2}{5 d}\right )+\frac {(7 a B+4 A b) \tan (c+d x) \sec ^5(c+d x) (a+b \cos (c+d x))^3}{6 d}\right )+\frac {A \tan (c+d x) \sec ^6(c+d x) (a+b \cos (c+d x))^4}{7 d}\)

\(\Big \downarrow \) 3042

\(\displaystyle \frac {1}{7} \left (\frac {1}{6} \left (\frac {1}{5} \int \frac {\left (a+b \sin \left (c+d x+\frac {\pi }{2}\right )\right ) \left (175 B a^3+(412 A b+504 C b) a^2+336 b^2 B a+24 A b^3+2 b \left (6 (6 A+7 C) a^2+98 b B a+b^2 (62 A+105 C)\right ) \sin \left (c+d x+\frac {\pi }{2}\right )^2+\left (24 (6 A+7 C) a^3+497 b B a^2+2 b^2 (244 A+315 C) a+210 b^3 B\right ) \sin \left (c+d x+\frac {\pi }{2}\right )\right )}{\sin \left (c+d x+\frac {\pi }{2}\right )^5}dx+\frac {3 \tan (c+d x) \sec ^4(c+d x) \left (2 a^2 (6 A+7 C)+21 a b B+4 A b^2\right ) (a+b \cos (c+d x))^2}{5 d}\right )+\frac {(7 a B+4 A b) \tan (c+d x) \sec ^5(c+d x) (a+b \cos (c+d x))^3}{6 d}\right )+\frac {A \tan (c+d x) \sec ^6(c+d x) (a+b \cos (c+d x))^4}{7 d}\)

\(\Big \downarrow \) 3510

\(\displaystyle \frac {1}{7} \left (\frac {1}{6} \left (\frac {1}{5} \left (\frac {a \tan (c+d x) \sec ^3(c+d x) \left (175 a^3 B+a^2 (412 A b+504 b C)+336 a b^2 B+24 A b^3\right )}{4 d}-\frac {1}{4} \int -\left (\left (8 b^2 \left (6 (6 A+7 C) a^2+98 b B a+b^2 (62 A+105 C)\right ) \cos ^2(c+d x)+105 \left (5 B a^4+4 b (5 A+6 C) a^3+36 b^2 B a^2+8 b^3 (3 A+4 C) a+8 b^4 B\right ) \cos (c+d x)+24 \left (4 (6 A+7 C) a^4+112 b B a^3+3 b^2 (50 A+63 C) a^2+91 b^3 B a+4 A b^4\right )\right ) \sec ^4(c+d x)\right )dx\right )+\frac {3 \tan (c+d x) \sec ^4(c+d x) \left (2 a^2 (6 A+7 C)+21 a b B+4 A b^2\right ) (a+b \cos (c+d x))^2}{5 d}\right )+\frac {(7 a B+4 A b) \tan (c+d x) \sec ^5(c+d x) (a+b \cos (c+d x))^3}{6 d}\right )+\frac {A \tan (c+d x) \sec ^6(c+d x) (a+b \cos (c+d x))^4}{7 d}\)

\(\Big \downarrow \) 25

\(\displaystyle \frac {1}{7} \left (\frac {1}{6} \left (\frac {1}{5} \left (\frac {1}{4} \int \left (8 b^2 \left (6 (6 A+7 C) a^2+98 b B a+b^2 (62 A+105 C)\right ) \cos ^2(c+d x)+105 \left (5 B a^4+4 b (5 A+6 C) a^3+36 b^2 B a^2+8 b^3 (3 A+4 C) a+8 b^4 B\right ) \cos (c+d x)+24 \left (4 (6 A+7 C) a^4+112 b B a^3+3 b^2 (50 A+63 C) a^2+91 b^3 B a+4 A b^4\right )\right ) \sec ^4(c+d x)dx+\frac {a \tan (c+d x) \sec ^3(c+d x) \left (175 a^3 B+a^2 (412 A b+504 b C)+336 a b^2 B+24 A b^3\right )}{4 d}\right )+\frac {3 \tan (c+d x) \sec ^4(c+d x) \left (2 a^2 (6 A+7 C)+21 a b B+4 A b^2\right ) (a+b \cos (c+d x))^2}{5 d}\right )+\frac {(7 a B+4 A b) \tan (c+d x) \sec ^5(c+d x) (a+b \cos (c+d x))^3}{6 d}\right )+\frac {A \tan (c+d x) \sec ^6(c+d x) (a+b \cos (c+d x))^4}{7 d}\)

\(\Big \downarrow \) 3042

\(\displaystyle \frac {1}{7} \left (\frac {1}{6} \left (\frac {1}{5} \left (\frac {1}{4} \int \frac {8 b^2 \left (6 (6 A+7 C) a^2+98 b B a+b^2 (62 A+105 C)\right ) \sin \left (c+d x+\frac {\pi }{2}\right )^2+105 \left (5 B a^4+4 b (5 A+6 C) a^3+36 b^2 B a^2+8 b^3 (3 A+4 C) a+8 b^4 B\right ) \sin \left (c+d x+\frac {\pi }{2}\right )+24 \left (4 (6 A+7 C) a^4+112 b B a^3+3 b^2 (50 A+63 C) a^2+91 b^3 B a+4 A b^4\right )}{\sin \left (c+d x+\frac {\pi }{2}\right )^4}dx+\frac {a \tan (c+d x) \sec ^3(c+d x) \left (175 a^3 B+a^2 (412 A b+504 b C)+336 a b^2 B+24 A b^3\right )}{4 d}\right )+\frac {3 \tan (c+d x) \sec ^4(c+d x) \left (2 a^2 (6 A+7 C)+21 a b B+4 A b^2\right ) (a+b \cos (c+d x))^2}{5 d}\right )+\frac {(7 a B+4 A b) \tan (c+d x) \sec ^5(c+d x) (a+b \cos (c+d x))^3}{6 d}\right )+\frac {A \tan (c+d x) \sec ^6(c+d x) (a+b \cos (c+d x))^4}{7 d}\)

\(\Big \downarrow \) 3500

\(\displaystyle \frac {1}{7} \left (\frac {1}{6} \left (\frac {1}{5} \left (\frac {1}{4} \left (\frac {1}{3} \int 3 \left (105 \left (5 B a^4+4 b (5 A+6 C) a^3+36 b^2 B a^2+8 b^3 (3 A+4 C) a+8 b^4 B\right )+8 \left (8 (6 A+7 C) a^4+224 b B a^3+84 b^2 (4 A+5 C) a^2+280 b^3 B a+35 b^4 (2 A+3 C)\right ) \cos (c+d x)\right ) \sec ^3(c+d x)dx+\frac {8 \tan (c+d x) \sec ^2(c+d x) \left (4 a^4 (6 A+7 C)+112 a^3 b B+3 a^2 b^2 (50 A+63 C)+91 a b^3 B+4 A b^4\right )}{d}\right )+\frac {a \tan (c+d x) \sec ^3(c+d x) \left (175 a^3 B+a^2 (412 A b+504 b C)+336 a b^2 B+24 A b^3\right )}{4 d}\right )+\frac {3 \tan (c+d x) \sec ^4(c+d x) \left (2 a^2 (6 A+7 C)+21 a b B+4 A b^2\right ) (a+b \cos (c+d x))^2}{5 d}\right )+\frac {(7 a B+4 A b) \tan (c+d x) \sec ^5(c+d x) (a+b \cos (c+d x))^3}{6 d}\right )+\frac {A \tan (c+d x) \sec ^6(c+d x) (a+b \cos (c+d x))^4}{7 d}\)

\(\Big \downarrow \) 27

\(\displaystyle \frac {1}{7} \left (\frac {1}{6} \left (\frac {1}{5} \left (\frac {1}{4} \left (\int \left (105 \left (5 B a^4+4 b (5 A+6 C) a^3+36 b^2 B a^2+8 b^3 (3 A+4 C) a+8 b^4 B\right )+8 \left (8 (6 A+7 C) a^4+224 b B a^3+84 b^2 (4 A+5 C) a^2+280 b^3 B a+35 b^4 (2 A+3 C)\right ) \cos (c+d x)\right ) \sec ^3(c+d x)dx+\frac {8 \tan (c+d x) \sec ^2(c+d x) \left (4 a^4 (6 A+7 C)+112 a^3 b B+3 a^2 b^2 (50 A+63 C)+91 a b^3 B+4 A b^4\right )}{d}\right )+\frac {a \tan (c+d x) \sec ^3(c+d x) \left (175 a^3 B+a^2 (412 A b+504 b C)+336 a b^2 B+24 A b^3\right )}{4 d}\right )+\frac {3 \tan (c+d x) \sec ^4(c+d x) \left (2 a^2 (6 A+7 C)+21 a b B+4 A b^2\right ) (a+b \cos (c+d x))^2}{5 d}\right )+\frac {(7 a B+4 A b) \tan (c+d x) \sec ^5(c+d x) (a+b \cos (c+d x))^3}{6 d}\right )+\frac {A \tan (c+d x) \sec ^6(c+d x) (a+b \cos (c+d x))^4}{7 d}\)

\(\Big \downarrow \) 3042

\(\displaystyle \frac {1}{7} \left (\frac {1}{6} \left (\frac {1}{5} \left (\frac {1}{4} \left (\int \frac {105 \left (5 B a^4+4 b (5 A+6 C) a^3+36 b^2 B a^2+8 b^3 (3 A+4 C) a+8 b^4 B\right )+8 \left (8 (6 A+7 C) a^4+224 b B a^3+84 b^2 (4 A+5 C) a^2+280 b^3 B a+35 b^4 (2 A+3 C)\right ) \sin \left (c+d x+\frac {\pi }{2}\right )}{\sin \left (c+d x+\frac {\pi }{2}\right )^3}dx+\frac {8 \tan (c+d x) \sec ^2(c+d x) \left (4 a^4 (6 A+7 C)+112 a^3 b B+3 a^2 b^2 (50 A+63 C)+91 a b^3 B+4 A b^4\right )}{d}\right )+\frac {a \tan (c+d x) \sec ^3(c+d x) \left (175 a^3 B+a^2 (412 A b+504 b C)+336 a b^2 B+24 A b^3\right )}{4 d}\right )+\frac {3 \tan (c+d x) \sec ^4(c+d x) \left (2 a^2 (6 A+7 C)+21 a b B+4 A b^2\right ) (a+b \cos (c+d x))^2}{5 d}\right )+\frac {(7 a B+4 A b) \tan (c+d x) \sec ^5(c+d x) (a+b \cos (c+d x))^3}{6 d}\right )+\frac {A \tan (c+d x) \sec ^6(c+d x) (a+b \cos (c+d x))^4}{7 d}\)

\(\Big \downarrow \) 3227

\(\displaystyle \frac {1}{7} \left (\frac {1}{6} \left (\frac {1}{5} \left (\frac {1}{4} \left (105 \left (5 a^4 B+4 a^3 b (5 A+6 C)+36 a^2 b^2 B+8 a b^3 (3 A+4 C)+8 b^4 B\right ) \int \sec ^3(c+d x)dx+8 \left (8 a^4 (6 A+7 C)+224 a^3 b B+84 a^2 b^2 (4 A+5 C)+280 a b^3 B+35 b^4 (2 A+3 C)\right ) \int \sec ^2(c+d x)dx+\frac {8 \tan (c+d x) \sec ^2(c+d x) \left (4 a^4 (6 A+7 C)+112 a^3 b B+3 a^2 b^2 (50 A+63 C)+91 a b^3 B+4 A b^4\right )}{d}\right )+\frac {a \tan (c+d x) \sec ^3(c+d x) \left (175 a^3 B+a^2 (412 A b+504 b C)+336 a b^2 B+24 A b^3\right )}{4 d}\right )+\frac {3 \tan (c+d x) \sec ^4(c+d x) \left (2 a^2 (6 A+7 C)+21 a b B+4 A b^2\right ) (a+b \cos (c+d x))^2}{5 d}\right )+\frac {(7 a B+4 A b) \tan (c+d x) \sec ^5(c+d x) (a+b \cos (c+d x))^3}{6 d}\right )+\frac {A \tan (c+d x) \sec ^6(c+d x) (a+b \cos (c+d x))^4}{7 d}\)

\(\Big \downarrow \) 3042

\(\displaystyle \frac {1}{7} \left (\frac {1}{6} \left (\frac {1}{5} \left (\frac {1}{4} \left (8 \left (8 a^4 (6 A+7 C)+224 a^3 b B+84 a^2 b^2 (4 A+5 C)+280 a b^3 B+35 b^4 (2 A+3 C)\right ) \int \csc \left (c+d x+\frac {\pi }{2}\right )^2dx+105 \left (5 a^4 B+4 a^3 b (5 A+6 C)+36 a^2 b^2 B+8 a b^3 (3 A+4 C)+8 b^4 B\right ) \int \csc \left (c+d x+\frac {\pi }{2}\right )^3dx+\frac {8 \tan (c+d x) \sec ^2(c+d x) \left (4 a^4 (6 A+7 C)+112 a^3 b B+3 a^2 b^2 (50 A+63 C)+91 a b^3 B+4 A b^4\right )}{d}\right )+\frac {a \tan (c+d x) \sec ^3(c+d x) \left (175 a^3 B+a^2 (412 A b+504 b C)+336 a b^2 B+24 A b^3\right )}{4 d}\right )+\frac {3 \tan (c+d x) \sec ^4(c+d x) \left (2 a^2 (6 A+7 C)+21 a b B+4 A b^2\right ) (a+b \cos (c+d x))^2}{5 d}\right )+\frac {(7 a B+4 A b) \tan (c+d x) \sec ^5(c+d x) (a+b \cos (c+d x))^3}{6 d}\right )+\frac {A \tan (c+d x) \sec ^6(c+d x) (a+b \cos (c+d x))^4}{7 d}\)

\(\Big \downarrow \) 4254

\(\displaystyle \frac {1}{7} \left (\frac {1}{6} \left (\frac {1}{5} \left (\frac {1}{4} \left (-\frac {8 \left (8 a^4 (6 A+7 C)+224 a^3 b B+84 a^2 b^2 (4 A+5 C)+280 a b^3 B+35 b^4 (2 A+3 C)\right ) \int 1d(-\tan (c+d x))}{d}+105 \left (5 a^4 B+4 a^3 b (5 A+6 C)+36 a^2 b^2 B+8 a b^3 (3 A+4 C)+8 b^4 B\right ) \int \csc \left (c+d x+\frac {\pi }{2}\right )^3dx+\frac {8 \tan (c+d x) \sec ^2(c+d x) \left (4 a^4 (6 A+7 C)+112 a^3 b B+3 a^2 b^2 (50 A+63 C)+91 a b^3 B+4 A b^4\right )}{d}\right )+\frac {a \tan (c+d x) \sec ^3(c+d x) \left (175 a^3 B+a^2 (412 A b+504 b C)+336 a b^2 B+24 A b^3\right )}{4 d}\right )+\frac {3 \tan (c+d x) \sec ^4(c+d x) \left (2 a^2 (6 A+7 C)+21 a b B+4 A b^2\right ) (a+b \cos (c+d x))^2}{5 d}\right )+\frac {(7 a B+4 A b) \tan (c+d x) \sec ^5(c+d x) (a+b \cos (c+d x))^3}{6 d}\right )+\frac {A \tan (c+d x) \sec ^6(c+d x) (a+b \cos (c+d x))^4}{7 d}\)

\(\Big \downarrow \) 24

\(\displaystyle \frac {1}{7} \left (\frac {1}{6} \left (\frac {1}{5} \left (\frac {1}{4} \left (105 \left (5 a^4 B+4 a^3 b (5 A+6 C)+36 a^2 b^2 B+8 a b^3 (3 A+4 C)+8 b^4 B\right ) \int \csc \left (c+d x+\frac {\pi }{2}\right )^3dx+\frac {8 \tan (c+d x) \left (8 a^4 (6 A+7 C)+224 a^3 b B+84 a^2 b^2 (4 A+5 C)+280 a b^3 B+35 b^4 (2 A+3 C)\right )}{d}+\frac {8 \tan (c+d x) \sec ^2(c+d x) \left (4 a^4 (6 A+7 C)+112 a^3 b B+3 a^2 b^2 (50 A+63 C)+91 a b^3 B+4 A b^4\right )}{d}\right )+\frac {a \tan (c+d x) \sec ^3(c+d x) \left (175 a^3 B+a^2 (412 A b+504 b C)+336 a b^2 B+24 A b^3\right )}{4 d}\right )+\frac {3 \tan (c+d x) \sec ^4(c+d x) \left (2 a^2 (6 A+7 C)+21 a b B+4 A b^2\right ) (a+b \cos (c+d x))^2}{5 d}\right )+\frac {(7 a B+4 A b) \tan (c+d x) \sec ^5(c+d x) (a+b \cos (c+d x))^3}{6 d}\right )+\frac {A \tan (c+d x) \sec ^6(c+d x) (a+b \cos (c+d x))^4}{7 d}\)

\(\Big \downarrow \) 4255

\(\displaystyle \frac {1}{7} \left (\frac {1}{6} \left (\frac {1}{5} \left (\frac {1}{4} \left (105 \left (5 a^4 B+4 a^3 b (5 A+6 C)+36 a^2 b^2 B+8 a b^3 (3 A+4 C)+8 b^4 B\right ) \left (\frac {1}{2} \int \sec (c+d x)dx+\frac {\tan (c+d x) \sec (c+d x)}{2 d}\right )+\frac {8 \tan (c+d x) \left (8 a^4 (6 A+7 C)+224 a^3 b B+84 a^2 b^2 (4 A+5 C)+280 a b^3 B+35 b^4 (2 A+3 C)\right )}{d}+\frac {8 \tan (c+d x) \sec ^2(c+d x) \left (4 a^4 (6 A+7 C)+112 a^3 b B+3 a^2 b^2 (50 A+63 C)+91 a b^3 B+4 A b^4\right )}{d}\right )+\frac {a \tan (c+d x) \sec ^3(c+d x) \left (175 a^3 B+a^2 (412 A b+504 b C)+336 a b^2 B+24 A b^3\right )}{4 d}\right )+\frac {3 \tan (c+d x) \sec ^4(c+d x) \left (2 a^2 (6 A+7 C)+21 a b B+4 A b^2\right ) (a+b \cos (c+d x))^2}{5 d}\right )+\frac {(7 a B+4 A b) \tan (c+d x) \sec ^5(c+d x) (a+b \cos (c+d x))^3}{6 d}\right )+\frac {A \tan (c+d x) \sec ^6(c+d x) (a+b \cos (c+d x))^4}{7 d}\)

\(\Big \downarrow \) 3042

\(\displaystyle \frac {1}{7} \left (\frac {1}{6} \left (\frac {1}{5} \left (\frac {1}{4} \left (105 \left (5 a^4 B+4 a^3 b (5 A+6 C)+36 a^2 b^2 B+8 a b^3 (3 A+4 C)+8 b^4 B\right ) \left (\frac {1}{2} \int \csc \left (c+d x+\frac {\pi }{2}\right )dx+\frac {\tan (c+d x) \sec (c+d x)}{2 d}\right )+\frac {8 \tan (c+d x) \left (8 a^4 (6 A+7 C)+224 a^3 b B+84 a^2 b^2 (4 A+5 C)+280 a b^3 B+35 b^4 (2 A+3 C)\right )}{d}+\frac {8 \tan (c+d x) \sec ^2(c+d x) \left (4 a^4 (6 A+7 C)+112 a^3 b B+3 a^2 b^2 (50 A+63 C)+91 a b^3 B+4 A b^4\right )}{d}\right )+\frac {a \tan (c+d x) \sec ^3(c+d x) \left (175 a^3 B+a^2 (412 A b+504 b C)+336 a b^2 B+24 A b^3\right )}{4 d}\right )+\frac {3 \tan (c+d x) \sec ^4(c+d x) \left (2 a^2 (6 A+7 C)+21 a b B+4 A b^2\right ) (a+b \cos (c+d x))^2}{5 d}\right )+\frac {(7 a B+4 A b) \tan (c+d x) \sec ^5(c+d x) (a+b \cos (c+d x))^3}{6 d}\right )+\frac {A \tan (c+d x) \sec ^6(c+d x) (a+b \cos (c+d x))^4}{7 d}\)

\(\Big \downarrow \) 4257

\(\displaystyle \frac {1}{7} \left (\frac {1}{6} \left (\frac {3 \tan (c+d x) \sec ^4(c+d x) \left (2 a^2 (6 A+7 C)+21 a b B+4 A b^2\right ) (a+b \cos (c+d x))^2}{5 d}+\frac {1}{5} \left (\frac {a \tan (c+d x) \sec ^3(c+d x) \left (175 a^3 B+a^2 (412 A b+504 b C)+336 a b^2 B+24 A b^3\right )}{4 d}+\frac {1}{4} \left (105 \left (5 a^4 B+4 a^3 b (5 A+6 C)+36 a^2 b^2 B+8 a b^3 (3 A+4 C)+8 b^4 B\right ) \left (\frac {\text {arctanh}(\sin (c+d x))}{2 d}+\frac {\tan (c+d x) \sec (c+d x)}{2 d}\right )+\frac {8 \tan (c+d x) \left (8 a^4 (6 A+7 C)+224 a^3 b B+84 a^2 b^2 (4 A+5 C)+280 a b^3 B+35 b^4 (2 A+3 C)\right )}{d}+\frac {8 \tan (c+d x) \sec ^2(c+d x) \left (4 a^4 (6 A+7 C)+112 a^3 b B+3 a^2 b^2 (50 A+63 C)+91 a b^3 B+4 A b^4\right )}{d}\right )\right )\right )+\frac {(7 a B+4 A b) \tan (c+d x) \sec ^5(c+d x) (a+b \cos (c+d x))^3}{6 d}\right )+\frac {A \tan (c+d x) \sec ^6(c+d x) (a+b \cos (c+d x))^4}{7 d}\)

Input:

Int[(a + b*Cos[c + d*x])^4*(A + B*Cos[c + d*x] + C*Cos[c + d*x]^2)*Sec[c + 
 d*x]^8,x]
 

Output:

(A*(a + b*Cos[c + d*x])^4*Sec[c + d*x]^6*Tan[c + d*x])/(7*d) + (((4*A*b + 
7*a*B)*(a + b*Cos[c + d*x])^3*Sec[c + d*x]^5*Tan[c + d*x])/(6*d) + ((3*(4* 
A*b^2 + 21*a*b*B + 2*a^2*(6*A + 7*C))*(a + b*Cos[c + d*x])^2*Sec[c + d*x]^ 
4*Tan[c + d*x])/(5*d) + ((a*(24*A*b^3 + 175*a^3*B + 336*a*b^2*B + a^2*(412 
*A*b + 504*b*C))*Sec[c + d*x]^3*Tan[c + d*x])/(4*d) + ((8*(224*a^3*b*B + 2 
80*a*b^3*B + 35*b^4*(2*A + 3*C) + 84*a^2*b^2*(4*A + 5*C) + 8*a^4*(6*A + 7* 
C))*Tan[c + d*x])/d + (8*(4*A*b^4 + 112*a^3*b*B + 91*a*b^3*B + 4*a^4*(6*A 
+ 7*C) + 3*a^2*b^2*(50*A + 63*C))*Sec[c + d*x]^2*Tan[c + d*x])/d + 105*(5* 
a^4*B + 36*a^2*b^2*B + 8*b^4*B + 8*a*b^3*(3*A + 4*C) + 4*a^3*b*(5*A + 6*C) 
)*(ArcTanh[Sin[c + d*x]]/(2*d) + (Sec[c + d*x]*Tan[c + d*x])/(2*d)))/4)/5) 
/6)/7
 

Defintions of rubi rules used

rule 24
Int[a_, x_Symbol] :> Simp[a*x, x] /; FreeQ[a, x]
 

rule 25
Int[-(Fx_), x_Symbol] :> Simp[Identity[-1]   Int[Fx, x], x]
 

rule 27
Int[(a_)*(Fx_), x_Symbol] :> Simp[a   Int[Fx, x], x] /; FreeQ[a, x] &&  !Ma 
tchQ[Fx, (b_)*(Gx_) /; FreeQ[b, x]]
 

rule 3042
Int[u_, x_Symbol] :> Int[DeactivateTrig[u, x], x] /; FunctionOfTrigOfLinear 
Q[u, x]
 

rule 3227
Int[((b_.)*sin[(e_.) + (f_.)*(x_)])^(m_)*((c_) + (d_.)*sin[(e_.) + (f_.)*(x 
_)]), x_Symbol] :> Simp[c   Int[(b*Sin[e + f*x])^m, x], x] + Simp[d/b   Int 
[(b*Sin[e + f*x])^(m + 1), x], x] /; FreeQ[{b, c, d, e, f, m}, x]
 

rule 3500
Int[((a_.) + (b_.)*sin[(e_.) + (f_.)*(x_)])^(m_)*((A_.) + (B_.)*sin[(e_.) + 
 (f_.)*(x_)] + (C_.)*sin[(e_.) + (f_.)*(x_)]^2), x_Symbol] :> Simp[(-(A*b^2 
 - a*b*B + a^2*C))*Cos[e + f*x]*((a + b*Sin[e + f*x])^(m + 1)/(b*f*(m + 1)* 
(a^2 - b^2))), x] + Simp[1/(b*(m + 1)*(a^2 - b^2))   Int[(a + b*Sin[e + f*x 
])^(m + 1)*Simp[b*(a*A - b*B + a*C)*(m + 1) - (A*b^2 - a*b*B + a^2*C + b*(A 
*b - a*B + b*C)*(m + 1))*Sin[e + f*x], x], x], x] /; FreeQ[{a, b, e, f, A, 
B, C}, x] && LtQ[m, -1] && NeQ[a^2 - b^2, 0]
 

rule 3510
Int[((a_.) + (b_.)*sin[(e_.) + (f_.)*(x_)])^(m_)*((c_.) + (d_.)*sin[(e_.) + 
 (f_.)*(x_)])*((A_.) + (B_.)*sin[(e_.) + (f_.)*(x_)] + (C_.)*sin[(e_.) + (f 
_.)*(x_)]^2), x_Symbol] :> Simp[(-(b*c - a*d))*(A*b^2 - a*b*B + a^2*C)*Cos[ 
e + f*x]*((a + b*Sin[e + f*x])^(m + 1)/(b^2*f*(m + 1)*(a^2 - b^2))), x] - S 
imp[1/(b^2*(m + 1)*(a^2 - b^2))   Int[(a + b*Sin[e + f*x])^(m + 1)*Simp[b*( 
m + 1)*((b*B - a*C)*(b*c - a*d) - A*b*(a*c - b*d)) + (b*B*(a^2*d + b^2*d*(m 
 + 1) - a*b*c*(m + 2)) + (b*c - a*d)*(A*b^2*(m + 2) + C*(a^2 + b^2*(m + 1)) 
))*Sin[e + f*x] - b*C*d*(m + 1)*(a^2 - b^2)*Sin[e + f*x]^2, x], x], x] /; F 
reeQ[{a, b, c, d, e, f, A, B, C}, x] && NeQ[b*c - a*d, 0] && NeQ[a^2 - b^2, 
 0] && LtQ[m, -1]
 

rule 3526
Int[((a_.) + (b_.)*sin[(e_.) + (f_.)*(x_)])^(m_)*((c_.) + (d_.)*sin[(e_.) + 
 (f_.)*(x_)])^(n_)*((A_.) + (B_.)*sin[(e_.) + (f_.)*(x_)] + (C_.)*sin[(e_.) 
 + (f_.)*(x_)]^2), x_Symbol] :> Simp[(-(c^2*C - B*c*d + A*d^2))*Cos[e + f*x 
]*(a + b*Sin[e + f*x])^m*((c + d*Sin[e + f*x])^(n + 1)/(d*f*(n + 1)*(c^2 - 
d^2))), x] + Simp[1/(d*(n + 1)*(c^2 - d^2))   Int[(a + b*Sin[e + f*x])^(m - 
 1)*(c + d*Sin[e + f*x])^(n + 1)*Simp[A*d*(b*d*m + a*c*(n + 1)) + (c*C - B* 
d)*(b*c*m + a*d*(n + 1)) - (d*(A*(a*d*(n + 2) - b*c*(n + 1)) + B*(b*d*(n + 
1) - a*c*(n + 2))) - C*(b*c*d*(n + 1) - a*(c^2 + d^2*(n + 1))))*Sin[e + f*x 
] + b*(d*(B*c - A*d)*(m + n + 2) - C*(c^2*(m + 1) + d^2*(n + 1)))*Sin[e + f 
*x]^2, x], x], x] /; FreeQ[{a, b, c, d, e, f, A, B, C}, x] && NeQ[b*c - a*d 
, 0] && NeQ[a^2 - b^2, 0] && NeQ[c^2 - d^2, 0] && GtQ[m, 0] && LtQ[n, -1]
 

rule 4254
Int[csc[(c_.) + (d_.)*(x_)]^(n_), x_Symbol] :> Simp[-d^(-1)   Subst[Int[Exp 
andIntegrand[(1 + x^2)^(n/2 - 1), x], x], x, Cot[c + d*x]], x] /; FreeQ[{c, 
 d}, x] && IGtQ[n/2, 0]
 

rule 4255
Int[(csc[(c_.) + (d_.)*(x_)]*(b_.))^(n_), x_Symbol] :> Simp[(-b)*Cos[c + d* 
x]*((b*Csc[c + d*x])^(n - 1)/(d*(n - 1))), x] + Simp[b^2*((n - 2)/(n - 1)) 
  Int[(b*Csc[c + d*x])^(n - 2), x], x] /; FreeQ[{b, c, d}, x] && GtQ[n, 1] 
&& IntegerQ[2*n]
 

rule 4257
Int[csc[(c_.) + (d_.)*(x_)], x_Symbol] :> Simp[-ArcTanh[Cos[c + d*x]]/d, x] 
 /; FreeQ[{c, d}, x]
 
Maple [A] (verified)

Time = 13.39 (sec) , antiderivative size = 354, normalized size of antiderivative = 0.78

method result size
parts \(-\frac {A \,a^{4} \left (-\frac {16}{35}-\frac {\sec \left (d x +c \right )^{6}}{7}-\frac {6 \sec \left (d x +c \right )^{4}}{35}-\frac {8 \sec \left (d x +c \right )^{2}}{35}\right ) \tan \left (d x +c \right )}{d}+\frac {\left (4 A \,a^{3} b +B \,a^{4}\right ) \left (-\left (-\frac {\sec \left (d x +c \right )^{5}}{6}-\frac {5 \sec \left (d x +c \right )^{3}}{24}-\frac {5 \sec \left (d x +c \right )}{16}\right ) \tan \left (d x +c \right )+\frac {5 \ln \left (\sec \left (d x +c \right )+\tan \left (d x +c \right )\right )}{16}\right )}{d}+\frac {\left (B \,b^{4}+4 C a \,b^{3}\right ) \left (\frac {\sec \left (d x +c \right ) \tan \left (d x +c \right )}{2}+\frac {\ln \left (\sec \left (d x +c \right )+\tan \left (d x +c \right )\right )}{2}\right )}{d}-\frac {\left (A \,b^{4}+4 B a \,b^{3}+6 C \,a^{2} b^{2}\right ) \left (-\frac {2}{3}-\frac {\sec \left (d x +c \right )^{2}}{3}\right ) \tan \left (d x +c \right )}{d}+\frac {\left (4 a A \,b^{3}+6 B \,a^{2} b^{2}+4 a^{3} b C \right ) \left (-\left (-\frac {\sec \left (d x +c \right )^{3}}{4}-\frac {3 \sec \left (d x +c \right )}{8}\right ) \tan \left (d x +c \right )+\frac {3 \ln \left (\sec \left (d x +c \right )+\tan \left (d x +c \right )\right )}{8}\right )}{d}-\frac {\left (6 A \,a^{2} b^{2}+4 B \,a^{3} b +a^{4} C \right ) \left (-\frac {8}{15}-\frac {\sec \left (d x +c \right )^{4}}{5}-\frac {4 \sec \left (d x +c \right )^{2}}{15}\right ) \tan \left (d x +c \right )}{d}+\frac {C \,b^{4} \tan \left (d x +c \right )}{d}\) \(354\)
derivativedivides \(\frac {-A \,a^{4} \left (-\frac {16}{35}-\frac {\sec \left (d x +c \right )^{6}}{7}-\frac {6 \sec \left (d x +c \right )^{4}}{35}-\frac {8 \sec \left (d x +c \right )^{2}}{35}\right ) \tan \left (d x +c \right )+B \,a^{4} \left (-\left (-\frac {\sec \left (d x +c \right )^{5}}{6}-\frac {5 \sec \left (d x +c \right )^{3}}{24}-\frac {5 \sec \left (d x +c \right )}{16}\right ) \tan \left (d x +c \right )+\frac {5 \ln \left (\sec \left (d x +c \right )+\tan \left (d x +c \right )\right )}{16}\right )-a^{4} C \left (-\frac {8}{15}-\frac {\sec \left (d x +c \right )^{4}}{5}-\frac {4 \sec \left (d x +c \right )^{2}}{15}\right ) \tan \left (d x +c \right )+4 A \,a^{3} b \left (-\left (-\frac {\sec \left (d x +c \right )^{5}}{6}-\frac {5 \sec \left (d x +c \right )^{3}}{24}-\frac {5 \sec \left (d x +c \right )}{16}\right ) \tan \left (d x +c \right )+\frac {5 \ln \left (\sec \left (d x +c \right )+\tan \left (d x +c \right )\right )}{16}\right )-4 B \,a^{3} b \left (-\frac {8}{15}-\frac {\sec \left (d x +c \right )^{4}}{5}-\frac {4 \sec \left (d x +c \right )^{2}}{15}\right ) \tan \left (d x +c \right )+4 a^{3} b C \left (-\left (-\frac {\sec \left (d x +c \right )^{3}}{4}-\frac {3 \sec \left (d x +c \right )}{8}\right ) \tan \left (d x +c \right )+\frac {3 \ln \left (\sec \left (d x +c \right )+\tan \left (d x +c \right )\right )}{8}\right )-6 A \,a^{2} b^{2} \left (-\frac {8}{15}-\frac {\sec \left (d x +c \right )^{4}}{5}-\frac {4 \sec \left (d x +c \right )^{2}}{15}\right ) \tan \left (d x +c \right )+6 B \,a^{2} b^{2} \left (-\left (-\frac {\sec \left (d x +c \right )^{3}}{4}-\frac {3 \sec \left (d x +c \right )}{8}\right ) \tan \left (d x +c \right )+\frac {3 \ln \left (\sec \left (d x +c \right )+\tan \left (d x +c \right )\right )}{8}\right )-6 C \,a^{2} b^{2} \left (-\frac {2}{3}-\frac {\sec \left (d x +c \right )^{2}}{3}\right ) \tan \left (d x +c \right )+4 a A \,b^{3} \left (-\left (-\frac {\sec \left (d x +c \right )^{3}}{4}-\frac {3 \sec \left (d x +c \right )}{8}\right ) \tan \left (d x +c \right )+\frac {3 \ln \left (\sec \left (d x +c \right )+\tan \left (d x +c \right )\right )}{8}\right )-4 B a \,b^{3} \left (-\frac {2}{3}-\frac {\sec \left (d x +c \right )^{2}}{3}\right ) \tan \left (d x +c \right )+4 C a \,b^{3} \left (\frac {\sec \left (d x +c \right ) \tan \left (d x +c \right )}{2}+\frac {\ln \left (\sec \left (d x +c \right )+\tan \left (d x +c \right )\right )}{2}\right )-A \,b^{4} \left (-\frac {2}{3}-\frac {\sec \left (d x +c \right )^{2}}{3}\right ) \tan \left (d x +c \right )+B \,b^{4} \left (\frac {\sec \left (d x +c \right ) \tan \left (d x +c \right )}{2}+\frac {\ln \left (\sec \left (d x +c \right )+\tan \left (d x +c \right )\right )}{2}\right )+C \tan \left (d x +c \right ) b^{4}}{d}\) \(592\)
default \(\frac {-A \,a^{4} \left (-\frac {16}{35}-\frac {\sec \left (d x +c \right )^{6}}{7}-\frac {6 \sec \left (d x +c \right )^{4}}{35}-\frac {8 \sec \left (d x +c \right )^{2}}{35}\right ) \tan \left (d x +c \right )+B \,a^{4} \left (-\left (-\frac {\sec \left (d x +c \right )^{5}}{6}-\frac {5 \sec \left (d x +c \right )^{3}}{24}-\frac {5 \sec \left (d x +c \right )}{16}\right ) \tan \left (d x +c \right )+\frac {5 \ln \left (\sec \left (d x +c \right )+\tan \left (d x +c \right )\right )}{16}\right )-a^{4} C \left (-\frac {8}{15}-\frac {\sec \left (d x +c \right )^{4}}{5}-\frac {4 \sec \left (d x +c \right )^{2}}{15}\right ) \tan \left (d x +c \right )+4 A \,a^{3} b \left (-\left (-\frac {\sec \left (d x +c \right )^{5}}{6}-\frac {5 \sec \left (d x +c \right )^{3}}{24}-\frac {5 \sec \left (d x +c \right )}{16}\right ) \tan \left (d x +c \right )+\frac {5 \ln \left (\sec \left (d x +c \right )+\tan \left (d x +c \right )\right )}{16}\right )-4 B \,a^{3} b \left (-\frac {8}{15}-\frac {\sec \left (d x +c \right )^{4}}{5}-\frac {4 \sec \left (d x +c \right )^{2}}{15}\right ) \tan \left (d x +c \right )+4 a^{3} b C \left (-\left (-\frac {\sec \left (d x +c \right )^{3}}{4}-\frac {3 \sec \left (d x +c \right )}{8}\right ) \tan \left (d x +c \right )+\frac {3 \ln \left (\sec \left (d x +c \right )+\tan \left (d x +c \right )\right )}{8}\right )-6 A \,a^{2} b^{2} \left (-\frac {8}{15}-\frac {\sec \left (d x +c \right )^{4}}{5}-\frac {4 \sec \left (d x +c \right )^{2}}{15}\right ) \tan \left (d x +c \right )+6 B \,a^{2} b^{2} \left (-\left (-\frac {\sec \left (d x +c \right )^{3}}{4}-\frac {3 \sec \left (d x +c \right )}{8}\right ) \tan \left (d x +c \right )+\frac {3 \ln \left (\sec \left (d x +c \right )+\tan \left (d x +c \right )\right )}{8}\right )-6 C \,a^{2} b^{2} \left (-\frac {2}{3}-\frac {\sec \left (d x +c \right )^{2}}{3}\right ) \tan \left (d x +c \right )+4 a A \,b^{3} \left (-\left (-\frac {\sec \left (d x +c \right )^{3}}{4}-\frac {3 \sec \left (d x +c \right )}{8}\right ) \tan \left (d x +c \right )+\frac {3 \ln \left (\sec \left (d x +c \right )+\tan \left (d x +c \right )\right )}{8}\right )-4 B a \,b^{3} \left (-\frac {2}{3}-\frac {\sec \left (d x +c \right )^{2}}{3}\right ) \tan \left (d x +c \right )+4 C a \,b^{3} \left (\frac {\sec \left (d x +c \right ) \tan \left (d x +c \right )}{2}+\frac {\ln \left (\sec \left (d x +c \right )+\tan \left (d x +c \right )\right )}{2}\right )-A \,b^{4} \left (-\frac {2}{3}-\frac {\sec \left (d x +c \right )^{2}}{3}\right ) \tan \left (d x +c \right )+B \,b^{4} \left (\frac {\sec \left (d x +c \right ) \tan \left (d x +c \right )}{2}+\frac {\ln \left (\sec \left (d x +c \right )+\tan \left (d x +c \right )\right )}{2}\right )+C \tan \left (d x +c \right ) b^{4}}{d}\) \(592\)
parallelrisch \(\frac {-44100 \left (\frac {\cos \left (7 d x +7 c \right )}{21}+\frac {\cos \left (5 d x +5 c \right )}{3}+\cos \left (3 d x +3 c \right )+\frac {5 \cos \left (d x +c \right )}{3}\right ) \left (\frac {B \,a^{4}}{4}+b \left (A +\frac {6 C}{5}\right ) a^{3}+\frac {9 B \,a^{2} b^{2}}{5}+\frac {6 b^{3} \left (A +\frac {4 C}{3}\right ) a}{5}+\frac {2 B \,b^{4}}{5}\right ) \ln \left (\tan \left (\frac {d x}{2}+\frac {c}{2}\right )-1\right )+44100 \left (\frac {\cos \left (7 d x +7 c \right )}{21}+\frac {\cos \left (5 d x +5 c \right )}{3}+\cos \left (3 d x +3 c \right )+\frac {5 \cos \left (d x +c \right )}{3}\right ) \left (\frac {B \,a^{4}}{4}+b \left (A +\frac {6 C}{5}\right ) a^{3}+\frac {9 B \,a^{2} b^{2}}{5}+\frac {6 b^{3} \left (A +\frac {4 C}{3}\right ) a}{5}+\frac {2 B \,b^{4}}{5}\right ) \ln \left (\tan \left (\frac {d x}{2}+\frac {c}{2}\right )+1\right )+\left (\left (16128 A +18816 C \right ) a^{4}+75264 B \,a^{3} b +112896 \left (A +\frac {25 C}{28}\right ) b^{2} a^{2}+67200 B a \,b^{3}+16800 b^{4} \left (A +\frac {9 C}{10}\right )\right ) \sin \left (3 d x +3 c \right )+\left (\left (5376 A +6272 C \right ) a^{4}+25088 B \,a^{3} b +37632 \left (A +\frac {5 C}{4}\right ) a^{2} b^{2}+31360 B a \,b^{3}+7840 b^{4} \left (A +\frac {15 C}{14}\right )\right ) \sin \left (5 d x +5 c \right )+\left (\left (768 A +896 C \right ) a^{4}+3584 B \,a^{3} b +5376 \left (A +\frac {5 C}{4}\right ) a^{2} b^{2}+4480 B a \,b^{3}+1120 b^{4} \left (A +\frac {3 C}{2}\right )\right ) \sin \left (7 d x +7 c \right )+\left (19810 B \,a^{4}+79240 b \left (A +\frac {186 C}{283}\right ) a^{3}+78120 B \,a^{2} b^{2}+52080 \left (A +\frac {20 C}{31}\right ) b^{3} a +8400 B \,b^{4}\right ) \sin \left (2 d x +2 c \right )+\left (7000 B \,a^{4}+28000 b \left (A +\frac {6 C}{5}\right ) a^{3}+50400 B \,a^{2} b^{2}+33600 \left (A +\frac {4 C}{5}\right ) b^{3} a +6720 B \,b^{4}\right ) \sin \left (4 d x +4 c \right )+\left (1050 B \,a^{4}+4200 b \left (A +\frac {6 C}{5}\right ) a^{3}+7560 B \,a^{2} b^{2}+5040 b^{3} \left (A +\frac {4 C}{3}\right ) a +1680 B \,b^{4}\right ) \sin \left (6 d x +6 c \right )+26880 \sin \left (d x +c \right ) \left (a^{4} \left (A +\frac {C}{2}\right )+2 B \,a^{3} b +3 \left (A +\frac {3 C}{4}\right ) a^{2} b^{2}+\frac {3 B a \,b^{3}}{2}+\frac {3 b^{4} \left (A +\frac {5 C}{6}\right )}{8}\right )}{1680 d \left (\cos \left (7 d x +7 c \right )+7 \cos \left (5 d x +5 c \right )+21 \cos \left (3 d x +3 c \right )+35 \cos \left (d x +c \right )\right )}\) \(637\)
risch \(\text {Expression too large to display}\) \(1648\)

Input:

int((a+b*cos(d*x+c))^4*(A+B*cos(d*x+c)+C*cos(d*x+c)^2)*sec(d*x+c)^8,x,meth 
od=_RETURNVERBOSE)
 

Output:

-A*a^4/d*(-16/35-1/7*sec(d*x+c)^6-6/35*sec(d*x+c)^4-8/35*sec(d*x+c)^2)*tan 
(d*x+c)+(4*A*a^3*b+B*a^4)/d*(-(-1/6*sec(d*x+c)^5-5/24*sec(d*x+c)^3-5/16*se 
c(d*x+c))*tan(d*x+c)+5/16*ln(sec(d*x+c)+tan(d*x+c)))+(B*b^4+4*C*a*b^3)/d*( 
1/2*sec(d*x+c)*tan(d*x+c)+1/2*ln(sec(d*x+c)+tan(d*x+c)))-(A*b^4+4*B*a*b^3+ 
6*C*a^2*b^2)/d*(-2/3-1/3*sec(d*x+c)^2)*tan(d*x+c)+(4*A*a*b^3+6*B*a^2*b^2+4 
*C*a^3*b)/d*(-(-1/4*sec(d*x+c)^3-3/8*sec(d*x+c))*tan(d*x+c)+3/8*ln(sec(d*x 
+c)+tan(d*x+c)))-(6*A*a^2*b^2+4*B*a^3*b+C*a^4)/d*(-8/15-1/5*sec(d*x+c)^4-4 
/15*sec(d*x+c)^2)*tan(d*x+c)+C*b^4/d*tan(d*x+c)
 

Fricas [A] (verification not implemented)

Time = 0.12 (sec) , antiderivative size = 450, normalized size of antiderivative = 0.99 \[ \int (a+b \cos (c+d x))^4 \left (A+B \cos (c+d x)+C \cos ^2(c+d x)\right ) \sec ^8(c+d x) \, dx=\frac {105 \, {\left (5 \, B a^{4} + 4 \, {\left (5 \, A + 6 \, C\right )} a^{3} b + 36 \, B a^{2} b^{2} + 8 \, {\left (3 \, A + 4 \, C\right )} a b^{3} + 8 \, B b^{4}\right )} \cos \left (d x + c\right )^{7} \log \left (\sin \left (d x + c\right ) + 1\right ) - 105 \, {\left (5 \, B a^{4} + 4 \, {\left (5 \, A + 6 \, C\right )} a^{3} b + 36 \, B a^{2} b^{2} + 8 \, {\left (3 \, A + 4 \, C\right )} a b^{3} + 8 \, B b^{4}\right )} \cos \left (d x + c\right )^{7} \log \left (-\sin \left (d x + c\right ) + 1\right ) + 2 \, {\left (16 \, {\left (8 \, {\left (6 \, A + 7 \, C\right )} a^{4} + 224 \, B a^{3} b + 84 \, {\left (4 \, A + 5 \, C\right )} a^{2} b^{2} + 280 \, B a b^{3} + 35 \, {\left (2 \, A + 3 \, C\right )} b^{4}\right )} \cos \left (d x + c\right )^{6} + 105 \, {\left (5 \, B a^{4} + 4 \, {\left (5 \, A + 6 \, C\right )} a^{3} b + 36 \, B a^{2} b^{2} + 8 \, {\left (3 \, A + 4 \, C\right )} a b^{3} + 8 \, B b^{4}\right )} \cos \left (d x + c\right )^{5} + 240 \, A a^{4} + 16 \, {\left (4 \, {\left (6 \, A + 7 \, C\right )} a^{4} + 112 \, B a^{3} b + 42 \, {\left (4 \, A + 5 \, C\right )} a^{2} b^{2} + 140 \, B a b^{3} + 35 \, A b^{4}\right )} \cos \left (d x + c\right )^{4} + 70 \, {\left (5 \, B a^{4} + 4 \, {\left (5 \, A + 6 \, C\right )} a^{3} b + 36 \, B a^{2} b^{2} + 24 \, A a b^{3}\right )} \cos \left (d x + c\right )^{3} + 48 \, {\left ({\left (6 \, A + 7 \, C\right )} a^{4} + 28 \, B a^{3} b + 42 \, A a^{2} b^{2}\right )} \cos \left (d x + c\right )^{2} + 280 \, {\left (B a^{4} + 4 \, A a^{3} b\right )} \cos \left (d x + c\right )\right )} \sin \left (d x + c\right )}{3360 \, d \cos \left (d x + c\right )^{7}} \] Input:

integrate((a+b*cos(d*x+c))^4*(A+B*cos(d*x+c)+C*cos(d*x+c)^2)*sec(d*x+c)^8, 
x, algorithm="fricas")
 

Output:

1/3360*(105*(5*B*a^4 + 4*(5*A + 6*C)*a^3*b + 36*B*a^2*b^2 + 8*(3*A + 4*C)* 
a*b^3 + 8*B*b^4)*cos(d*x + c)^7*log(sin(d*x + c) + 1) - 105*(5*B*a^4 + 4*( 
5*A + 6*C)*a^3*b + 36*B*a^2*b^2 + 8*(3*A + 4*C)*a*b^3 + 8*B*b^4)*cos(d*x + 
 c)^7*log(-sin(d*x + c) + 1) + 2*(16*(8*(6*A + 7*C)*a^4 + 224*B*a^3*b + 84 
*(4*A + 5*C)*a^2*b^2 + 280*B*a*b^3 + 35*(2*A + 3*C)*b^4)*cos(d*x + c)^6 + 
105*(5*B*a^4 + 4*(5*A + 6*C)*a^3*b + 36*B*a^2*b^2 + 8*(3*A + 4*C)*a*b^3 + 
8*B*b^4)*cos(d*x + c)^5 + 240*A*a^4 + 16*(4*(6*A + 7*C)*a^4 + 112*B*a^3*b 
+ 42*(4*A + 5*C)*a^2*b^2 + 140*B*a*b^3 + 35*A*b^4)*cos(d*x + c)^4 + 70*(5* 
B*a^4 + 4*(5*A + 6*C)*a^3*b + 36*B*a^2*b^2 + 24*A*a*b^3)*cos(d*x + c)^3 + 
48*((6*A + 7*C)*a^4 + 28*B*a^3*b + 42*A*a^2*b^2)*cos(d*x + c)^2 + 280*(B*a 
^4 + 4*A*a^3*b)*cos(d*x + c))*sin(d*x + c))/(d*cos(d*x + c)^7)
 

Sympy [F(-1)]

Timed out. \[ \int (a+b \cos (c+d x))^4 \left (A+B \cos (c+d x)+C \cos ^2(c+d x)\right ) \sec ^8(c+d x) \, dx=\text {Timed out} \] Input:

integrate((a+b*cos(d*x+c))**4*(A+B*cos(d*x+c)+C*cos(d*x+c)**2)*sec(d*x+c)* 
*8,x)
 

Output:

Timed out
 

Maxima [A] (verification not implemented)

Time = 0.05 (sec) , antiderivative size = 746, normalized size of antiderivative = 1.64 \[ \int (a+b \cos (c+d x))^4 \left (A+B \cos (c+d x)+C \cos ^2(c+d x)\right ) \sec ^8(c+d x) \, dx=\text {Too large to display} \] Input:

integrate((a+b*cos(d*x+c))^4*(A+B*cos(d*x+c)+C*cos(d*x+c)^2)*sec(d*x+c)^8, 
x, algorithm="maxima")
 

Output:

1/3360*(96*(5*tan(d*x + c)^7 + 21*tan(d*x + c)^5 + 35*tan(d*x + c)^3 + 35* 
tan(d*x + c))*A*a^4 + 224*(3*tan(d*x + c)^5 + 10*tan(d*x + c)^3 + 15*tan(d 
*x + c))*C*a^4 + 896*(3*tan(d*x + c)^5 + 10*tan(d*x + c)^3 + 15*tan(d*x + 
c))*B*a^3*b + 1344*(3*tan(d*x + c)^5 + 10*tan(d*x + c)^3 + 15*tan(d*x + c) 
)*A*a^2*b^2 + 6720*(tan(d*x + c)^3 + 3*tan(d*x + c))*C*a^2*b^2 + 4480*(tan 
(d*x + c)^3 + 3*tan(d*x + c))*B*a*b^3 + 1120*(tan(d*x + c)^3 + 3*tan(d*x + 
 c))*A*b^4 - 35*B*a^4*(2*(15*sin(d*x + c)^5 - 40*sin(d*x + c)^3 + 33*sin(d 
*x + c))/(sin(d*x + c)^6 - 3*sin(d*x + c)^4 + 3*sin(d*x + c)^2 - 1) - 15*l 
og(sin(d*x + c) + 1) + 15*log(sin(d*x + c) - 1)) - 140*A*a^3*b*(2*(15*sin( 
d*x + c)^5 - 40*sin(d*x + c)^3 + 33*sin(d*x + c))/(sin(d*x + c)^6 - 3*sin( 
d*x + c)^4 + 3*sin(d*x + c)^2 - 1) - 15*log(sin(d*x + c) + 1) + 15*log(sin 
(d*x + c) - 1)) - 840*C*a^3*b*(2*(3*sin(d*x + c)^3 - 5*sin(d*x + c))/(sin( 
d*x + c)^4 - 2*sin(d*x + c)^2 + 1) - 3*log(sin(d*x + c) + 1) + 3*log(sin(d 
*x + c) - 1)) - 1260*B*a^2*b^2*(2*(3*sin(d*x + c)^3 - 5*sin(d*x + c))/(sin 
(d*x + c)^4 - 2*sin(d*x + c)^2 + 1) - 3*log(sin(d*x + c) + 1) + 3*log(sin( 
d*x + c) - 1)) - 840*A*a*b^3*(2*(3*sin(d*x + c)^3 - 5*sin(d*x + c))/(sin(d 
*x + c)^4 - 2*sin(d*x + c)^2 + 1) - 3*log(sin(d*x + c) + 1) + 3*log(sin(d* 
x + c) - 1)) - 3360*C*a*b^3*(2*sin(d*x + c)/(sin(d*x + c)^2 - 1) - log(sin 
(d*x + c) + 1) + log(sin(d*x + c) - 1)) - 840*B*b^4*(2*sin(d*x + c)/(sin(d 
*x + c)^2 - 1) - log(sin(d*x + c) + 1) + log(sin(d*x + c) - 1)) + 3360*...
 

Giac [B] (verification not implemented)

Leaf count of result is larger than twice the leaf count of optimal. 1888 vs. \(2 (439) = 878\).

Time = 0.24 (sec) , antiderivative size = 1888, normalized size of antiderivative = 4.16 \[ \int (a+b \cos (c+d x))^4 \left (A+B \cos (c+d x)+C \cos ^2(c+d x)\right ) \sec ^8(c+d x) \, dx=\text {Too large to display} \] Input:

integrate((a+b*cos(d*x+c))^4*(A+B*cos(d*x+c)+C*cos(d*x+c)^2)*sec(d*x+c)^8, 
x, algorithm="giac")
 

Output:

1/1680*(105*(5*B*a^4 + 20*A*a^3*b + 24*C*a^3*b + 36*B*a^2*b^2 + 24*A*a*b^3 
 + 32*C*a*b^3 + 8*B*b^4)*log(abs(tan(1/2*d*x + 1/2*c) + 1)) - 105*(5*B*a^4 
 + 20*A*a^3*b + 24*C*a^3*b + 36*B*a^2*b^2 + 24*A*a*b^3 + 32*C*a*b^3 + 8*B* 
b^4)*log(abs(tan(1/2*d*x + 1/2*c) - 1)) - 2*(1680*A*a^4*tan(1/2*d*x + 1/2* 
c)^13 - 1155*B*a^4*tan(1/2*d*x + 1/2*c)^13 + 1680*C*a^4*tan(1/2*d*x + 1/2* 
c)^13 - 4620*A*a^3*b*tan(1/2*d*x + 1/2*c)^13 + 6720*B*a^3*b*tan(1/2*d*x + 
1/2*c)^13 - 4200*C*a^3*b*tan(1/2*d*x + 1/2*c)^13 + 10080*A*a^2*b^2*tan(1/2 
*d*x + 1/2*c)^13 - 6300*B*a^2*b^2*tan(1/2*d*x + 1/2*c)^13 + 10080*C*a^2*b^ 
2*tan(1/2*d*x + 1/2*c)^13 - 4200*A*a*b^3*tan(1/2*d*x + 1/2*c)^13 + 6720*B* 
a*b^3*tan(1/2*d*x + 1/2*c)^13 - 3360*C*a*b^3*tan(1/2*d*x + 1/2*c)^13 + 168 
0*A*b^4*tan(1/2*d*x + 1/2*c)^13 - 840*B*b^4*tan(1/2*d*x + 1/2*c)^13 + 1680 
*C*b^4*tan(1/2*d*x + 1/2*c)^13 - 3360*A*a^4*tan(1/2*d*x + 1/2*c)^11 + 980* 
B*a^4*tan(1/2*d*x + 1/2*c)^11 - 5600*C*a^4*tan(1/2*d*x + 1/2*c)^11 + 3920* 
A*a^3*b*tan(1/2*d*x + 1/2*c)^11 - 22400*B*a^3*b*tan(1/2*d*x + 1/2*c)^11 + 
10080*C*a^3*b*tan(1/2*d*x + 1/2*c)^11 - 33600*A*a^2*b^2*tan(1/2*d*x + 1/2* 
c)^11 + 15120*B*a^2*b^2*tan(1/2*d*x + 1/2*c)^11 - 47040*C*a^2*b^2*tan(1/2* 
d*x + 1/2*c)^11 + 10080*A*a*b^3*tan(1/2*d*x + 1/2*c)^11 - 31360*B*a*b^3*ta 
n(1/2*d*x + 1/2*c)^11 + 13440*C*a*b^3*tan(1/2*d*x + 1/2*c)^11 - 7840*A*b^4 
*tan(1/2*d*x + 1/2*c)^11 + 3360*B*b^4*tan(1/2*d*x + 1/2*c)^11 - 10080*C*b^ 
4*tan(1/2*d*x + 1/2*c)^11 + 14448*A*a^4*tan(1/2*d*x + 1/2*c)^9 - 2975*B...
 

Mupad [B] (verification not implemented)

Time = 2.25 (sec) , antiderivative size = 1044, normalized size of antiderivative = 2.30 \[ \int (a+b \cos (c+d x))^4 \left (A+B \cos (c+d x)+C \cos ^2(c+d x)\right ) \sec ^8(c+d x) \, dx=\text {Too large to display} \] Input:

int(((a + b*cos(c + d*x))^4*(A + B*cos(c + d*x) + C*cos(c + d*x)^2))/cos(c 
 + d*x)^8,x)
 

Output:

(atanh((4*tan(c/2 + (d*x)/2)*((5*B*a^4)/16 + (B*b^4)/2 + (9*B*a^2*b^2)/4 + 
 (3*A*a*b^3)/2 + (5*A*a^3*b)/4 + 2*C*a*b^3 + (3*C*a^3*b)/2))/((5*B*a^4)/4 
+ 2*B*b^4 + 9*B*a^2*b^2 + 6*A*a*b^3 + 5*A*a^3*b + 8*C*a*b^3 + 6*C*a^3*b))* 
((5*B*a^4)/8 + B*b^4 + (9*B*a^2*b^2)/2 + 3*A*a*b^3 + (5*A*a^3*b)/2 + 4*C*a 
*b^3 + 3*C*a^3*b))/d - (tan(c/2 + (d*x)/2)^13*(2*A*a^4 + 2*A*b^4 - (11*B*a 
^4)/8 - B*b^4 + 2*C*a^4 + 2*C*b^4 + 12*A*a^2*b^2 - (15*B*a^2*b^2)/2 + 12*C 
*a^2*b^2 - 5*A*a*b^3 - (11*A*a^3*b)/2 + 8*B*a*b^3 + 8*B*a^3*b - 4*C*a*b^3 
- 5*C*a^3*b) - tan(c/2 + (d*x)/2)^3*(4*A*a^4 + (28*A*b^4)/3 + (7*B*a^4)/6 
+ 4*B*b^4 + (20*C*a^4)/3 + 12*C*b^4 + 40*A*a^2*b^2 + 18*B*a^2*b^2 + 56*C*a 
^2*b^2 + 12*A*a*b^3 + (14*A*a^3*b)/3 + (112*B*a*b^3)/3 + (80*B*a^3*b)/3 + 
16*C*a*b^3 + 12*C*a^3*b) - tan(c/2 + (d*x)/2)^11*(4*A*a^4 + (28*A*b^4)/3 - 
 (7*B*a^4)/6 - 4*B*b^4 + (20*C*a^4)/3 + 12*C*b^4 + 40*A*a^2*b^2 - 18*B*a^2 
*b^2 + 56*C*a^2*b^2 - 12*A*a*b^3 - (14*A*a^3*b)/3 + (112*B*a*b^3)/3 + (80* 
B*a^3*b)/3 - 16*C*a*b^3 - 12*C*a^3*b) + tan(c/2 + (d*x)/2)^5*((86*A*a^4)/5 
 + (58*A*b^4)/3 + (85*B*a^4)/24 + 5*B*b^4 + (226*C*a^4)/15 + 30*C*b^4 + (4 
52*A*a^2*b^2)/5 + (27*B*a^2*b^2)/2 + 116*C*a^2*b^2 + 9*A*a*b^3 + (85*A*a^3 
*b)/6 + (232*B*a*b^3)/3 + (904*B*a^3*b)/15 + 20*C*a*b^3 + 9*C*a^3*b) + tan 
(c/2 + (d*x)/2)^9*((86*A*a^4)/5 + (58*A*b^4)/3 - (85*B*a^4)/24 - 5*B*b^4 + 
 (226*C*a^4)/15 + 30*C*b^4 + (452*A*a^2*b^2)/5 - (27*B*a^2*b^2)/2 + 116*C* 
a^2*b^2 - 9*A*a*b^3 - (85*A*a^3*b)/6 + (232*B*a*b^3)/3 + (904*B*a^3*b)/...
 

Reduce [B] (verification not implemented)

Time = 0.20 (sec) , antiderivative size = 1909, normalized size of antiderivative = 4.20 \[ \int (a+b \cos (c+d x))^4 \left (A+B \cos (c+d x)+C \cos ^2(c+d x)\right ) \sec ^8(c+d x) \, dx =\text {Too large to display} \] Input:

int((a+b*cos(d*x+c))^4*(A+B*cos(d*x+c)+C*cos(d*x+c)^2)*sec(d*x+c)^8,x)
 

Output:

( - 2625*cos(c + d*x)*log(tan((c + d*x)/2) - 1)*sin(c + d*x)**6*a**4*b - 2 
520*cos(c + d*x)*log(tan((c + d*x)/2) - 1)*sin(c + d*x)**6*a**3*b*c - 6300 
*cos(c + d*x)*log(tan((c + d*x)/2) - 1)*sin(c + d*x)**6*a**2*b**3 - 3360*c 
os(c + d*x)*log(tan((c + d*x)/2) - 1)*sin(c + d*x)**6*a*b**3*c - 840*cos(c 
 + d*x)*log(tan((c + d*x)/2) - 1)*sin(c + d*x)**6*b**5 + 7875*cos(c + d*x) 
*log(tan((c + d*x)/2) - 1)*sin(c + d*x)**4*a**4*b + 7560*cos(c + d*x)*log( 
tan((c + d*x)/2) - 1)*sin(c + d*x)**4*a**3*b*c + 18900*cos(c + d*x)*log(ta 
n((c + d*x)/2) - 1)*sin(c + d*x)**4*a**2*b**3 + 10080*cos(c + d*x)*log(tan 
((c + d*x)/2) - 1)*sin(c + d*x)**4*a*b**3*c + 2520*cos(c + d*x)*log(tan((c 
 + d*x)/2) - 1)*sin(c + d*x)**4*b**5 - 7875*cos(c + d*x)*log(tan((c + d*x) 
/2) - 1)*sin(c + d*x)**2*a**4*b - 7560*cos(c + d*x)*log(tan((c + d*x)/2) - 
 1)*sin(c + d*x)**2*a**3*b*c - 18900*cos(c + d*x)*log(tan((c + d*x)/2) - 1 
)*sin(c + d*x)**2*a**2*b**3 - 10080*cos(c + d*x)*log(tan((c + d*x)/2) - 1) 
*sin(c + d*x)**2*a*b**3*c - 2520*cos(c + d*x)*log(tan((c + d*x)/2) - 1)*si 
n(c + d*x)**2*b**5 + 2625*cos(c + d*x)*log(tan((c + d*x)/2) - 1)*a**4*b + 
2520*cos(c + d*x)*log(tan((c + d*x)/2) - 1)*a**3*b*c + 6300*cos(c + d*x)*l 
og(tan((c + d*x)/2) - 1)*a**2*b**3 + 3360*cos(c + d*x)*log(tan((c + d*x)/2 
) - 1)*a*b**3*c + 840*cos(c + d*x)*log(tan((c + d*x)/2) - 1)*b**5 + 2625*c 
os(c + d*x)*log(tan((c + d*x)/2) + 1)*sin(c + d*x)**6*a**4*b + 2520*cos(c 
+ d*x)*log(tan((c + d*x)/2) + 1)*sin(c + d*x)**6*a**3*b*c + 6300*cos(c ...