\(\int x^2 \sin (a+b x) \text {Si}(a+b x) \, dx\) [55]

Optimal result

Integrand size = 16, antiderivative size = 187 \[ \int x^2 \sin (a+b x) \text {Si}(a+b x) \, dx=-\frac {x}{b^2}+\frac {a \cos (2 a+2 b x)}{4 b^3}-\frac {x \cos (2 a+2 b x)}{4 b^2}-\frac {a \operatorname {CosIntegral}(2 a+2 b x)}{b^3}+\frac {a \log (a+b x)}{b^3}+\frac {\cos (a+b x) \sin (a+b x)}{b^3}+\frac {\sin (2 a+2 b x)}{8 b^3}+\frac {2 \cos (a+b x) \text {Si}(a+b x)}{b^3}-\frac {x^2 \cos (a+b x) \text {Si}(a+b x)}{b}+\frac {2 x \sin (a+b x) \text {Si}(a+b x)}{b^2}-\frac {\text {Si}(2 a+2 b x)}{b^3}+\frac {a^2 \text {Si}(2 a+2 b x)}{2 b^3} \]

[Out]

-x/b^2-a*Ci(2*b*x+2*a)/b^3+1/4*a*cos(2*b*x+2*a)/b^3-1/4*x*cos(2*b*x+2*a)/b^2+a*ln(b*x+a)/b^3+2*cos(b*x+a)*Si(b
*x+a)/b^3-x^2*cos(b*x+a)*Si(b*x+a)/b-Si(2*b*x+2*a)/b^3+1/2*a^2*Si(2*b*x+2*a)/b^3+cos(b*x+a)*sin(b*x+a)/b^3+2*x
*Si(b*x+a)*sin(b*x+a)/b^2+1/8*sin(2*b*x+2*a)/b^3

Rubi [A] (verified)

Time = 0.47 (sec) , antiderivative size = 187, normalized size of antiderivative = 1.00, number of steps used = 21, number of rules used = 16, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 1.000, Rules used = {6648, 4669, 6873, 6874, 2718, 3377, 2717, 3380, 6654, 2715, 8, 3393, 3383, 6646, 4491, 12} \[ \int x^2 \sin (a+b x) \text {Si}(a+b x) \, dx=\frac {a^2 \text {Si}(2 a+2 b x)}{2 b^3}-\frac {a \operatorname {CosIntegral}(2 a+2 b x)}{b^3}-\frac {\text {Si}(2 a+2 b x)}{b^3}+\frac {2 \text {Si}(a+b x) \cos (a+b x)}{b^3}+\frac {a \log (a+b x)}{b^3}+\frac {\sin (2 a+2 b x)}{8 b^3}+\frac {a \cos (2 a+2 b x)}{4 b^3}+\frac {\sin (a+b x) \cos (a+b x)}{b^3}+\frac {2 x \text {Si}(a+b x) \sin (a+b x)}{b^2}-\frac {x \cos (2 a+2 b x)}{4 b^2}-\frac {x^2 \text {Si}(a+b x) \cos (a+b x)}{b}-\frac {x}{b^2} \]

[In]

Int[x^2*Sin[a + b*x]*SinIntegral[a + b*x],x]

[Out]

-(x/b^2) + (a*Cos[2*a + 2*b*x])/(4*b^3) - (x*Cos[2*a + 2*b*x])/(4*b^2) - (a*CosIntegral[2*a + 2*b*x])/b^3 + (a
*Log[a + b*x])/b^3 + (Cos[a + b*x]*Sin[a + b*x])/b^3 + Sin[2*a + 2*b*x]/(8*b^3) + (2*Cos[a + b*x]*SinIntegral[
a + b*x])/b^3 - (x^2*Cos[a + b*x]*SinIntegral[a + b*x])/b + (2*x*Sin[a + b*x]*SinIntegral[a + b*x])/b^2 - SinI
ntegral[2*a + 2*b*x]/b^3 + (a^2*SinIntegral[2*a + 2*b*x])/(2*b^3)

Rule 8

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

Rule 12

