\(\int (b \cos (c+d x))^{3/2} (A+B \cos (c+d x)) \sec ^3(c+d x) \, dx\) [810]

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

Optimal result

Integrand size = 31, antiderivative size = 110 \[ \int (b \cos (c+d x))^{3/2} (A+B \cos (c+d x)) \sec ^3(c+d x) \, dx=-\frac {2 A b \sqrt {b \cos (c+d x)} E\left (\left .\frac {1}{2} (c+d x)\right |2\right )}{d \sqrt {\cos (c+d x)}}+\frac {2 b^2 B \sqrt {\cos (c+d x)} \operatorname {EllipticF}\left (\frac {1}{2} (c+d x),2\right )}{d \sqrt {b \cos (c+d x)}}+\frac {2 A b^2 \sin (c+d x)}{d \sqrt {b \cos (c+d x)}} \] Output:

-2*A*b*(b*cos(d*x+c))^(1/2)*EllipticE(sin(1/2*d*x+1/2*c),2^(1/2))/d/cos(d* 
x+c)^(1/2)+2*b^2*B*cos(d*x+c)^(1/2)*InverseJacobiAM(1/2*d*x+1/2*c,2^(1/2)) 
/d/(b*cos(d*x+c))^(1/2)+2*A*b^2*sin(d*x+c)/d/(b*cos(d*x+c))^(1/2)
 

Mathematica [A] (verified)

Time = 0.51 (sec) , antiderivative size = 73, normalized size of antiderivative = 0.66 \[ \int (b \cos (c+d x))^{3/2} (A+B \cos (c+d x)) \sec ^3(c+d x) \, dx=\frac {2 (b \cos (c+d x))^{3/2} \left (-A E\left (\left .\frac {1}{2} (c+d x)\right |2\right )+B \operatorname {EllipticF}\left (\frac {1}{2} (c+d x),2\right )+\frac {A \sin (c+d x)}{\sqrt {\cos (c+d x)}}\right )}{d \cos ^{\frac {3}{2}}(c+d x)} \] Input:

Integrate[(b*Cos[c + d*x])^(3/2)*(A + B*Cos[c + d*x])*Sec[c + d*x]^3,x]
 

Output:

(2*(b*Cos[c + d*x])^(3/2)*(-(A*EllipticE[(c + d*x)/2, 2]) + B*EllipticF[(c 
 + d*x)/2, 2] + (A*Sin[c + d*x])/Sqrt[Cos[c + d*x]]))/(d*Cos[c + d*x]^(3/2 
))
 

Rubi [A] (verified)

Time = 0.57 (sec) , antiderivative size = 117, normalized size of antiderivative = 1.06, number of steps used = 10, number of rules used = 10, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.323, Rules used = {3042, 2030, 3227, 3042, 3116, 3042, 3121, 3042, 3119, 3120}

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 ^3(c+d x) (b \cos (c+d x))^{3/2} (A+B \cos (c+d x)) \, dx\)

\(\Big \downarrow \) 3042

\(\displaystyle \int \frac {\left (b \sin \left (c+d x+\frac {\pi }{2}\right )\right )^{3/2} \left (A+B \sin \left (c+d x+\frac {\pi }{2}\right )\right )}{\sin \left (c+d x+\frac {\pi }{2}\right )^3}dx\)

\(\Big \downarrow \) 2030

\(\displaystyle b^3 \int \frac {A+B \sin \left (\frac {1}{2} (2 c+\pi )+d x\right )}{\left (b \sin \left (\frac {1}{2} (2 c+\pi )+d x\right )\right )^{3/2}}dx\)

\(\Big \downarrow \) 3227

\(\displaystyle b^3 \left (A \int \frac {1}{(b \cos (c+d x))^{3/2}}dx+\frac {B \int \frac {1}{\sqrt {b \cos (c+d x)}}dx}{b}\right )\)

\(\Big \downarrow \) 3042

\(\displaystyle b^3 \left (A \int \frac {1}{\left (b \sin \left (c+d x+\frac {\pi }{2}\right )\right )^{3/2}}dx+\frac {B \int \frac {1}{\sqrt {b \sin \left (c+d x+\frac {\pi }{2}\right )}}dx}{b}\right )\)

\(\Big \downarrow \) 3116

