3.2.7 \(\int (d x)^m \operatorname {PolyLog}(4,a x^2) \, dx\) [107]

3.2.7.1 Optimal result
3.2.7.2 Mathematica [C] (verified)
3.2.7.3 Rubi [A] (verified)
3.2.7.4 Maple [C] (verified)
3.2.7.5 Fricas [F]
3.2.7.6 Sympy [F]
3.2.7.7 Maxima [F]
3.2.7.8 Giac [F]
3.2.7.9 Mupad [F(-1)]

3.2.7.1 Optimal result

Integrand size = 13, antiderivative size = 142 \[ \int (d x)^m \operatorname {PolyLog}\left (4,a x^2\right ) \, dx=\frac {16 a (d x)^{3+m} \operatorname {Hypergeometric2F1}\left (1,\frac {3+m}{2},\frac {5+m}{2},a x^2\right )}{d^3 (1+m)^4 (3+m)}+\frac {8 (d x)^{1+m} \log \left (1-a x^2\right )}{d (1+m)^4}+\frac {4 (d x)^{1+m} \operatorname {PolyLog}\left (2,a x^2\right )}{d (1+m)^3}-\frac {2 (d x)^{1+m} \operatorname {PolyLog}\left (3,a x^2\right )}{d (1+m)^2}+\frac {(d x)^{1+m} \operatorname {PolyLog}\left (4,a x^2\right )}{d (1+m)} \]

output
16*a*(d*x)^(3+m)*hypergeom([1, 3/2+1/2*m],[5/2+1/2*m],a*x^2)/d^3/(1+m)^4/( 
3+m)+8*(d*x)^(1+m)*ln(-a*x^2+1)/d/(1+m)^4+4*(d*x)^(1+m)*polylog(2,a*x^2)/d 
/(1+m)^3-2*(d*x)^(1+m)*polylog(3,a*x^2)/d/(1+m)^2+(d*x)^(1+m)*polylog(4,a* 
x^2)/d/(1+m)
 
3.2.7.2 Mathematica [C] (verified)

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

Time = 0.08 (sec) , antiderivative size = 166, normalized size of antiderivative = 1.17 \[ \int (d x)^m \operatorname {PolyLog}\left (4,a x^2\right ) \, dx=\frac {2 x (d x)^m \operatorname {Gamma}\left (\frac {3+m}{2}\right ) \left (4 a (1+m) x^2 \operatorname {Gamma}\left (\frac {1+m}{2}\right ) \, _2\tilde {F}_1\left (1,\frac {3+m}{2};\frac {5+m}{2};a x^2\right )+8 \log \left (1-a x^2\right )+4 (1+m) \operatorname {PolyLog}\left (2,a x^2\right )-2 \operatorname {PolyLog}\left (3,a x^2\right )-4 m \operatorname {PolyLog}\left (3,a x^2\right )-2 m^2 \operatorname {PolyLog}\left (3,a x^2\right )+\operatorname {PolyLog}\left (4,a x^2\right )+3 m \operatorname {PolyLog}\left (4,a x^2\right )+3 m^2 \operatorname {PolyLog}\left (4,a x^2\right )+m^3 \operatorname {PolyLog}\left (4,a x^2\right )\right )}{(1+m)^5 \operatorname {Gamma}\left (\frac {1+m}{2}\right )} \]

input
Integrate[(d*x)^m*PolyLog[4, a*x^2],x]
 
output
(2*x*(d*x)^m*Gamma[(3 + m)/2]*(4*a*(1 + m)*x^2*Gamma[(1 + m)/2]*Hypergeome 
tricPFQRegularized[{1, (3 + m)/2}, {(5 + m)/2}, a*x^2] + 8*Log[1 - a*x^2] 
+ 4*(1 + m)*PolyLog[2, a*x^2] - 2*PolyLog[3, a*x^2] - 4*m*PolyLog[3, a*x^2 
] - 2*m^2*PolyLog[3, a*x^2] + PolyLog[4, a*x^2] + 3*m*PolyLog[4, a*x^2] + 
3*m^2*PolyLog[4, a*x^2] + m^3*PolyLog[4, a*x^2]))/((1 + m)^5*Gamma[(1 + m) 
/2])
 
3.2.7.3 Rubi [A] (verified)

Time = 0.48 (sec) , antiderivative size = 163, normalized size of antiderivative = 1.15, number of steps used = 7, number of rules used = 7, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.538, Rules used = {7145, 7145, 7145, 25, 2905, 8, 278}

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 \operatorname {PolyLog}\left (4,a x^2\right ) (d x)^m \, dx\)