Int[(a_)*(u_), x_Symbol] :> Dist[a, Int[u, x], x] /; FreeQ[a, x] &&  !MatchQ[u, (b_)*(v_) /; FreeQ[b, x]]

Rule 2715

Int[((b_.)*sin[(c_.) + (d_.)*(x_)])^(n_), x_Symbol] :> Simp[(-b)*Cos[c + d*x]*((b*Sin[c + d*x])^(n - 1)/(d*n))
, x] + Dist[b^2*((n - 1)/n), Int[(b*Sin[c + d*x])^(n - 2), x], x] /; FreeQ[{b, c, d}, x] && GtQ[n, 1] && Integ
erQ[2*n]

Rule 2717

Int[sin[Pi/2 + (c_.) + (d_.)*(x_)], x_Symbol] :> Simp[Sin[c + d*x]/d, x] /; FreeQ[{c, d}, x]

Rule 2718

Int[sin[(c_.) + (d_.)*(x_)], x_Symbol] :> Simp[-Cos[c + d*x]/d, x] /; FreeQ[{c, d}, x]

Rule 3377

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

Rule 3380

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]

Rule 3383

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

Rule 3393

Int[((c_.) + (d_.)*(x_))^(m_)*sin[(e_.) + (f_.)*(x_)]^(n_), x_Symbol] :> Int[ExpandTrigReduce[(c + d*x)^m, Sin
[e + f*x]^n, x], x] /; FreeQ[{c, d, e, f, m}, x] && IGtQ[n, 1] && ( !RationalQ[m] || (GeQ[m, -1] && LtQ[m, 1])
)

Rule 4491

Int[Cos[(a_.) + (b_.)*(x_)]^(p_.)*((c_.) + (d_.)*(x_))^(m_.)*Sin[(a_.) + (b_.)*(x_)]^(n_.), x_Symbol] :> Int[E
xpandTrigReduce[(c + d*x)^m, Sin[a + b*x]^n*Cos[a + b*x]^p, x], x] /; FreeQ[{a, b, c, d, m}, x] && IGtQ[n, 0]
&& IGtQ[p, 0]

Rule 4669

Int[Cos[w_]^(p_.)*(u_.)*Sin[v_]^(p_.), x_Symbol] :> Dist[1/2^p, Int[u*Sin[2*v]^p, x], x] /; EqQ[w, v] && Integ
erQ[p]

Rule 6646

Int[Sin[(a_.) + (b_.)*(x_)]*SinIntegral[(c_.) + (d_.)*(x_)], x_Symbol] :> Simp[(-Cos[a + b*x])*(SinIntegral[c
+ d*x]/b), x] + Dist[d/b, Int[Cos[a + b*x]*(Sin[c + d*x]/(c + d*x)), x], x] /; FreeQ[{a, b, c, d}, x]

Rule 6648

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

Rule 6654

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

Rule 6873

Int[u_, x_Symbol] :> With[{v = NormalizeIntegrand[u, x]}, Int[v, x] /; v =!= u]

Rule 6874

Int[u_, x_Symbol] :> With[{v = ExpandIntegrand[u, x]}, Int[v, x] /; SumQ[v]]