\(\displaystyle b^3 \left (A \left (\frac {2 \sin (c+d x)}{b d \sqrt {b \cos (c+d x)}}-\frac {\int \sqrt {b \cos (c+d x)}dx}{b^2}\right )+\frac {B \int \frac {1}{\sqrt {b \sin \left (c+d x+\frac {\pi }{2}\right )}}dx}{b}\right )\)

\(\Big \downarrow \) 3042

\(\displaystyle b^3 \left (A \left (\frac {2 \sin (c+d x)}{b d \sqrt {b \cos (c+d x)}}-\frac {\int \sqrt {b \sin \left (c+d x+\frac {\pi }{2}\right )}dx}{b^2}\right )+\frac {B \int \frac {1}{\sqrt {b \sin \left (c+d x+\frac {\pi }{2}\right )}}dx}{b}\right )\)

\(\Big \downarrow \) 3121

\(\displaystyle b^3 \left (A \left (\frac {2 \sin (c+d x)}{b d \sqrt {b \cos (c+d x)}}-\frac {\sqrt {b \cos (c+d x)} \int \sqrt {\cos (c+d x)}dx}{b^2 \sqrt {\cos (c+d x)}}\right )+\frac {B \sqrt {\cos (c+d x)} \int \frac {1}{\sqrt {\cos (c+d x)}}dx}{b \sqrt {b \cos (c+d x)}}\right )\)

\(\Big \downarrow \) 3042

\(\displaystyle b^3 \left (A \left (\frac {2 \sin (c+d x)}{b d \sqrt {b \cos (c+d x)}}-\frac {\sqrt {b \cos (c+d x)} \int \sqrt {\sin \left (c+d x+\frac {\pi }{2}\right )}dx}{b^2 \sqrt {\cos (c+d x)}}\right )+\frac {B \sqrt {\cos (c+d x)} \int \frac {1}{\sqrt {\sin \left (c+d x+\frac {\pi }{2}\right )}}dx}{b \sqrt {b \cos (c+d x)}}\right )\)

\(\Big \downarrow \) 3119

\(\displaystyle b^3 \left (\frac {B \sqrt {\cos (c+d x)} \int \frac {1}{\sqrt {\sin \left (c+d x+\frac {\pi }{2}\right )}}dx}{b \sqrt {b \cos (c+d x)}}+A \left (\frac {2 \sin (c+d x)}{b d \sqrt {b \cos (c+d x)}}-\frac {2 E\left (\left .\frac {1}{2} (c+d x)\right |2\right ) \sqrt {b \cos (c+d x)}}{b^2 d \sqrt {\cos (c+d x)}}\right )\right )\)

\(\Big \downarrow \) 3120

\(\displaystyle b^3 \left (A \left (\frac {2 \sin (c+d x)}{b d \sqrt {b \cos (c+d x)}}-\frac {2 E\left (\left .\frac {1}{2} (c+d x)\right |2\right ) \sqrt {b \cos (c+d x)}}{b^2 d \sqrt {\cos (c+d x)}}\right )+\frac {2 B \sqrt {\cos (c+d x)} \operatorname {EllipticF}\left (\frac {1}{2} (c+d x),2\right )}{b d \sqrt {b \cos (c+d x)}}\right )\)

Input:

Int[(b*Cos[c + d*x])^(3/2)*(A + B*Cos[c + d*x])*Sec[c + d*x]^3,x]
 

Output:

b^3*((2*B*Sqrt[Cos[c + d*x]]*EllipticF[(c + d*x)/2, 2])/(b*d*Sqrt[b*Cos[c 
+ d*x]]) + A*((-2*Sqrt[b*Cos[c + d*x]]*EllipticE[(c + d*x)/2, 2])/(b^2*d*S 
qrt[Cos[c + d*x]]) + (2*Sin[c + d*x])/(b*d*Sqrt[b*Cos[c + d*x]])))
 

Defintions of rubi rules used

rule 2030
Int[(Fx_.)*(v_)^(m_.)*((b_)*(v_))^(n_), x_Symbol] :> Simp[1/b^m   Int[(b*v) 
^(m + n)*Fx, x], x] /; FreeQ[{b, n}, x] && IntegerQ[m]
 

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

