3.588 \(\int \frac{\sin ^3(a x)}{x (a x \cos (a x)-\sin (a x))^2} \, dx\)

Optimal. Leaf size=56 \[ \frac{\sin (a x)}{a^2 x^2}+\frac{\sin ^2(a x)}{a^2 x^2 (a x \cos (a x)-\sin (a x))}+\text{Si}(a x)+\frac{\cos (a x)}{a x} \]

[Out]

Cos[a*x]/(a*x) + Sin[a*x]/(a^2*x^2) + Sin[a*x]^2/(a^2*x^2*(a*x*Cos[a*x] - Sin[a*x])) + SinIntegral[a*x]

________________________________________________________________________________________

Rubi [A]  time = 0.101601, antiderivative size = 56, normalized size of antiderivative = 1., number of steps used = 4, number of rules used = 3, integrand size = 26, \(\frac{\text{number of rules}}{\text{integrand size}}\) = 0.115, Rules used = {4598, 3297, 3299} \[ \frac{\sin (a x)}{a^2 x^2}+\frac{\sin ^2(a x)}{a^2 x^2 (a x \cos (a x)-\sin (a x))}+\text{Si}(a x)+\frac{\cos (a x)}{a x} \]

Antiderivative was successfully verified.

[In]

Int[Sin[a*x]^3/(x*(a*x*Cos[a*x] - Sin[a*x])^2),x]

[Out]

Cos[a*x]/(a*x) + Sin[a*x]/(a^2*x^2) + Sin[a*x]^2/(a^2*x^2*(a*x*Cos[a*x] - Sin[a*x])) + SinIntegral[a*x]

Rule 4598

Int[(((b_.)*(x_))^(m_)*Sin[(a_.)*(x_)]^(n_))/(Cos[(a_.)*(x_)]*(d_.)*(x_) + (c_.)*Sin[(a_.)*(x_)])^2, x_Symbol]
:> Simp[(b*(b*x)^(m - 1)*Sin[a*x]^(n - 1))/(a*d*(c*Sin[a*x] + d*x*Cos[a*x])), x] - Dist[(b^2*(n - 1))/d^2, In
t[(b*x)^(m - 2)*Sin[a*x]^(n - 2), x], x] /; FreeQ[{a, b, c, d, m, n}, x] && EqQ[a*c + d, 0] && EqQ[m, 2 - n]

Rule 3297

Int[((c_.) + (d_.)*(x_))^(m_)*sin[(e_.) + (f_.)*(x_)], x_Symbol] :> Simp[((c + d*x)^(m + 1)*Sin[e + f*x])/(d*(
m + 1)), x] - Dist[f/(d*(m + 1)), Int[(c + d*x)^(m + 1)*Cos[e + f*x], x], x] /; FreeQ[{c, d, e, f}, x] && LtQ[
m, -1]

Rule 3299

Int[sin[(e_.) + (f_.)*(x_)]/((c_.) + (d_.)*(x_)), x_Symbol] :> Simp[SinIntegral[e + f*x]/d, x] /; FreeQ[{c, d,
e, f}, x] && EqQ[d*e - c*f, 0]

Rubi steps

\begin{align*} \int \frac{\sin ^3(a x)}{x (a x \cos (a x)-\sin (a x))^2} \, dx &=\frac{\sin ^2(a x)}{a^2 x^2 (a x \cos (a x)-\sin (a x))}-\frac{2 \int \frac{\sin (a x)}{x^3} \, dx}{a^2}\\ &=\frac{\sin (a x)}{a^2 x^2}+\frac{\sin ^2(a x)}{a^2 x^2 (a x \cos (a x)-\sin (a x))}-\frac{\int \frac{\cos (a x)}{x^2} \, dx}{a}\\ &=\frac{\cos (a x)}{a x}+\frac{\sin (a x)}{a^2 x^2}+\frac{\sin ^2(a x)}{a^2 x^2 (a x \cos (a x)-\sin (a x))}+\int \frac{\sin (a x)}{x} \, dx\\ &=\frac{\cos (a x)}{a x}+\frac{\sin (a x)}{a^2 x^2}+\frac{\sin ^2(a x)}{a^2 x^2 (a x \cos (a x)-\sin (a x))}+\text{Si}(a x)\\ \end{align*}