Rubi steps \begin{align*} \text {integral}& = -\frac {x^2 \cos (a+b x) \text {Si}(a+b x)}{b}+\frac {2 \int x \cos (a+b x) \text {Si}(a+b x) \, dx}{b}+\int \frac {x^2 \cos (a+b x) \sin (a+b x)}{a+b x} \, dx \\ & = -\frac {x^2 \cos (a+b x) \text {Si}(a+b x)}{b}+\frac {2 x \sin (a+b x) \text {Si}(a+b x)}{b^2}+\frac {1}{2} \int \frac {x^2 \sin (2 (a+b x))}{a+b x} \, dx-\frac {2 \int \sin (a+b x) \text {Si}(a+b x) \, dx}{b^2}-\frac {2 \int \frac {x \sin ^2(a+b x)}{a+b x} \, dx}{b} \\ & = \frac {2 \cos (a+b x) \text {Si}(a+b x)}{b^3}-\frac {x^2 \cos (a+b x) \text {Si}(a+b x)}{b}+\frac {2 x \sin (a+b x) \text {Si}(a+b x)}{b^2}+\frac {1}{2} \int \frac {x^2 \sin (2 a+2 b x)}{a+b x} \, dx-\frac {2 \int \frac {\cos (a+b x) \sin (a+b x)}{a+b x} \, dx}{b^2}-\frac {2 \int \left (\frac {\sin ^2(a+b x)}{b}-\frac {a \sin ^2(a+b x)}{b (a+b x)}\right ) \, dx}{b} \\ & = \frac {2 \cos (a+b x) \text {Si}(a+b x)}{b^3}-\frac {x^2 \cos (a+b x) \text {Si}(a+b x)}{b}+\frac {2 x \sin (a+b x) \text {Si}(a+b x)}{b^2}+\frac {1}{2} \int \left (-\frac {a \sin (2 a+2 b x)}{b^2}+\frac {x \sin (2 a+2 b x)}{b}+\frac {a^2 \sin (2 a+2 b x)}{b^2 (a+b x)}\right ) \, dx-\frac {2 \int \sin ^2(a+b x) \, dx}{b^2}-\frac {2 \int \frac {\sin (2 a+2 b x)}{2 (a+b x)} \, dx}{b^2}+\frac {(2 a) \int \frac {\sin ^2(a+b x)}{a+b x} \, dx}{b^2} \\ & = \frac {\cos (a+b x) \sin (a+b x)}{b^3}+\frac {2 \cos (a+b x) \text {Si}(a+b x)}{b^3}-\frac {x^2 \cos (a+b x) \text {Si}(a+b x)}{b}+\frac {2 x \sin (a+b x) \text {Si}(a+b x)}{b^2}-\frac {\int 1 \, dx}{b^2}-\frac {\int \frac {\sin (2 a+2 b x)}{a+b x} \, dx}{b^2}-\frac {a \int \sin (2 a+2 b x) \, dx}{2 b^2}+\frac {(2 a) \int \left (\frac {1}{2 (a+b x)}-\frac {\cos (2 a+2 b x)}{2 (a+b x)}\right ) \, dx}{b^2}+\frac {a^2 \int \frac {\sin (2 a+2 b x)}{a+b x} \, dx}{2 b^2}+\frac {\int x \sin (2 a+2 b x) \, dx}{2 b} \\ & = -\frac {x}{b^2}+\frac {a \cos (2 a+2 b x)}{4 b^3}-\frac {x \cos (2 a+2 b x)}{4 b^2}+\frac {a \log (a+b x)}{b^3}+\frac {\cos (a+b x) \sin (a+b x)}{b^3}+\frac {2 \cos (a+b x) \text {Si}(a+b x)}{b^3}-\frac {x^2 \cos (a+b x) \text {Si}(a+b x)}{b}+\frac {2 x \sin (a+b x) \text {Si}(a+b x)}{b^2}-\frac {\text {Si}(2 a+2 b x)}{b^3}+\frac {a^2 \text {Si}(2 a+2 b x)}{2 b^3}+\frac {\int \cos (2 a+2 b x) \, dx}{4 b^2}-\frac {a \int \frac {\cos (2 a+2 b x)}{a+b x} \, dx}{b^2} \\ & = -\frac {x}{b^2}+\frac {a \cos (2 a+2 b x)}{4 b^3}-\frac {x \cos (2 a+2 b x)}{4 b^2}-\frac {a \operatorname {CosIntegral}(2 a+2 b x)}{b^3}+\frac {a \log (a+b x)}{b^3}+\frac {\cos (a+b x) \sin (a+b x)}{b^3}+\frac {\sin (2 a+2 b x)}{8 b^3}+\frac {2 \cos (a+b x) \text {Si}(a+b x)}{b^3}-\frac {x^2 \cos (a+b x) \text {Si}(a+b x)}{b}+\frac {2 x \sin (a+b x) \text {Si}(a+b x)}{b^2}-\frac {\text {Si}(2 a+2 b x)}{b^3}+\frac {a^2 \text {Si}(2 a+2 b x)}{2 b^3} \\ \end{align*}

Mathematica [A] (verified)