rule 3116
Int[((b_.)*sin[(c_.) + (d_.)*(x_)])^(n_), x_Symbol] :> Simp[Cos[c + d*x]*(( 
b*Sin[c + d*x])^(n + 1)/(b*d*(n + 1))), x] + Simp[(n + 2)/(b^2*(n + 1))   I 
nt[(b*Sin[c + d*x])^(n + 2), x], x] /; FreeQ[{b, c, d}, x] && LtQ[n, -1] && 
 IntegerQ[2*n]
 

rule 3119
Int[Sqrt[sin[(c_.) + (d_.)*(x_)]], x_Symbol] :> Simp[(2/d)*EllipticE[(1/2)* 
(c - Pi/2 + d*x), 2], x] /; FreeQ[{c, d}, x]
 

rule 3120
Int[1/Sqrt[sin[(c_.) + (d_.)*(x_)]], x_Symbol] :> Simp[(2/d)*EllipticF[(1/2 
)*(c - Pi/2 + d*x), 2], x] /; FreeQ[{c, d}, x]
 

rule 3121
Int[((b_)*sin[(c_.) + (d_.)*(x_)])^(n_), x_Symbol] :> Simp[(b*Sin[c + d*x]) 
^n/Sin[c + d*x]^n   Int[Sin[c + d*x]^n, x], x] /; FreeQ[{b, c, d}, x] && Lt 
Q[-1, n, 1] && IntegerQ[2*n]
 

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]
 
Maple [B] (verified)

Leaf count of result is larger than twice the leaf count of optimal. \(216\) vs. \(2(103)=206\).

Time = 8.17 (sec) , antiderivative size = 217, normalized size of antiderivative = 1.97

method result size
default \(\frac {2 b^{2} \sqrt {-2 b \sin \left (\frac {d x}{2}+\frac {c}{2}\right )^{4}+\sin \left (\frac {d x}{2}+\frac {c}{2}\right )^{2} b}\, \left (2 A \cos \left (\frac {d x}{2}+\frac {c}{2}\right ) \sin \left (\frac {d x}{2}+\frac {c}{2}\right )^{2}-A \sqrt {\frac {1}{2}-\frac {\cos \left (d x +c \right )}{2}}\, \sqrt {2 \sin \left (\frac {d x}{2}+\frac {c}{2}\right )^{2}-1}\, \operatorname {EllipticE}\left (\cos \left (\frac {d x}{2}+\frac {c}{2}\right ), \sqrt {2}\right )-B \sqrt {\frac {1}{2}-\frac {\cos \left (d x +c \right )}{2}}\, \sqrt {2 \sin \left (\frac {d x}{2}+\frac {c}{2}\right )^{2}-1}\, \operatorname {EllipticF}\left (\cos \left (\frac {d x}{2}+\frac {c}{2}\right ), \sqrt {2}\right )\right )}{\sqrt {-b \left (2 \sin \left (\frac {d x}{2}+\frac {c}{2}\right )^{4}-\sin \left (\frac {d x}{2}+\frac {c}{2}\right )^{2}\right )}\, \sin \left (\frac {d x}{2}+\frac {c}{2}\right ) \sqrt {\left (2 \cos \left (\frac {d x}{2}+\frac {c}{2}\right )^{2}-1\right ) b}\, d}\) \(217\)
parts \(-\frac {2 A \,b^{2} \left (-2 \cos \left (\frac {d x}{2}+\frac {c}{2}\right ) \sqrt {-2 b \sin \left (\frac {d x}{2}+\frac {c}{2}\right )^{4}+\sin \left (\frac {d x}{2}+\frac {c}{2}\right )^{2} b}\, \sin \left (\frac {d x}{2}+\frac {c}{2}\right )^{2}+\sqrt {\frac {1}{2}-\frac {\cos \left (d x +c \right )}{2}}\, \sqrt {2 \sin \left (\frac {d x}{2}+\frac {c}{2}\right )^{2}-1}\, \sqrt {-2 b \sin \left (\frac {d x}{2}+\frac {c}{2}\right )^{4}+\sin \left (\frac {d x}{2}+\frac {c}{2}\right )^{2} b}\, \operatorname {EllipticE}\left (\cos \left (\frac {d x}{2}+\frac {c}{2}\right ), \sqrt {2}\right )\right )}{\sqrt {-b \left (2 \sin \left (\frac {d x}{2}+\frac {c}{2}\right )^{4}-\sin \left (\frac {d x}{2}+\frac {c}{2}\right )^{2}\right )}\, \sin \left (\frac {d x}{2}+\frac {c}{2}\right ) \sqrt {\left (2 \cos \left (\frac {d x}{2}+\frac {c}{2}\right )^{2}-1\right ) b}\, d}-\frac {2 B \sqrt {\left (2 \cos \left (\frac {d x}{2}+\frac {c}{2}\right )^{2}-1\right ) b \sin \left (\frac {d x}{2}+\frac {c}{2}\right )^{2}}\, b^{2} \sqrt {\frac {1}{2}-\frac {\cos \left (d x +c \right )}{2}}\, \sqrt {-2 \cos \left (\frac {d x}{2}+\frac {c}{2}\right )^{2}+1}\, \operatorname {EllipticF}\left (\cos \left (\frac {d x}{2}+\frac {c}{2}\right ), \sqrt {2}\right )}{\sqrt {-b \left (2 \sin \left (\frac {d x}{2}+\frac {c}{2}\right )^{4}-\sin \left (\frac {d x}{2}+\frac {c}{2}\right )^{2}\right )}\, \sin \left (\frac {d x}{2}+\frac {c}{2}\right ) \sqrt {\left (2 \cos \left (\frac {d x}{2}+\frac {c}{2}\right )^{2}-1\right ) b}\, d}\) \(344\)