\(\Big \downarrow \) 7145

\(\displaystyle \frac {\operatorname {PolyLog}\left (4,a x^2\right ) (d x)^{m+1}}{d (m+1)}-\frac {2 \int (d x)^m \operatorname {PolyLog}\left (3,a x^2\right )dx}{m+1}\)

\(\Big \downarrow \) 7145

\(\displaystyle \frac {\operatorname {PolyLog}\left (4,a x^2\right ) (d x)^{m+1}}{d (m+1)}-\frac {2 \left (\frac {\operatorname {PolyLog}\left (3,a x^2\right ) (d x)^{m+1}}{d (m+1)}-\frac {2 \int (d x)^m \operatorname {PolyLog}\left (2,a x^2\right )dx}{m+1}\right )}{m+1}\)

\(\Big \downarrow \) 7145

\(\displaystyle \frac {\operatorname {PolyLog}\left (4,a x^2\right ) (d x)^{m+1}}{d (m+1)}-\frac {2 \left (\frac {\operatorname {PolyLog}\left (3,a x^2\right ) (d x)^{m+1}}{d (m+1)}-\frac {2 \left (\frac {\operatorname {PolyLog}\left (2,a x^2\right ) (d x)^{m+1}}{d (m+1)}-\frac {2 \int -(d x)^m \log \left (1-a x^2\right )dx}{m+1}\right )}{m+1}\right )}{m+1}\)

\(\Big \downarrow \) 25

\(\displaystyle \frac {\operatorname {PolyLog}\left (4,a x^2\right ) (d x)^{m+1}}{d (m+1)}-\frac {2 \left (\frac {\operatorname {PolyLog}\left (3,a x^2\right ) (d x)^{m+1}}{d (m+1)}-\frac {2 \left (\frac {2 \int (d x)^m \log \left (1-a x^2\right )dx}{m+1}+\frac {\operatorname {PolyLog}\left (2,a x^2\right ) (d x)^{m+1}}{d (m+1)}\right )}{m+1}\right )}{m+1}\)

\(\Big \downarrow \) 2905

\(\displaystyle \frac {\operatorname {PolyLog}\left (4,a x^2\right ) (d x)^{m+1}}{d (m+1)}-\frac {2 \left (\frac {\operatorname {PolyLog}\left (3,a x^2\right ) (d x)^{m+1}}{d (m+1)}-\frac {2 \left (\frac {2 \left (\frac {2 a \int \frac {x (d x)^{m+1}}{1-a x^2}dx}{d (m+1)}+\frac {\log \left (1-a x^2\right ) (d x)^{m+1}}{d (m+1)}\right )}{m+1}+\frac {\operatorname {PolyLog}\left (2,a x^2\right ) (d x)^{m+1}}{d (m+1)}\right )}{m+1}\right )}{m+1}\)

\(\Big \downarrow \) 8

\(\displaystyle \frac {\operatorname {PolyLog}\left (4,a x^2\right ) (d x)^{m+1}}{d (m+1)}-\frac {2 \left (\frac {\operatorname {PolyLog}\left (3,a x^2\right ) (d x)^{m+1}}{d (m+1)}-\frac {2 \left (\frac {2 \left (\frac {2 a \int \frac {(d x)^{m+2}}{1-a x^2}dx}{d^2 (m+1)}+\frac {\log \left (1-a x^2\right ) (d x)^{m+1}}{d (m+1)}\right )}{m+1}+\frac {\operatorname {PolyLog}\left (2,a x^2\right ) (d x)^{m+1}}{d (m+1)}\right )}{m+1}\right )}{m+1}\)

\(\Big \downarrow \) 278

\(\displaystyle \frac {\operatorname {PolyLog}\left (4,a x^2\right ) (d x)^{m+1}}{d (m+1)}-\frac {2 \left (\frac {\operatorname {PolyLog}\left (3,a x^2\right ) (d x)^{m+1}}{d (m+1)}-\frac {2 \left (\frac {2 \left (\frac {2 a (d x)^{m+3} \operatorname {Hypergeometric2F1}\left (1,\frac {m+3}{2},\frac {m+5}{2},a x^2\right )}{d^3 (m+1) (m+3)}+\frac {\log \left (1-a x^2\right ) (d x)^{m+1}}{d (m+1)}\right )}{m+1}+\frac {\operatorname {PolyLog}\left (2,a x^2\right ) (d x)^{m+1}}{d (m+1)}\right )}{m+1}\right )}{m+1}\)