Time = 0.25 (sec) , antiderivative size = 123, normalized size of antiderivative = 0.66 \[ \int x^2 \sin (a+b x) \text {Si}(a+b x) \, dx=\frac {-8 b x+2 a \cos (2 (a+b x))-2 b x \cos (2 (a+b x))-8 a \operatorname {CosIntegral}(2 (a+b x))+8 a \log (a+b x)+5 \sin (2 (a+b x))-8 \left (\left (-2+b^2 x^2\right ) \cos (a+b x)-2 b x \sin (a+b x)\right ) \text {Si}(a+b x)-8 \text {Si}(2 (a+b x))+4 a^2 \text {Si}(2 (a+b x))}{8 b^3} \]

[In]

Integrate[x^2*Sin[a + b*x]*SinIntegral[a + b*x],x]

[Out]

(-8*b*x + 2*a*Cos[2*(a + b*x)] - 2*b*x*Cos[2*(a + b*x)] - 8*a*CosIntegral[2*(a + b*x)] + 8*a*Log[a + b*x] + 5*
Sin[2*(a + b*x)] - 8*((-2 + b^2*x^2)*Cos[a + b*x] - 2*b*x*Sin[a + b*x])*SinIntegral[a + b*x] - 8*SinIntegral[2
*(a + b*x)] + 4*a^2*SinIntegral[2*(a + b*x)])/(8*b^3)

Maple [A] (verified)

Time = 2.51 (sec) , antiderivative size = 175, normalized size of antiderivative = 0.94

method	result	size
derivativedivides	\(\frac {\operatorname {Si}\left (b x +a \right ) \left (-a^{2} \cos \left (b x +a \right )-2 a \left (\sin \left (b x +a \right )-\left (b x +a \right ) \cos \left (b x +a \right )\right )-\left (b x +a \right )^{2} \cos \left (b x +a \right )+2 \cos \left (b x +a \right )+2 \left (b x +a \right ) \sin \left (b x +a \right )\right )+\frac {a^{2} \operatorname {Si}\left (2 b x +2 a \right )}{2}+a \cos \left (b x +a \right )^{2}-\frac {\left (b x +a \right ) \cos \left (b x +a \right )^{2}}{2}+\frac {5 \sin \left (b x +a \right ) \cos \left (b x +a \right )}{4}-\frac {3 b x}{4}-\frac {3 a}{4}+a \ln \left (b x +a \right )-a \,\operatorname {Ci}\left (2 b x +2 a \right )-\operatorname {Si}\left (2 b x +2 a \right )}{b^{3}}\)	\(175\)
default	\(\frac {\operatorname {Si}\left (b x +a \right ) \left (-a^{2} \cos \left (b x +a \right )-2 a \left (\sin \left (b x +a \right )-\left (b x +a \right ) \cos \left (b x +a \right )\right )-\left (b x +a \right )^{2} \cos \left (b x +a \right )+2 \cos \left (b x +a \right )+2 \left (b x +a \right ) \sin \left (b x +a \right )\right )+\frac {a^{2} \operatorname {Si}\left (2 b x +2 a \right )}{2}+a \cos \left (b x +a \right )^{2}-\frac {\left (b x +a \right ) \cos \left (b x +a \right )^{2}}{2}+\frac {5 \sin \left (b x +a \right ) \cos \left (b x +a \right )}{4}-\frac {3 b x}{4}-\frac {3 a}{4}+a \ln \left (b x +a \right )-a \,\operatorname {Ci}\left (2 b x +2 a \right )-\operatorname {Si}\left (2 b x +2 a \right )}{b^{3}}\)	\(175\)

[In]

int(x^2*Si(b*x+a)*sin(b*x+a),x,method=_RETURNVERBOSE)

[Out]

1/b^3*(Si(b*x+a)*(-a^2*cos(b*x+a)-2*a*(sin(b*x+a)-(b*x+a)*cos(b*x+a))-(b*x+a)^2*cos(b*x+a)+2*cos(b*x+a)+2*(b*x
+a)*sin(b*x+a))+1/2*a^2*Si(2*b*x+2*a)+a*cos(b*x+a)^2-1/2*(b*x+a)*cos(b*x+a)^2+5/4*sin(b*x+a)*cos(b*x+a)-3/4*b*
x-3/4*a+a*ln(b*x+a)-a*Ci(2*b*x+2*a)-Si(2*b*x+2*a))