Input:

int((cos(d*x+c)*b)^(3/2)*(A+B*cos(d*x+c))*sec(d*x+c)^3,x,method=_RETURNVER 
BOSE)
 

Output:

2*b^2*(-2*b*sin(1/2*d*x+1/2*c)^4+sin(1/2*d*x+1/2*c)^2*b)^(1/2)*(2*A*cos(1/ 
2*d*x+1/2*c)*sin(1/2*d*x+1/2*c)^2-A*(sin(1/2*d*x+1/2*c)^2)^(1/2)*(2*sin(1/ 
2*d*x+1/2*c)^2-1)^(1/2)*EllipticE(cos(1/2*d*x+1/2*c),2^(1/2))-B*(sin(1/2*d 
*x+1/2*c)^2)^(1/2)*(2*sin(1/2*d*x+1/2*c)^2-1)^(1/2)*EllipticF(cos(1/2*d*x+ 
1/2*c),2^(1/2)))/(-b*(2*sin(1/2*d*x+1/2*c)^4-sin(1/2*d*x+1/2*c)^2))^(1/2)/ 
sin(1/2*d*x+1/2*c)/((2*cos(1/2*d*x+1/2*c)^2-1)*b)^(1/2)/d
 

Fricas [C] (verification not implemented)

Result contains complex when optimal does not.

Time = 0.09 (sec) , antiderivative size = 172, normalized size of antiderivative = 1.56 \[ \int (b \cos (c+d x))^{3/2} (A+B \cos (c+d x)) \sec ^3(c+d x) \, dx=-\frac {2 \, {\left (i \, \sqrt {\frac {1}{2}} B b^{\frac {3}{2}} \cos \left (d x + c\right ) {\rm weierstrassPInverse}\left (-4, 0, \cos \left (d x + c\right ) + i \, \sin \left (d x + c\right )\right ) - i \, \sqrt {\frac {1}{2}} B b^{\frac {3}{2}} \cos \left (d x + c\right ) {\rm weierstrassPInverse}\left (-4, 0, \cos \left (d x + c\right ) - i \, \sin \left (d x + c\right )\right ) + i \, \sqrt {\frac {1}{2}} A b^{\frac {3}{2}} \cos \left (d x + c\right ) {\rm weierstrassZeta}\left (-4, 0, {\rm weierstrassPInverse}\left (-4, 0, \cos \left (d x + c\right ) + i \, \sin \left (d x + c\right )\right )\right ) - i \, \sqrt {\frac {1}{2}} A b^{\frac {3}{2}} \cos \left (d x + c\right ) {\rm weierstrassZeta}\left (-4, 0, {\rm weierstrassPInverse}\left (-4, 0, \cos \left (d x + c\right ) - i \, \sin \left (d x + c\right )\right )\right ) - \sqrt {b \cos \left (d x + c\right )} A b \sin \left (d x + c\right )\right )}}{d \cos \left (d x + c\right )} \] Input:

integrate((b*cos(d*x+c))^(3/2)*(A+B*cos(d*x+c))*sec(d*x+c)^3,x, algorithm= 
"fricas")
 

Output:

