\(\int \frac {\operatorname {FresnelC}(b x)^2}{x^5} \, dx\) [152]

Optimal result

Integrand size = 10, antiderivative size = 127 \[ \int \frac {\operatorname {FresnelC}(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 \cos \left (\frac {1}{2} b^2 \pi x^2\right ) \operatorname {FresnelC}(b x)}{6 x^3}-\frac {1}{12} b^4 \pi ^2 \operatorname {FresnelC}(b x)^2-\frac {\operatorname {FresnelC}(b x)^2}{4 x^4}+\frac {b^3 \pi \operatorname {FresnelC}(b x) \sin \left (\frac {1}{2} b^2 \pi x^2\right )}{6 x}-\frac {1}{12} b^4 \pi \text {Si}\left (b^2 \pi x^2\right ) \]

[Out]

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

Rubi [A] (verified)

Time = 0.10 (sec) , antiderivative size = 127, normalized size of antiderivative = 1.00, number of steps used = 9, number of rules used = 9, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.900, Rules used = {6566, 6592, 6600, 6576, 30, 3456, 3461, 3378, 3380} \[ \int \frac {\operatorname {FresnelC}(b x)^2}{x^5} \, dx=-\frac {1}{12} \pi ^2 b^4 \operatorname {FresnelC}(b x)^2-\frac {b \operatorname {FresnelC}(b x) \cos \left (\frac {1}{2} \pi b^2 x^2\right )}{6 x^3}-\frac {b^2}{24 x^2}-\frac {b^2 \cos \left (\pi b^2 x^2\right )}{24 x^2}-\frac {1}{12} \pi b^4 \text {Si}\left (b^2 \pi x^2\right )+\frac {\pi b^3 \operatorname {FresnelC}(b x) \sin \left (\frac {1}{2} \pi b^2 x^2\right )}{6 x}-\frac {\operatorname {FresnelC}(b x)^2}{4 x^4} \]

[In]

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

[Out]

-1/24*b^2/x^2 - (b^2*Cos[b^2*Pi*x^2])/(24*x^2) - (b*Cos[(b^2*Pi*x^2)/2]*FresnelC[b*x])/(6*x^3) - (b^4*Pi^2*Fre
snelC[b*x]^2)/12 - FresnelC[b*x]^2/(4*x^4) + (b^3*Pi*FresnelC[b*x]*Sin[(b^2*Pi*x^2)/2])/(6*x) - (b^4*Pi*SinInt
egral[b^2*Pi*x^2])/12

Rule 30

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

Rule 3378

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 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 3456

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

Rule 3461

Int[((a_.) + Cos[(c_.) + (d_.)*(x_)^(n_)]*(b_.))^(p_.)*(x_)^(m_.), x_Symbol] :> Dist[1/n, Subst[Int[x^(Simplif
y[(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[Simpl
ify[(m + 1)/n]] && (EqQ[p, 1] || EqQ[m, n - 1] || (IntegerQ[p] && GtQ[Simplify[(m + 1)/n], 0]))

Rule 6566

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

Rule 6576

Int[Cos[(d_.)*(x_)^2]*FresnelC[(b_.)*(x_)]^(n_.), x_Symbol] :> Dist[Pi*(b/(2*d)), Subst[Int[x^n, x], x, Fresne
lC[b*x]], x] /; FreeQ[{b, d, n}, x] && EqQ[d^2, (Pi^2/4)*b^4]

Rule 6592

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

Rule 6600

Int[FresnelC[(b_.)*(x_)]*(x_)^(m_)*Sin[(d_.)*(x_)^2], x_Symbol] :> Simp[x^(m + 1)*Sin[d*x^2]*(FresnelC[b*x]/(m
+ 1)), x] + (-Dist[2*(d/(m + 1)), Int[x^(m + 2)*Cos[d*x^2]*FresnelC[b*x], x], x] - Dist[b/(2*(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]

Rubi steps \begin{align*} \text {integral}& = -\frac {\operatorname {FresnelC}(b x)^2}{4 x^4}+\frac {1}{2} b \int \frac {\cos \left (\frac {1}{2} b^2 \pi x^2\right ) \operatorname {FresnelC}(b x)}{x^4} \, dx \\ & = -\frac {b^2}{24 x^2}-\frac {b \cos \left (\frac {1}{2} b^2 \pi x^2\right ) \operatorname {FresnelC}(b x)}{6 x^3}-\frac {\operatorname {FresnelC}(b x)^2}{4 x^4}+\frac {1}{12} b^2 \int \frac {\cos \left (b^2 \pi x^2\right )}{x^3} \, dx-\frac {1}{6} \left (b^3 \pi \right ) \int \frac {\operatorname {FresnelC}(b x) \sin \left (\frac {1}{2} b^2 \pi x^2\right )}{x^2} \, dx \\ & = -\frac {b^2}{24 x^2}-\frac {b \cos \left (\frac {1}{2} b^2 \pi x^2\right ) \operatorname {FresnelC}(b x)}{6 x^3}-\frac {\operatorname {FresnelC}(b x)^2}{4 x^4}+\frac {b^3 \pi \operatorname {FresnelC}(b x) \sin \left (\frac {1}{2} b^2 \pi x^2\right )}{6 x}+\frac {1}{24} b^2 \text {Subst}\left (\int \frac {\cos \left (b^2 \pi x\right )}{x^2} \, dx,x,x^2\right )-\frac {1}{12} \left (b^4 \pi \right ) \int \frac {\sin \left (b^2 \pi x^2\right )}{x} \, dx-\frac {1}{6} \left (b^5 \pi ^2\right ) \int \cos \left (\frac {1}{2} b^2 \pi x^2\right ) \operatorname {FresnelC}(b x) \, dx \\ & = -\frac {b^2}{24 x^2}-\frac {b^2 \cos \left (b^2 \pi x^2\right )}{24 x^2}-\frac {b \cos \left (\frac {1}{2} b^2 \pi x^2\right ) \operatorname {FresnelC}(b x)}{6 x^3}-\frac {\operatorname {FresnelC}(b x)^2}{4 x^4}+\frac {b^3 \pi \operatorname {FresnelC}(b x) \sin \left (\frac {1}{2} b^2 \pi x^2\right )}{6 x}-\frac {1}{24} b^4 \pi \text {Si}\left (b^2 \pi x^2\right )-\frac {1}{24} \left (b^4 \pi \right ) \text {Subst}\left (\int \frac {\sin \left (b^2 \pi x\right )}{x} \, dx,x,x^2\right )-\frac {1}{6} \left (b^4 \pi ^2\right ) \text {Subst}(\int x \, dx,x,\operatorname {FresnelC}(b x)) \\ & = -\frac {b^2}{24 x^2}-\frac {b^2 \cos \left (b^2 \pi x^2\right )}{24 x^2}-\frac {b \cos \left (\frac {1}{2} b^2 \pi x^2\right ) \operatorname {FresnelC}(b x)}{6 x^3}-\frac {1}{12} b^4 \pi ^2 \operatorname {FresnelC}(b x)^2-\frac {\operatorname {FresnelC}(b x)^2}{4 x^4}+\frac {b^3 \pi \operatorname {FresnelC}(b x) \sin \left (\frac {1}{2} b^2 \pi x^2\right )}{6 x}-\frac {1}{12} b^4 \pi \text {Si}\left (b^2 \pi x^2\right ) \\ \end{align*}

Mathematica [A] (verified)

Time = 0.01 (sec) , antiderivative size = 127, normalized size of antiderivative = 1.00 \[ \int \frac {\operatorname {FresnelC}(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 \cos \left (\frac {1}{2} b^2 \pi x^2\right ) \operatorname {FresnelC}(b x)}{6 x^3}-\frac {1}{12} b^4 \pi ^2 \operatorname {FresnelC}(b x)^2-\frac {\operatorname {FresnelC}(b x)^2}{4 x^4}+\frac {b^3 \pi \operatorname {FresnelC}(b x) \sin \left (\frac {1}{2} b^2 \pi x^2\right )}{6 x}-\frac {1}{12} b^4 \pi \text {Si}\left (b^2 \pi x^2\right ) \]

[In]

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

[Out]

-1/24*b^2/x^2 - (b^2*Cos[b^2*Pi*x^2])/(24*x^2) - (b*Cos[(b^2*Pi*x^2)/2]*FresnelC[b*x])/(6*x^3) - (b^4*Pi^2*Fre
snelC[b*x]^2)/12 - FresnelC[b*x]^2/(4*x^4) + (b^3*Pi*FresnelC[b*x]*Sin[(b^2*Pi*x^2)/2])/(6*x) - (b^4*Pi*SinInt
egral[b^2*Pi*x^2])/12

Maple [F]

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

[In]

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

[Out]

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

Fricas [A] (verification not implemented)

none

Time = 0.26 (sec) , antiderivative size = 102, normalized size of antiderivative = 0.80 \[ \int \frac {\operatorname {FresnelC}(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} \operatorname {C}\left (b x\right ) \sin \left (\frac {1}{2} \, \pi b^{2} x^{2}\right ) + b^{2} x^{2} \cos \left (\frac {1}{2} \, \pi b^{2} x^{2}\right )^{2} + 2 \, b x \cos \left (\frac {1}{2} \, \pi b^{2} x^{2}\right ) \operatorname {C}\left (b x\right ) + {\left (\pi ^{2} b^{4} x^{4} + 3\right )} \operatorname {C}\left (b x\right )^{2}}{12 \, x^{4}} \]

[In]

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

[Out]

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

Sympy [F]

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

[In]

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

[Out]

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

Maxima [F]

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

[In]

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

[Out]

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

Giac [F]

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

[In]

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

[Out]

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

Mupad [F(-1)]

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

[In]

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

[Out]

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