Fricas [A] (verification not implemented)

none

Time = 0.26 (sec) , antiderivative size = 114, normalized size of antiderivative = 0.61 \[ \int x^2 \sin (a+b x) \text {Si}(a+b x) \, dx=-\frac {2 \, {\left (b x - a\right )} \cos \left (b x + a\right )^{2} + 4 \, {\left (b^{2} x^{2} - 2\right )} \cos \left (b x + a\right ) \operatorname {Si}\left (b x + a\right ) + 3 \, b x + 4 \, a \operatorname {Ci}\left (2 \, b x + 2 \, a\right ) - 4 \, a \log \left (b x + a\right ) - {\left (8 \, b x \operatorname {Si}\left (b x + a\right ) + 5 \, \cos \left (b x + a\right )\right )} \sin \left (b x + a\right ) - 2 \, {\left (a^{2} - 2\right )} \operatorname {Si}\left (2 \, b x + 2 \, a\right )}{4 \, b^{3}} \]

[In]

integrate(x^2*sin_integral(b*x+a)*sin(b*x+a),x, algorithm="fricas")

[Out]

-1/4*(2*(b*x - a)*cos(b*x + a)^2 + 4*(b^2*x^2 - 2)*cos(b*x + a)*sin_integral(b*x + a) + 3*b*x + 4*a*cos_integr
al(2*b*x + 2*a) - 4*a*log(b*x + a) - (8*b*x*sin_integral(b*x + a) + 5*cos(b*x + a))*sin(b*x + a) - 2*(a^2 - 2)
*sin_integral(2*b*x + 2*a))/b^3

Sympy [F]

\[ \int x^2 \sin (a+b x) \text {Si}(a+b x) \, dx=\int x^{2} \sin {\left (a + b x \right )} \operatorname {Si}{\left (a + b x \right )}\, dx \]

[In]

integrate(x**2*Si(b*x+a)*sin(b*x+a),x)

[Out]

Integral(x**2*sin(a + b*x)*Si(a + b*x), x)

Maxima [F]

\[ \int x^2 \sin (a+b x) \text {Si}(a+b x) \, dx=\int { x^{2} \sin \left (b x + a\right ) \operatorname {Si}\left (b x + a\right ) \,d x } \]

[In]

integrate(x^2*sin_integral(b*x+a)*sin(b*x+a),x, algorithm="maxima")

[Out]

integrate(x^2*sin(b*x + a)*sin_integral(b*x + a), x)

Giac [C] (verification not implemented)

Result contains higher order function than in optimal. Order 9 vs. order 4.

