3.1.43 \(\int \frac {\operatorname {FresnelS}(b x)^2}{x^5} \, dx\) [43]

3.1.43.1 Optimal result

Integrand size = 10, antiderivative size = 127 \[ \int \frac {\operatorname {FresnelS}(b x)^2}{x^5} \, dx=-\frac {b^2}{24 x^2}+\frac {b^2 \cos \left (b^2 \pi x^2\right )}{24 x^2}-\frac {b^3 \pi \cos \left (\frac {1}{2} b^2 \pi x^2\right ) \operatorname {FresnelS}(b x)}{6 x}-\frac {1}{12} b^4 \pi ^2 \operatorname {FresnelS}(b x)^2-\frac {\operatorname {FresnelS}(b x)^2}{4 x^4}-\frac {b \operatorname {FresnelS}(b x) \sin \left (\frac {1}{2} b^2 \pi x^2\right )}{6 x^3}+\frac {1}{12} b^4 \pi \text {Si}\left (b^2 \pi x^2\right ) \]

output

-1/24*b^2/x^2+1/24*b^2*cos(b^2*Pi*x^2)/x^2-1/6*b^3*Pi*cos(1/2*b^2*Pi*x^2)* 
FresnelS(b*x)/x-1/12*b^4*Pi^2*FresnelS(b*x)^2-1/4*FresnelS(b*x)^2/x^4+1/12 
*b^4*Pi*Si(b^2*Pi*x^2)-1/6*b*FresnelS(b*x)*sin(1/2*b^2*Pi*x^2)/x^3
 

3.1.43.2 Mathematica [A] (verified)

Time = 0.01 (sec) , antiderivative size = 127, normalized size of antiderivative = 1.00 \[ \int \frac {\operatorname {FresnelS}(b x)^2}{x^5} \, dx=-\frac {b^2}{24 x^2}+\frac {b^2 \cos \left (b^2 \pi x^2\right )}{24 x^2}-\frac {b^3 \pi \cos \left (\frac {1}{2} b^2 \pi x^2\right ) \operatorname {FresnelS}(b x)}{6 x}-\frac {1}{12} b^4 \pi ^2 \operatorname {FresnelS}(b x)^2-\frac {\operatorname {FresnelS}(b x)^2}{4 x^4}-\frac {b \operatorname {FresnelS}(b x) \sin \left (\frac {1}{2} b^2 \pi x^2\right )}{6 x^3}+\frac {1}{12} b^4 \pi \text {Si}\left (b^2 \pi x^2\right ) \]

input

Integrate[FresnelS[b*x]^2/x^5,x]
 

output

-1/24*b^2/x^2 + (b^2*Cos[b^2*Pi*x^2])/(24*x^2) - (b^3*Pi*Cos[(b^2*Pi*x^2)/ 
2]*FresnelS[b*x])/(6*x) - (b^4*Pi^2*FresnelS[b*x]^2)/12 - FresnelS[b*x]^2/ 
(4*x^4) - (b*FresnelS[b*x]*Sin[(b^2*Pi*x^2)/2])/(6*x^3) + (b^4*Pi*SinInteg 
ral[b^2*Pi*x^2])/12
 

3.1.43.3 Rubi [A] (verified)

Time = 0.88 (sec) , antiderivative size = 142, normalized size of antiderivative = 1.12, number of steps used = 13, number of rules used = 12, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 1.200, Rules used = {6984, 7010, 3861, 3042, 3778, 25, 3042, 3780, 7018, 3856, 6994, 15}

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 \frac {\operatorname {FresnelS}(b x)^2}{x^5} \, dx\)

\(\Big \downarrow \) 6984

\(\displaystyle \frac {1}{2} b \int \frac {\operatorname {FresnelS}(b x) \sin \left (\frac {1}{2} b^2 \pi x^2\right )}{x^4}dx-\frac {\operatorname {FresnelS}(b x)^2}{4 x^4}\)

\(\Big \downarrow \) 7010