Mathematica [C]  time = 7.54569, size = 242, normalized size = 4.32 \[ \frac{-i e \text{CosIntegral}(-a x+i) (a x \cos (a x)-\sin (a x))+i e \text{CosIntegral}(a x+i) (a x \cos (a x)-\sin (a x))-i e \text{ExpIntegralEi}(-1-i a x) \sin (a x)+i e \text{ExpIntegralEi}(-1+i a x) \sin (a x)+i e a x \text{ExpIntegralEi}(-1-i a x) \cos (a x)-i e a x \text{ExpIntegralEi}(-1+i a x) \cos (a x)-2 \text{Si}(a x) \sin (a x)-e \text{Si}(i-a x) \sin (a x)+e \text{Si}(a x+i) \sin (a x)+2 a x \text{Si}(a x) \cos (a x)+e a x \text{Si}(i-a x) \cos (a x)-e a x \text{Si}(a x+i) \cos (a x)+\cos (2 a x)+1}{2 a x \cos (a x)-2 \sin (a x)} \]

Warning: Unable to verify antiderivative.

[In]

Integrate[Sin[a*x]^3/(x*(a*x*Cos[a*x] - Sin[a*x])^2),x]

[Out]

(1 + Cos[2*a*x] + I*a*E*x*Cos[a*x]*ExpIntegralEi[-1 - I*a*x] - I*a*E*x*Cos[a*x]*ExpIntegralEi[-1 + I*a*x] - I*
E*CosIntegral[I - a*x]*(a*x*Cos[a*x] - Sin[a*x]) + I*E*CosIntegral[I + a*x]*(a*x*Cos[a*x] - Sin[a*x]) - I*E*Ex
pIntegralEi[-1 - I*a*x]*Sin[a*x] + I*E*ExpIntegralEi[-1 + I*a*x]*Sin[a*x] + 2*a*x*Cos[a*x]*SinIntegral[a*x] -
2*Sin[a*x]*SinIntegral[a*x] + a*E*x*Cos[a*x]*SinIntegral[I - a*x] - E*Sin[a*x]*SinIntegral[I - a*x] - a*E*x*Co
s[a*x]*SinIntegral[I + a*x] + E*Sin[a*x]*SinIntegral[I + a*x])/(2*a*x*Cos[a*x] - 2*Sin[a*x])

________________________________________________________________________________________

Maple [C]  time = 0.666, size = 108, normalized size = 1.9 \begin{align*}{\frac{{\frac{i}{2}}{{\rm e}^{iax}}}{-1+iax}}+{\frac{i}{2}}{\it Ei} \left ( 1,-iax \right ) +{\frac{{\frac{i}{2}}{{\rm e}^{-iax}}}{1+iax}}-{\frac{i}{2}}{\it Ei} \left ( 1,iax \right ) +2\,{\frac{{{\rm e}^{iax}}}{ \left ( ax+i \right ) \left ( ax-i \right ) \left ( ax{{\rm e}^{2\,iax}}+i{{\rm e}^{2\,iax}}+ax-i \right ) }} \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

int(sin(a*x)^3/x/(a*x*cos(a*x)-sin(a*x))^2,x)

[Out]

1/2*I*exp(I*a*x)/(-1+I*a*x)+1/2*I*Ei(1,-I*a*x)+1/2*I*exp(-I*a*x)/(1+I*a*x)-1/2*I*Ei(1,I*a*x)+2*exp(I*a*x)/(a*x
+I)/(a*x-I)/(a*x*exp(2*I*a*x)+I*exp(2*I*a*x)+a*x-I)

________________________________________________________________________________________

Maxima [F(-2)]  time = 0., size = 0, normalized size = 0. \begin{align*} \text{Exception raised: RuntimeError} \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate(sin(a*x)^3/x/(a*x*cos(a*x)-sin(a*x))^2,x, algorithm="maxima")

