\(\int x^6 \operatorname {FresnelS}(b x) \, dx\) [2]

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

Optimal result

Integrand size = 8, antiderivative size = 109 \[ \int x^6 \operatorname {FresnelS}(b x) \, dx=-\frac {24 x^2 \cos \left (\frac {1}{2} b^2 \pi x^2\right )}{7 b^5 \pi ^3}+\frac {x^6 \cos \left (\frac {1}{2} b^2 \pi x^2\right )}{7 b \pi }+\frac {1}{7} x^7 \operatorname {FresnelS}(b x)+\frac {48 \sin \left (\frac {1}{2} b^2 \pi x^2\right )}{7 b^7 \pi ^4}-\frac {6 x^4 \sin \left (\frac {1}{2} b^2 \pi x^2\right )}{7 b^3 \pi ^2} \] Output:

-24/7*x^2*cos(1/2*b^2*Pi*x^2)/b^5/Pi^3+1/7*x^6*cos(1/2*b^2*Pi*x^2)/b/Pi+1/ 
7*x^7*FresnelS(b*x)+48/7*sin(1/2*b^2*Pi*x^2)/b^7/Pi^4-6/7*x^4*sin(1/2*b^2* 
Pi*x^2)/b^3/Pi^2
 

Mathematica [A] (verified)

Time = 0.04 (sec) , antiderivative size = 83, normalized size of antiderivative = 0.76 \[ \int x^6 \operatorname {FresnelS}(b x) \, dx=\frac {x^2 \left (-24+b^4 \pi ^2 x^4\right ) \cos \left (\frac {1}{2} b^2 \pi x^2\right )}{7 b^5 \pi ^3}+\frac {1}{7} x^7 \operatorname {FresnelS}(b x)-\frac {6 \left (-8+b^4 \pi ^2 x^4\right ) \sin \left (\frac {1}{2} b^2 \pi x^2\right )}{7 b^7 \pi ^4} \] Input:

Integrate[x^6*FresnelS[b*x],x]
 

Output:

(x^2*(-24 + b^4*Pi^2*x^4)*Cos[(b^2*Pi*x^2)/2])/(7*b^5*Pi^3) + (x^7*Fresnel 
S[b*x])/7 - (6*(-8 + b^4*Pi^2*x^4)*Sin[(b^2*Pi*x^2)/2])/(7*b^7*Pi^4)
 

Rubi [A] (verified)

Time = 0.60 (sec) , antiderivative size = 125, normalized size of antiderivative = 1.15, number of steps used = 12, number of rules used = 11, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 1.375, Rules used = {6980, 3860, 3042, 3777, 3042, 3777, 25, 3042, 3777, 3042, 3117}

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

\(\Big \downarrow \) 6980

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

\(\Big \downarrow \) 3860

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

\(\Big \downarrow \) 3042

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

\(\Big \downarrow \) 3777

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

\(\Big \downarrow \) 3042

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

\(\Big \downarrow \) 3777

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

\(\Big \downarrow \) 25

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

\(\Big \downarrow \) 3042

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

\(\Big \downarrow \) 3777

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

\(\Big \downarrow \) 3042

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

\(\Big \downarrow \) 3117

\(\displaystyle \frac {1}{7} x^7 \operatorname {FresnelS}(b x)-\frac {1}{14} b \left (\frac {6 \left (\frac {2 x^4 \sin \left (\frac {1}{2} \pi b^2 x^2\right )}{\pi b^2}-\frac {4 \left (\frac {4 \sin \left (\frac {1}{2} \pi b^2 x^2\right )}{\pi ^2 b^4}-\frac {2 x^2 \cos \left (\frac {1}{2} \pi b^2 x^2\right )}{\pi b^2}\right )}{\pi b^2}\right )}{\pi b^2}-\frac {2 x^6 \cos \left (\frac {1}{2} \pi b^2 x^2\right )}{\pi b^2}\right )\)

Input:

Int[x^6*FresnelS[b*x],x]
 

Output:

(x^7*FresnelS[b*x])/7 - (b*((-2*x^6*Cos[(b^2*Pi*x^2)/2])/(b^2*Pi) + (6*((2 
*x^4*Sin[(b^2*Pi*x^2)/2])/(b^2*Pi) - (4*((-2*x^2*Cos[(b^2*Pi*x^2)/2])/(b^2 
*Pi) + (4*Sin[(b^2*Pi*x^2)/2])/(b^4*Pi^2)))/(b^2*Pi)))/(b^2*Pi)))/14
 

Defintions of rubi rules used

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 3117
Int[sin[Pi/2 + (c_.) + (d_.)*(x_)], x_Symbol] :> Simp[Sin[c + d*x]/d, x] /; 
 FreeQ[{c, d}, x]
 

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