\(\displaystyle \frac {1}{2} b \left (\frac {1}{3} \pi b^2 \int \frac {\cos \left (\frac {1}{2} b^2 \pi x^2\right ) \operatorname {FresnelS}(b x)}{x^2}dx-\frac {1}{6} b \int \frac {\cos \left (b^2 \pi x^2\right )}{x^3}dx-\frac {\operatorname {FresnelS}(b x) \sin \left (\frac {1}{2} \pi b^2 x^2\right )}{3 x^3}-\frac {b}{12 x^2}\right )-\frac {\operatorname {FresnelS}(b x)^2}{4 x^4}\)

\(\Big \downarrow \) 3861

\(\displaystyle \frac {1}{2} b \left (\frac {1}{3} \pi b^2 \int \frac {\cos \left (\frac {1}{2} b^2 \pi x^2\right ) \operatorname {FresnelS}(b x)}{x^2}dx-\frac {1}{12} b \int \frac {\cos \left (b^2 \pi x^2\right )}{x^4}dx^2-\frac {\operatorname {FresnelS}(b x) \sin \left (\frac {1}{2} \pi b^2 x^2\right )}{3 x^3}-\frac {b}{12 x^2}\right )-\frac {\operatorname {FresnelS}(b x)^2}{4 x^4}\)

\(\Big \downarrow \) 3042

\(\displaystyle \frac {1}{2} b \left (\frac {1}{3} \pi b^2 \int \frac {\cos \left (\frac {1}{2} b^2 \pi x^2\right ) \operatorname {FresnelS}(b x)}{x^2}dx-\frac {1}{12} b \int \frac {\sin \left (b^2 \pi x^2+\frac {\pi }{2}\right )}{x^4}dx^2-\frac {\operatorname {FresnelS}(b x) \sin \left (\frac {1}{2} \pi b^2 x^2\right )}{3 x^3}-\frac {b}{12 x^2}\right )-\frac {\operatorname {FresnelS}(b x)^2}{4 x^4}\)

\(\Big \downarrow \) 3778

\(\displaystyle \frac {1}{2} b \left (\frac {1}{3} \pi b^2 \int \frac {\cos \left (\frac {1}{2} b^2 \pi x^2\right ) \operatorname {FresnelS}(b x)}{x^2}dx-\frac {1}{12} b \left (\pi b^2 \int -\frac {\sin \left (b^2 \pi x^2\right )}{x^2}dx^2-\frac {\cos \left (\pi b^2 x^2\right )}{x^2}\right )-\frac {\operatorname {FresnelS}(b x) \sin \left (\frac {1}{2} \pi b^2 x^2\right )}{3 x^3}-\frac {b}{12 x^2}\right )-\frac {\operatorname {FresnelS}(b x)^2}{4 x^4}\)

\(\Big \downarrow \) 25

\(\displaystyle \frac {1}{2} b \left (\frac {1}{3} \pi b^2 \int \frac {\cos \left (\frac {1}{2} b^2 \pi x^2\right ) \operatorname {FresnelS}(b x)}{x^2}dx-\frac {1}{12} b \left (-\pi b^2 \int \frac {\sin \left (b^2 \pi x^2\right )}{x^2}dx^2-\frac {\cos \left (\pi b^2 x^2\right )}{x^2}\right )-\frac {\operatorname {FresnelS}(b x) \sin \left (\frac {1}{2} \pi b^2 x^2\right )}{3 x^3}-\frac {b}{12 x^2}\right )-\frac {\operatorname {FresnelS}(b x)^2}{4 x^4}\)

\(\Big \downarrow \) 3042

\(\displaystyle \frac {1}{2} b \left (\frac {1}{3} \pi b^2 \int \frac {\cos \left (\frac {1}{2} b^2 \pi x^2\right ) \operatorname {FresnelS}(b x)}{x^2}dx-\frac {1}{12} b \left (-\pi b^2 \int \frac {\sin \left (b^2 \pi x^2\right )}{x^2}dx^2-\frac {\cos \left (\pi b^2 x^2\right )}{x^2}\right )-\frac {\operatorname {FresnelS}(b x) \sin \left (\frac {1}{2} \pi b^2 x^2\right )}{3 x^3}-\frac {b}{12 x^2}\right )-\frac {\operatorname {FresnelS}(b x)^2}{4 x^4}\)