[Out]

Exception raised: RuntimeError

________________________________________________________________________________________

Fricas [A]  time = 2.08375, size = 138, normalized size = 2.46 \begin{align*} \frac{a x \cos \left (a x\right ) \operatorname{Si}\left (a x\right ) + \cos \left (a x\right )^{2} - \sin \left (a x\right ) \operatorname{Si}\left (a x\right )}{a x \cos \left (a x\right ) - \sin \left (a x\right )} \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate(sin(a*x)^3/x/(a*x*cos(a*x)-sin(a*x))^2,x, algorithm="fricas")

[Out]

(a*x*cos(a*x)*sin_integral(a*x) + cos(a*x)^2 - sin(a*x)*sin_integral(a*x))/(a*x*cos(a*x) - sin(a*x))

________________________________________________________________________________________

Sympy [F(-1)]  time = 0., size = 0, normalized size = 0. \begin{align*} \text{Timed out} \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate(sin(a*x)**3/x/(a*x*cos(a*x)-sin(a*x))**2,x)

[Out]

Timed out

________________________________________________________________________________________

Giac [C]  time = 1.28642, size = 670, normalized size = 11.96 \begin{align*} \text{result too large to display} \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate(sin(a*x)^3/x/(a*x*cos(a*x)-sin(a*x))^2,x, algorithm="giac")

[Out]

1/2*(a^3*x^3*imag_part(cos_integral(a*x))*tan(1/2*a*x)^4 - a^3*x^3*imag_part(cos_integral(-a*x))*tan(1/2*a*x)^
4 + 2*a^3*x^3*sin_integral(a*x)*tan(1/2*a*x)^4 + 2*a^2*x^2*imag_part(cos_integral(a*x))*tan(1/2*a*x)^3 - 2*a^2
*x^2*imag_part(cos_integral(-a*x))*tan(1/2*a*x)^3 + 4*a^2*x^2*sin_integral(a*x)*tan(1/2*a*x)^3 - 2*a^2*x^2*tan
(1/2*a*x)^4 - a^3*x^3*imag_part(cos_integral(a*x)) + a^3*x^3*imag_part(cos_integral(-a*x)) - 2*a^3*x^3*sin_int
egral(a*x) + a*x*imag_part(cos_integral(a*x))*tan(1/2*a*x)^4 - a*x*imag_part(cos_integral(-a*x))*tan(1/2*a*x)^
4 + 2*a*x*sin_integral(a*x)*tan(1/2*a*x)^4 + 2*a^2*x^2*imag_part(cos_integral(a*x))*tan(1/2*a*x) - 2*a^2*x^2*i
mag_part(cos_integral(-a*x))*tan(1/2*a*x) + 4*a^2*x^2*sin_integral(a*x)*tan(1/2*a*x) + 4*a^2*x^2*tan(1/2*a*x)^
2 - 2*a^2*x^2 + 2*imag_part(cos_integral(a*x))*tan(1/2*a*x)^3 - 2*imag_part(cos_integral(-a*x))*tan(1/2*a*x)^3
+ 4*sin_integral(a*x)*tan(1/2*a*x)^3 - 4*tan(1/2*a*x)^4 - a*x*imag_part(cos_integral(a*x)) + a*x*imag_part(co
s_integral(-a*x)) - 2*a*x*sin_integral(a*x) + 2*imag_part(cos_integral(a*x))*tan(1/2*a*x) - 2*imag_part(cos_in
tegral(-a*x))*tan(1/2*a*x) + 4*sin_integral(a*x)*tan(1/2*a*x) - 4)/(a^3*x^3*tan(1/2*a*x)^4 + 2*a^2*x^2*tan(1/2
*a*x)^3 - a^3*x^3 + a*x*tan(1/2*a*x)^4 + 2*a^2*x^2*tan(1/2*a*x) + 2*tan(1/2*a*x)^3 - a*x + 2*tan(1/2*a*x))