-2*(I*sqrt(1/2)*B*b^(3/2)*cos(d*x + c)*weierstrassPInverse(-4, 0, cos(d*x 
+ c) + I*sin(d*x + c)) - I*sqrt(1/2)*B*b^(3/2)*cos(d*x + c)*weierstrassPIn 
verse(-4, 0, cos(d*x + c) - I*sin(d*x + c)) + I*sqrt(1/2)*A*b^(3/2)*cos(d* 
x + c)*weierstrassZeta(-4, 0, weierstrassPInverse(-4, 0, cos(d*x + c) + I* 
sin(d*x + c))) - I*sqrt(1/2)*A*b^(3/2)*cos(d*x + c)*weierstrassZeta(-4, 0, 
 weierstrassPInverse(-4, 0, cos(d*x + c) - I*sin(d*x + c))) - sqrt(b*cos(d 
*x + c))*A*b*sin(d*x + c))/(d*cos(d*x + c))
 

Sympy [F(-1)]

Timed out. \[ \int (b \cos (c+d x))^{3/2} (A+B \cos (c+d x)) \sec ^3(c+d x) \, dx=\text {Timed out} \] Input:

integrate((b*cos(d*x+c))**(3/2)*(A+B*cos(d*x+c))*sec(d*x+c)**3,x)
 

Output:

Timed out
 

Maxima [F]

\[ \int (b \cos (c+d x))^{3/2} (A+B \cos (c+d x)) \sec ^3(c+d x) \, dx=\int { {\left (B \cos \left (d x + c\right ) + A\right )} \left (b \cos \left (d x + c\right )\right )^{\frac {3}{2}} \sec \left (d x + c\right )^{3} \,d x } \] Input:

integrate((b*cos(d*x+c))^(3/2)*(A+B*cos(d*x+c))*sec(d*x+c)^3,x, algorithm= 
"maxima")
 

Output:

integrate((B*cos(d*x + c) + A)*(b*cos(d*x + c))^(3/2)*sec(d*x + c)^3, x)
 

Giac [F]

\[ \int (b \cos (c+d x))^{3/2} (A+B \cos (c+d x)) \sec ^3(c+d x) \, dx=\int { {\left (B \cos \left (d x + c\right ) + A\right )} \left (b \cos \left (d x + c\right )\right )^{\frac {3}{2}} \sec \left (d x + c\right )^{3} \,d x } \] Input:

integrate((b*cos(d*x+c))^(3/2)*(A+B*cos(d*x+c))*sec(d*x+c)^3,x, algorithm= 
"giac")
 

Output:

integrate((B*cos(d*x + c) + A)*(b*cos(d*x + c))^(3/2)*sec(d*x + c)^3, x)
 

Mupad [F(-1)]

Timed out. \[ \int (b \cos (c+d x))^{3/2} (A+B \cos (c+d x)) \sec ^3(c+d x) \, dx=\int \frac {{\left (b\,\cos \left (c+d\,x\right )\right )}^{3/2}\,\left (A+B\,\cos \left (c+d\,x\right )\right )}{{\cos \left (c+d\,x\right )}^3} \,d x \] Input:

int(((b*cos(c + d*x))^(3/2)*(A + B*cos(c + d*x)))/cos(c + d*x)^3,x)
                                                                                    
                                                                                    
 

Output:

int(((b*cos(c + d*x))^(3/2)*(A + B*cos(c + d*x)))/cos(c + d*x)^3, x)
 

Reduce [F]

\[ \int (b \cos (c+d x))^{3/2} (A+B \cos (c+d x)) \sec ^3(c+d x) \, dx=\sqrt {b}\, b \left (\left (\int \sqrt {\cos \left (d x +c \right )}\, \cos \left (d x +c \right ) \sec \left (d x +c \right )^{3}d x \right ) a +\left (\int \sqrt {\cos \left (d x +c \right )}\, \cos \left (d x +c \right )^{2} \sec \left (d x +c \right )^{3}d x \right ) b \right ) \] Input:

int((b*cos(d*x+c))^(3/2)*(A+B*cos(d*x+c))*sec(d*x+c)^3,x)
 

Output:

sqrt(b)*b*(int(sqrt(cos(c + d*x))*cos(c + d*x)*sec(c + d*x)**3,x)*a + int( 
sqrt(cos(c + d*x))*cos(c + d*x)**2*sec(c + d*x)**3,x)*b)