rule 3860
Int[(x_)^(m_.)*((a_.) + (b_.)*Sin[(c_.) + (d_.)*(x_)^(n_)])^(p_.), x_Symbol 
] :> Simp[1/n   Subst[Int[x^(Simplify[(m + 1)/n] - 1)*(a + b*Sin[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 6980
Int[FresnelS[(b_.)*(x_)]*((d_.)*(x_))^(m_.), x_Symbol] :> Simp[(d*x)^(m + 1 
)*(FresnelS[b*x]/(d*(m + 1))), x] - Simp[b/(d*(m + 1))   Int[(d*x)^(m + 1)* 
Sin[(Pi/2)*b^2*x^2], x], x] /; FreeQ[{b, d, m}, x] && NeQ[m, -1]
 
Maple [C] (verified)

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

Time = 0.48 (sec) , antiderivative size = 29, normalized size of antiderivative = 0.27

method result size
meijerg \(\frac {\pi \,b^{3} x^{10} \operatorname {hypergeom}\left (\left [\frac {3}{4}, \frac {5}{2}\right ], \left [\frac {3}{2}, \frac {7}{4}, \frac {7}{2}\right ], -\frac {x^{4} \pi ^{2} b^{4}}{16}\right )}{60}\) \(29\)
derivativedivides \(\frac {\frac {\operatorname {FresnelS}\left (b x \right ) b^{7} x^{7}}{7}+\frac {b^{6} x^{6} \cos \left (\frac {b^{2} \pi \,x^{2}}{2}\right )}{7 \pi }-\frac {6 \left (\frac {b^{4} x^{4} \sin \left (\frac {b^{2} \pi \,x^{2}}{2}\right )}{\pi }-\frac {4 \left (-\frac {b^{2} x^{2} \cos \left (\frac {b^{2} \pi \,x^{2}}{2}\right )}{\pi }+\frac {2 \sin \left (\frac {b^{2} \pi \,x^{2}}{2}\right )}{\pi ^{2}}\right )}{\pi }\right )}{7 \pi }}{b^{7}}\) \(107\)
default \(\frac {\frac {\operatorname {FresnelS}\left (b x \right ) b^{7} x^{7}}{7}+\frac {b^{6} x^{6} \cos \left (\frac {b^{2} \pi \,x^{2}}{2}\right )}{7 \pi }-\frac {6 \left (\frac {b^{4} x^{4} \sin \left (\frac {b^{2} \pi \,x^{2}}{2}\right )}{\pi }-\frac {4 \left (-\frac {b^{2} x^{2} \cos \left (\frac {b^{2} \pi \,x^{2}}{2}\right )}{\pi }+\frac {2 \sin \left (\frac {b^{2} \pi \,x^{2}}{2}\right )}{\pi ^{2}}\right )}{\pi }\right )}{7 \pi }}{b^{7}}\) \(107\)
parts \(\frac {x^{7} \operatorname {FresnelS}\left (b x \right )}{7}-\frac {b \left (-\frac {x^{6} \cos \left (\frac {b^{2} \pi \,x^{2}}{2}\right )}{b^{2} \pi }+\frac {\frac {6 x^{4} \sin \left (\frac {b^{2} \pi \,x^{2}}{2}\right )}{b^{2} \pi }-\frac {24 \left (-\frac {x^{2} \cos \left (\frac {b^{2} \pi \,x^{2}}{2}\right )}{b^{2} \pi }+\frac {2 \sin \left (\frac {b^{2} \pi \,x^{2}}{2}\right )}{b^{4} \pi ^{2}}\right )}{b^{2} \pi }}{b^{2} \pi }\right )}{7}\) \(113\)

Input:

int(x^6*FresnelS(b*x),x,method=_RETURNVERBOSE)
 

Output:

1/60*Pi*b^3*x^10*hypergeom([3/4,5/2],[3/2,7/4,7/2],-1/16*x^4*Pi^2*b^4)
 

Fricas [A] (verification not implemented)

Time = 0.09 (sec) , antiderivative size = 78, normalized size of antiderivative = 0.72 \[ \int x^6 \operatorname {FresnelS}(b x) \, dx=\frac {\pi ^{4} b^{7} x^{7} \operatorname {S}\left (b x\right ) + {\left (\pi ^{3} b^{6} x^{6} - 24 \, \pi b^{2} x^{2}\right )} \cos \left (\frac {1}{2} \, \pi b^{2} x^{2}\right ) - 6 \, {\left (\pi ^{2} b^{4} x^{4} - 8\right )} \sin \left (\frac {1}{2} \, \pi b^{2} x^{2}\right )}{7 \, \pi ^{4} b^{7}} \] Input:

integrate(x^6*fresnel_sin(b*x),x, algorithm="fricas")
 

Output:

1/7*(pi^4*b^7*x^7*fresnel_sin(b*x) + (pi^3*b^6*x^6 - 24*pi*b^2*x^2)*cos(1/ 
2*pi*b^2*x^2) - 6*(pi^2*b^4*x^4 - 8)*sin(1/2*pi*b^2*x^2))/(pi^4*b^7)
 

Sympy [A] (verification not implemented)

Time = 0.99 (sec) , antiderivative size = 156, normalized size of antiderivative = 1.43 \[ \int x^6 \operatorname {FresnelS}(b x) \, dx=\frac {3 x^{7} S\left (b x\right ) \Gamma \left (\frac {3}{4}\right )}{28 \Gamma \left (\frac {7}{4}\right )} + \frac {3 x^{6} \cos {\left (\frac {\pi b^{2} x^{2}}{2} \right )} \Gamma \left (\frac {3}{4}\right )}{28 \pi b \Gamma \left (\frac {7}{4}\right )} - \frac {9 x^{4} \sin {\left (\frac {\pi b^{2} x^{2}}{2} \right )} \Gamma \left (\frac {3}{4}\right )}{14 \pi ^{2} b^{3} \Gamma \left (\frac {7}{4}\right )} - \frac {18 x^{2} \cos {\left (\frac {\pi b^{2} x^{2}}{2} \right )} \Gamma \left (\frac {3}{4}\right )}{7 \pi ^{3} b^{5} \Gamma \left (\frac {7}{4}\right )} + \frac {36 \sin {\left (\frac {\pi b^{2} x^{2}}{2} \right )} \Gamma \left (\frac {3}{4}\right )}{7 \pi ^{4} b^{7} \Gamma \left (\frac {7}{4}\right )} \] Input:

integrate(x**6*fresnels(b*x),x)
 

Output:

3*x**7*fresnels(b*x)*gamma(3/4)/(28*gamma(7/4)) + 3*x**6*cos(pi*b**2*x**2/ 
2)*gamma(3/4)/(28*pi*b*gamma(7/4)) - 9*x**4*sin(pi*b**2*x**2/2)*gamma(3/4) 
/(14*pi**2*b**3*gamma(7/4)) - 18*x**2*cos(pi*b**2*x**2/2)*gamma(3/4)/(7*pi 
**3*b**5*gamma(7/4)) + 36*sin(pi*b**2*x**2/2)*gamma(3/4)/(7*pi**4*b**7*gam 
ma(7/4))
 

Maxima [A] (verification not implemented)

Time = 0.05 (sec) , antiderivative size = 74, normalized size of antiderivative = 0.68 \[ \int x^6 \operatorname {FresnelS}(b x) \, dx=\frac {1}{7} \, x^{7} \operatorname {S}\left (b x\right ) + \frac {{\left (\pi ^{3} b^{6} x^{6} - 24 \, \pi b^{2} x^{2}\right )} \cos \left (\frac {1}{2} \, \pi b^{2} x^{2}\right ) - 6 \, {\left (\pi ^{2} b^{4} x^{4} - 8\right )} \sin \left (\frac {1}{2} \, \pi b^{2} x^{2}\right )}{7 \, \pi ^{4} b^{7}} \] Input:

integrate(x^6*fresnel_sin(b*x),x, algorithm="maxima")
                                                                                    
                                                                                    
 

Output:

1/7*x^7*fresnel_sin(b*x) + 1/7*((pi^3*b^6*x^6 - 24*pi*b^2*x^2)*cos(1/2*pi* 
b^2*x^2) - 6*(pi^2*b^4*x^4 - 8)*sin(1/2*pi*b^2*x^2))/(pi^4*b^7)
 

Giac [F]

\[ \int x^6 \operatorname {FresnelS}(b x) \, dx=\int { x^{6} \operatorname {S}\left (b x\right ) \,d x } \] Input:

integrate(x^6*fresnel_sin(b*x),x, algorithm="giac")
 

Output:

integrate(x^6*fresnel_sin(b*x), x)
 

Mupad [F(-1)]

Timed out. \[ \int x^6 \operatorname {FresnelS}(b x) \, dx=\int x^6\,\mathrm {FresnelS}\left (b\,x\right ) \,d x \] Input:

int(x^6*FresnelS(b*x),x)
 

Output:

int(x^6*FresnelS(b*x), x)
 

Reduce [F]

\[ \int x^6 \operatorname {FresnelS}(b x) \, dx=\int x^{6} \mathrm {FresnelS}\left (b x \right )d x \] Input:

int(x^6*FresnelS(b*x),x)
 

Output:

int(x^6*FresnelS(b*x),x)