\(\Big \downarrow \) 3780

\(\displaystyle \frac {1}{2} b \left (\frac {1}{3} \pi b^2 \int \frac {\cos \left (\frac {1}{2} b^2 \pi x^2\right ) \operatorname {FresnelS}(b x)}{x^2}dx-\frac {\operatorname {FresnelS}(b x) \sin \left (\frac {1}{2} \pi b^2 x^2\right )}{3 x^3}-\frac {1}{12} b \left (-\pi b^2 \text {Si}\left (b^2 \pi x^2\right )-\frac {\cos \left (\pi b^2 x^2\right )}{x^2}\right )-\frac {b}{12 x^2}\right )-\frac {\operatorname {FresnelS}(b x)^2}{4 x^4}\)

\(\Big \downarrow \) 7018

\(\displaystyle \frac {1}{2} b \left (\frac {1}{3} \pi b^2 \left (-\pi b^2 \int \operatorname {FresnelS}(b x) \sin \left (\frac {1}{2} b^2 \pi x^2\right )dx+\frac {1}{2} b \int \frac {\sin \left (b^2 \pi x^2\right )}{x}dx-\frac {\operatorname {FresnelS}(b x) \cos \left (\frac {1}{2} \pi b^2 x^2\right )}{x}\right )-\frac {\operatorname {FresnelS}(b x) \sin \left (\frac {1}{2} \pi b^2 x^2\right )}{3 x^3}-\frac {1}{12} b \left (-\pi b^2 \text {Si}\left (b^2 \pi x^2\right )-\frac {\cos \left (\pi b^2 x^2\right )}{x^2}\right )-\frac {b}{12 x^2}\right )-\frac {\operatorname {FresnelS}(b x)^2}{4 x^4}\)

\(\Big \downarrow \) 3856

\(\displaystyle \frac {1}{2} b \left (\frac {1}{3} \pi b^2 \left (-\pi b^2 \int \operatorname {FresnelS}(b x) \sin \left (\frac {1}{2} b^2 \pi x^2\right )dx-\frac {\operatorname {FresnelS}(b x) \cos \left (\frac {1}{2} \pi b^2 x^2\right )}{x}+\frac {1}{4} b \text {Si}\left (b^2 \pi x^2\right )\right )-\frac {\operatorname {FresnelS}(b x) \sin \left (\frac {1}{2} \pi b^2 x^2\right )}{3 x^3}-\frac {1}{12} b \left (-\pi b^2 \text {Si}\left (b^2 \pi x^2\right )-\frac {\cos \left (\pi b^2 x^2\right )}{x^2}\right )-\frac {b}{12 x^2}\right )-\frac {\operatorname {FresnelS}(b x)^2}{4 x^4}\)

\(\Big \downarrow \) 6994

\(\displaystyle \frac {1}{2} b \left (\frac {1}{3} \pi b^2 \left (-\pi b \int \operatorname {FresnelS}(b x)d\operatorname {FresnelS}(b x)-\frac {\operatorname {FresnelS}(b x) \cos \left (\frac {1}{2} \pi b^2 x^2\right )}{x}+\frac {1}{4} b \text {Si}\left (b^2 \pi x^2\right )\right )-\frac {\operatorname {FresnelS}(b x) \sin \left (\frac {1}{2} \pi b^2 x^2\right )}{3 x^3}-\frac {1}{12} b \left (-\pi b^2 \text {Si}\left (b^2 \pi x^2\right )-\frac {\cos \left (\pi b^2 x^2\right )}{x^2}\right )-\frac {b}{12 x^2}\right )-\frac {\operatorname {FresnelS}(b x)^2}{4 x^4}\)

\(\Big \downarrow \) 15