Time = 0.33 (sec) , antiderivative size = 398, normalized size of antiderivative = 2.13 \[ \int x^2 \sin (a+b x) \text {Si}(a+b x) \, dx={\left (\frac {2 \, x \sin \left (b x + a\right )}{b^{2}} - \frac {{\left (b^{2} x^{2} - 2\right )} \cos \left (b x + a\right )}{b^{3}}\right )} \operatorname {Si}\left (b x + a\right ) + \frac {a^{2} \Im \left ( \operatorname {Ci}\left (2 \, b x + 2 \, a\right ) \right ) \tan \left (b x + a\right )^{2} - a^{2} \Im \left ( \operatorname {Ci}\left (-2 \, b x - 2 \, a\right ) \right ) \tan \left (b x + a\right )^{2} + 2 \, a^{2} \operatorname {Si}\left (2 \, b x + 2 \, a\right ) \tan \left (b x + a\right )^{2} - 3 \, b x \tan \left (b x + a\right )^{2} + 4 \, a \log \left ({\left | b x + a \right |}\right ) \tan \left (b x + a\right )^{2} - 2 \, a \Re \left ( \operatorname {Ci}\left (2 \, b x + 2 \, a\right ) \right ) \tan \left (b x + a\right )^{2} - 2 \, a \Re \left ( \operatorname {Ci}\left (-2 \, b x - 2 \, a\right ) \right ) \tan \left (b x + a\right )^{2} + a^{2} \Im \left ( \operatorname {Ci}\left (2 \, b x + 2 \, a\right ) \right ) - a^{2} \Im \left ( \operatorname {Ci}\left (-2 \, b x - 2 \, a\right ) \right ) + 2 \, a^{2} \operatorname {Si}\left (2 \, b x + 2 \, a\right ) - a \tan \left (b x + a\right )^{2} - 2 \, \Im \left ( \operatorname {Ci}\left (2 \, b x + 2 \, a\right ) \right ) \tan \left (b x + a\right )^{2} + 2 \, \Im \left ( \operatorname {Ci}\left (-2 \, b x - 2 \, a\right ) \right ) \tan \left (b x + a\right )^{2} - 4 \, \operatorname {Si}\left (2 \, b x + 2 \, a\right ) \tan \left (b x + a\right )^{2} - 5 \, b x + 4 \, a \log \left ({\left | b x + a \right |}\right ) - 2 \, a \Re \left ( \operatorname {Ci}\left (2 \, b x + 2 \, a\right ) \right ) - 2 \, a \Re \left ( \operatorname {Ci}\left (-2 \, b x - 2 \, a\right ) \right ) + a - 2 \, \Im \left ( \operatorname {Ci}\left (2 \, b x + 2 \, a\right ) \right ) + 2 \, \Im \left ( \operatorname {Ci}\left (-2 \, b x - 2 \, a\right ) \right ) - 4 \, \operatorname {Si}\left (2 \, b x + 2 \, a\right ) + 5 \, \tan \left (b x + a\right )}{4 \, {\left (b^{3} \tan \left (b x + a\right )^{2} + b^{3}\right )}} \]

[In]

integrate(x^2*sin_integral(b*x+a)*sin(b*x+a),x, algorithm="giac")

[Out]

(2*x*sin(b*x + a)/b^2 - (b^2*x^2 - 2)*cos(b*x + a)/b^3)*sin_integral(b*x + a) + 1/4*(a^2*imag_part(cos_integra
l(2*b*x + 2*a))*tan(b*x + a)^2 - a^2*imag_part(cos_integral(-2*b*x - 2*a))*tan(b*x + a)^2 + 2*a^2*sin_integral
(2*b*x + 2*a)*tan(b*x + a)^2 - 3*b*x*tan(b*x + a)^2 + 4*a*log(abs(b*x + a))*tan(b*x + a)^2 - 2*a*real_part(cos
_integral(2*b*x + 2*a))*tan(b*x + a)^2 - 2*a*real_part(cos_integral(-2*b*x - 2*a))*tan(b*x + a)^2 + a^2*imag_p
art(cos_integral(2*b*x + 2*a)) - a^2*imag_part(cos_integral(-2*b*x - 2*a)) + 2*a^2*sin_integral(2*b*x + 2*a) -
a*tan(b*x + a)^2 - 2*imag_part(cos_integral(2*b*x + 2*a))*tan(b*x + a)^2 + 2*imag_part(cos_integral(-2*b*x -
2*a))*tan(b*x + a)^2 - 4*sin_integral(2*b*x + 2*a)*tan(b*x + a)^2 - 5*b*x + 4*a*log(abs(b*x + a)) - 2*a*real_p
art(cos_integral(2*b*x + 2*a)) - 2*a*real_part(cos_integral(-2*b*x - 2*a)) + a - 2*imag_part(cos_integral(2*b*
x + 2*a)) + 2*imag_part(cos_integral(-2*b*x - 2*a)) - 4*sin_integral(2*b*x + 2*a) + 5*tan(b*x + a))/(b^3*tan(b
*x + a)^2 + b^3)

Mupad [F(-1)]

Timed out. \[ \int x^2 \sin (a+b x) \text {Si}(a+b x) \, dx=\int x^2\,\mathrm {sinint}\left (a+b\,x\right )\,\sin \left (a+b\,x\right ) \,d x \]

[In]

int(x^2*sinint(a + b*x)*sin(a + b*x),x)

[Out]

int(x^2*sinint(a + b*x)*sin(a + b*x), x)