input
Int[(d*x)^m*PolyLog[4, a*x^2],x]
 
output
(-2*((-2*((2*((2*a*(d*x)^(3 + m)*Hypergeometric2F1[1, (3 + m)/2, (5 + m)/2 
, a*x^2])/(d^3*(1 + m)*(3 + m)) + ((d*x)^(1 + m)*Log[1 - a*x^2])/(d*(1 + m 
))))/(1 + m) + ((d*x)^(1 + m)*PolyLog[2, a*x^2])/(d*(1 + m))))/(1 + m) + ( 
(d*x)^(1 + m)*PolyLog[3, a*x^2])/(d*(1 + m))))/(1 + m) + ((d*x)^(1 + m)*Po 
lyLog[4, a*x^2])/(d*(1 + m))
 

3.2.7.3.1 Defintions of rubi rules used

rule 8
Int[(u_.)*(x_)^(m_.)*((a_.)*(x_))^(p_), x_Symbol] :> Simp[1/a^m   Int[u*(a* 
x)^(m + p), x], x] /; FreeQ[{a, m, p}, x] && IntegerQ[m]
 

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

rule 278
Int[((c_.)*(x_))^(m_.)*((a_) + (b_.)*(x_)^2)^(p_), x_Symbol] :> Simp[a^p*(( 
c*x)^(m + 1)/(c*(m + 1)))*Hypergeometric2F1[-p, (m + 1)/2, (m + 1)/2 + 1, ( 
-b)*(x^2/a)], x] /; FreeQ[{a, b, c, m, p}, x] &&  !IGtQ[p, 0] && (ILtQ[p, 0 
] || GtQ[a, 0])
 

rule 2905
Int[((a_.) + Log[(c_.)*((d_) + (e_.)*(x_)^(n_))^(p_.)]*(b_.))*((f_.)*(x_))^ 
(m_.), x_Symbol] :> Simp[(f*x)^(m + 1)*((a + b*Log[c*(d + e*x^n)^p])/(f*(m 
+ 1))), x] - Simp[b*e*n*(p/(f*(m + 1)))   Int[x^(n - 1)*((f*x)^(m + 1)/(d + 
 e*x^n)), x], x] /; FreeQ[{a, b, c, d, e, f, m, n, p}, x] && NeQ[m, -1]
 

rule 7145
Int[((d_.)*(x_))^(m_.)*PolyLog[n_, (a_.)*((b_.)*(x_)^(p_.))^(q_.)], x_Symbo 
l] :> Simp[(d*x)^(m + 1)*(PolyLog[n, a*(b*x^p)^q]/(d*(m + 1))), x] - Simp[p 
*(q/(m + 1))   Int[(d*x)^m*PolyLog[n - 1, a*(b*x^p)^q], x], x] /; FreeQ[{a, 
 b, d, m, p, q}, x] && NeQ[m, -1] && GtQ[n, 0]
 
3.2.7.4 Maple [C] (verified)

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

Time = 0.68 (sec) , antiderivative size = 259, normalized size of antiderivative = 1.82

method result size
meijerg \(-\frac {\left (d x \right )^{m} x^{-m} \left (-a \right )^{-\frac {1}{2}-\frac {m}{2}} \left (\frac {2 x^{1+m} \left (-a \right )^{\frac {3}{2}+\frac {m}{2}} \left (-48-16 m \right )}{\left (3+m \right ) \left (1+m \right )^{5} a}-\frac {2 x^{1+m} \left (-a \right )^{\frac {3}{2}+\frac {m}{2}} \left (-24-8 m \right ) \ln \left (-a \,x^{2}+1\right )}{\left (3+m \right ) \left (1+m \right )^{4} a}+\frac {2 x^{1+m} \left (-a \right )^{\frac {3}{2}+\frac {m}{2}} \left (12+4 m \right ) \operatorname {polylog}\left (2, a \,x^{2}\right )}{\left (3+m \right ) \left (1+m \right )^{3} a}+\frac {2 x^{1+m} \left (-a \right )^{\frac {3}{2}+\frac {m}{2}} \left (-6-2 m \right ) \operatorname {polylog}\left (3, a \,x^{2}\right )}{\left (3+m \right ) \left (1+m \right )^{2} a}+\frac {2 x^{1+m} \left (-a \right )^{\frac {3}{2}+\frac {m}{2}} \operatorname {polylog}\left (4, a \,x^{2}\right )}{\left (1+m \right ) a}+\frac {2 x^{1+m} \left (-a \right )^{\frac {3}{2}+\frac {m}{2}} \left (24+8 m \right ) \operatorname {LerchPhi}\left (a \,x^{2}, 1, \frac {1}{2}+\frac {m}{2}\right )}{\left (3+m \right ) \left (1+m \right )^{4} a}\right )}{2}\) \(259\)