\(\displaystyle \frac {1}{2} b \left (\frac {1}{3} \pi b^2 \left (-\frac {\operatorname {FresnelS}(b x) \cos \left (\frac {1}{2} \pi b^2 x^2\right )}{x}+\frac {1}{4} b \text {Si}\left (b^2 \pi x^2\right )-\frac {1}{2} \pi b \operatorname {FresnelS}(b x)^2\right )-\frac {\operatorname {FresnelS}(b x) \sin \left (\frac {1}{2} \pi b^2 x^2\right )}{3 x^3}-\frac {1}{12} b \left (-\pi b^2 \text {Si}\left (b^2 \pi x^2\right )-\frac {\cos \left (\pi b^2 x^2\right )}{x^2}\right )-\frac {b}{12 x^2}\right )-\frac {\operatorname {FresnelS}(b x)^2}{4 x^4}\)

input

Int[FresnelS[b*x]^2/x^5,x]
 

output

-1/4*FresnelS[b*x]^2/x^4 + (b*(-1/12*b/x^2 - (FresnelS[b*x]*Sin[(b^2*Pi*x^ 
2)/2])/(3*x^3) + (b^2*Pi*(-((Cos[(b^2*Pi*x^2)/2]*FresnelS[b*x])/x) - (b*Pi 
*FresnelS[b*x]^2)/2 + (b*SinIntegral[b^2*Pi*x^2])/4))/3 - (b*(-(Cos[b^2*Pi 
*x^2]/x^2) - b^2*Pi*SinIntegral[b^2*Pi*x^2]))/12))/2
 

3.1.43.3.1 Defintions of rubi rules used

rule 15

Int[(a_.)*(x_)^(m_.), x_Symbol] :> Simp[a*(x^(m + 1)/(m + 1)), x] /; FreeQ[ 
{a, m}, x] && NeQ[m, -1]
 

rule 25

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

rule 3042

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

rule 3778

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] - Simp[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 3780

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

rule 3856

Int[Sin[(d_.)*(x_)^(n_)]/(x_), x_Symbol] :> Simp[SinIntegral[d*x^n]/n, x] / 
; FreeQ[{d, n}, x]
 

rule 3861

Int[((a_.) + Cos[(c_.) + (d_.)*(x_)^(n_)]*(b_.))^(p_.)*(x_)^(m_.), x_Symbol 
] :> Simp[1/n   Subst[Int[x^(Simplify[(m + 1)/n] - 1)*(a + b*Cos[c + d*x])^ 
p, x], x, x^n], x] /; FreeQ[{a, b, c, d, m, n, p}, x] && IntegerQ[Simplify[ 
(m + 1)/n]] && (EqQ[p, 1] || EqQ[m, n - 1] || (IntegerQ[p] && GtQ[Simplify[ 
(m + 1)/n], 0]))
 

rule 6984

Int[FresnelS[(b_.)*(x_)]^2*(x_)^(m_.), x_Symbol] :> Simp[x^(m + 1)*(Fresnel 
S[b*x]^2/(m + 1)), x] - Simp[2*(b/(m + 1))   Int[x^(m + 1)*Sin[(Pi/2)*b^2*x 
^2]*FresnelS[b*x], x], x] /; FreeQ[b, x] && IntegerQ[m] && NeQ[m, -1]
 

rule 6994

Int[FresnelS[(b_.)*(x_)]^(n_.)*Sin[(d_.)*(x_)^2], x_Symbol] :> Simp[Pi*(b/( 
2*d))   Subst[Int[x^n, x], x, FresnelS[b*x]], x] /; FreeQ[{b, d, n}, x] && 
EqQ[d^2, (Pi^2/4)*b^4]
 

rule 7010

Int[FresnelS[(b_.)*(x_)]*(x_)^(m_)*Sin[(d_.)*(x_)^2], x_Symbol] :> Simp[x^( 
m + 1)*Sin[d*x^2]*(FresnelS[b*x]/(m + 1)), x] + (-Simp[d*(x^(m + 2)/(Pi*b*( 
m + 1)*(m + 2))), x] - Simp[2*(d/(m + 1))   Int[x^(m + 2)*Cos[d*x^2]*Fresne 
lS[b*x], x], x] + Simp[d/(Pi*b*(m + 1))   Int[x^(m + 1)*Cos[2*d*x^2], x], x 
]) /; FreeQ[{b, d}, x] && EqQ[d^2, (Pi^2/4)*b^4] && ILtQ[m, -2]
 

rule 7018

Int[Cos[(d_.)*(x_)^2]*FresnelS[(b_.)*(x_)]*(x_)^(m_), x_Symbol] :> Simp[x^( 
m + 1)*Cos[d*x^2]*(FresnelS[b*x]/(m + 1)), x] + (Simp[2*(d/(m + 1))   Int[x 
^(m + 2)*Sin[d*x^2]*FresnelS[b*x], x], x] - Simp[d/(Pi*b*(m + 1))   Int[x^( 
m + 1)*Sin[2*d*x^2], x], x]) /; FreeQ[{b, d}, x] && EqQ[d^2, (Pi^2/4)*b^4] 
&& ILtQ[m, -1]
 

3.1.43.4 Maple [F]

\[\int \frac {\operatorname {FresnelS}\left (b x \right )^{2}}{x^{5}}d x\]

input

int(FresnelS(b*x)^2/x^5,x)
 

output

int(FresnelS(b*x)^2/x^5,x)
 

3.1.43.5 Fricas [A] (verification not implemented)

Time = 0.26 (sec) , antiderivative size = 111, normalized size of antiderivative = 0.87 \[ \int \frac {\operatorname {FresnelS}(b x)^2}{x^5} \, dx=\frac {\pi b^{4} x^{4} \operatorname {Si}\left (\pi b^{2} x^{2}\right ) - 2 \, \pi b^{3} x^{3} \cos \left (\frac {1}{2} \, \pi b^{2} x^{2}\right ) \operatorname {S}\left (b x\right ) + b^{2} x^{2} \cos \left (\frac {1}{2} \, \pi b^{2} x^{2}\right )^{2} - b^{2} x^{2} - 2 \, b x \operatorname {S}\left (b x\right ) \sin \left (\frac {1}{2} \, \pi b^{2} x^{2}\right ) - {\left (\pi ^{2} b^{4} x^{4} + 3\right )} \operatorname {S}\left (b x\right )^{2}}{12 \, x^{4}} \]

input

integrate(fresnel_sin(b*x)^2/x^5,x, algorithm="fricas")
 

output

1/12*(pi*b^4*x^4*sin_integral(pi*b^2*x^2) - 2*pi*b^3*x^3*cos(1/2*pi*b^2*x^ 
2)*fresnel_sin(b*x) + b^2*x^2*cos(1/2*pi*b^2*x^2)^2 - b^2*x^2 - 2*b*x*fres 
nel_sin(b*x)*sin(1/2*pi*b^2*x^2) - (pi^2*b^4*x^4 + 3)*fresnel_sin(b*x)^2)/ 
x^4
 

3.1.43.6 Sympy [F]

\[ \int \frac {\operatorname {FresnelS}(b x)^2}{x^5} \, dx=\int \frac {S^{2}\left (b x\right )}{x^{5}}\, dx \]

input

integrate(fresnels(b*x)**2/x**5,x)
 

output

Integral(fresnels(b*x)**2/x**5, x)
 

3.1.43.7 Maxima [F]

\[ \int \frac {\operatorname {FresnelS}(b x)^2}{x^5} \, dx=\int { \frac {\operatorname {S}\left (b x\right )^{2}}{x^{5}} \,d x } \]

input

integrate(fresnel_sin(b*x)^2/x^5,x, algorithm="maxima")
 

output

integrate(fresnel_sin(b*x)^2/x^5, x)
 

3.1.43.8 Giac [F]

\[ \int \frac {\operatorname {FresnelS}(b x)^2}{x^5} \, dx=\int { \frac {\operatorname {S}\left (b x\right )^{2}}{x^{5}} \,d x } \]

input

integrate(fresnel_sin(b*x)^2/x^5,x, algorithm="giac")
 

output

integrate(fresnel_sin(b*x)^2/x^5, x)
 

3.1.43.9 Mupad [F(-1)]

Timed out. \[ \int \frac {\operatorname {FresnelS}(b x)^2}{x^5} \, dx=\int \frac {{\mathrm {FresnelS}\left (b\,x\right )}^2}{x^5} \,d x \]

input

int(FresnelS(b*x)^2/x^5,x)
 

output

int(FresnelS(b*x)^2/x^5, x)