input
int((d*x)^m*polylog(4,a*x^2),x,method=_RETURNVERBOSE)
 
output
-1/2*(d*x)^m*x^(-m)*(-a)^(-1/2-1/2*m)*(2/(3+m)*x^(1+m)*(-a)^(3/2+1/2*m)*(- 
48-16*m)/(1+m)^5/a-2/(3+m)*x^(1+m)*(-a)^(3/2+1/2*m)*(-24-8*m)/(1+m)^4*ln(- 
a*x^2+1)/a+2/(3+m)*x^(1+m)*(-a)^(3/2+1/2*m)*(12+4*m)/(1+m)^3/a*polylog(2,a 
*x^2)+2/(3+m)*x^(1+m)*(-a)^(3/2+1/2*m)*(-6-2*m)/(1+m)^2/a*polylog(3,a*x^2) 
+2*x^(1+m)*(-a)^(3/2+1/2*m)/(1+m)/a*polylog(4,a*x^2)+2/(3+m)*x^(1+m)*(-a)^ 
(3/2+1/2*m)*(24+8*m)/(1+m)^4/a*LerchPhi(a*x^2,1,1/2+1/2*m))
 
3.2.7.5 Fricas [F]

\[ \int (d x)^m \operatorname {PolyLog}\left (4,a x^2\right ) \, dx=\int { \left (d x\right )^{m} {\rm Li}_{4}(a x^{2}) \,d x } \]

input
integrate((d*x)^m*polylog(4,a*x^2),x, algorithm="fricas")
 
output
integral((d*x)^m*polylog(4, a*x^2), x)
 
3.2.7.6 Sympy [F]

\[ \int (d x)^m \operatorname {PolyLog}\left (4,a x^2\right ) \, dx=\int \left (d x\right )^{m} \operatorname {Li}_{4}\left (a x^{2}\right )\, dx \]

input
integrate((d*x)**m*polylog(4,a*x**2),x)
 
output
Integral((d*x)**m*polylog(4, a*x**2), x)
 
3.2.7.7 Maxima [F]

\[ \int (d x)^m \operatorname {PolyLog}\left (4,a x^2\right ) \, dx=\int { \left (d x\right )^{m} {\rm Li}_{4}(a x^{2}) \,d x } \]

input
integrate((d*x)^m*polylog(4,a*x^2),x, algorithm="maxima")
 
output
-16*a*d^m*integrate(-x^2*x^m/(m^4 + 4*m^3 - (a*m^4 + 4*a*m^3 + 6*a*m^2 + 4 
*a*m + a)*x^2 + 6*m^2 + 4*m + 1), x) + (4*(d^m*m + d^m)*x*x^m*dilog(a*x^2) 
 + 8*d^m*x*x^m*log(-a*x^2 + 1) + (d^m*m^3 + 3*d^m*m^2 + 3*d^m*m + d^m)*x*x 
^m*polylog(4, a*x^2) - 2*(d^m*m^2 + 2*d^m*m + d^m)*x*x^m*polylog(3, a*x^2) 
)/(m^4 + 4*m^3 + 6*m^2 + 4*m + 1)
 
3.2.7.8 Giac [F]

\[ \int (d x)^m \operatorname {PolyLog}\left (4,a x^2\right ) \, dx=\int { \left (d x\right )^{m} {\rm Li}_{4}(a x^{2}) \,d x } \]

input
integrate((d*x)^m*polylog(4,a*x^2),x, algorithm="giac")
 
output
integrate((d*x)^m*polylog(4, a*x^2), x)
 
3.2.7.9 Mupad [F(-1)]

Timed out. \[ \int (d x)^m \operatorname {PolyLog}\left (4,a x^2\right ) \, dx=\int \mathrm {polylog}\left (4,a\,x^2\right )\,{\left (d\,x\right )}^m \,d x \]

input
int(polylog(4, a*x^2)*(d*x)^m,x)
 
output
int(polylog(4, a*x^2)*(d*x)^m, x)