\(\int \sqrt [3]{x} \sqrt {a+b x} \sqrt {c+d x} \, dx\) [476]

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

Optimal result

Integrand size = 24, antiderivative size = 78 \[ \int \sqrt [3]{x} \sqrt {a+b x} \sqrt {c+d x} \, dx=\frac {3 x^{4/3} \sqrt {a+b x} \sqrt {c+d x} \operatorname {AppellF1}\left (\frac {4}{3},-\frac {1}{2},-\frac {1}{2},\frac {7}{3},-\frac {b x}{a},-\frac {d x}{c}\right )}{4 \sqrt {1+\frac {b x}{a}} \sqrt {1+\frac {d x}{c}}} \] Output:

3/4*x^(4/3)*(b*x+a)^(1/2)*(d*x+c)^(1/2)*AppellF1(4/3,-1/2,-1/2,7/3,-b*x/a, 
-d*x/c)/(1+b*x/a)^(1/2)/(1+d*x/c)^(1/2)
 

Mathematica [B] (verified)

Leaf count is larger than twice the leaf count of optimal. \(199\) vs. \(2(78)=156\).

Time = 2.84 (sec) , antiderivative size = 199, normalized size of antiderivative = 2.55 \[ \int \sqrt [3]{x} \sqrt {a+b x} \sqrt {c+d x} \, dx=\frac {3 \sqrt [3]{x} \left (8 (a+b x) (c+d x) (3 b c+3 a d+8 b d x)-24 a c (b c+a d) \sqrt {1+\frac {b x}{a}} \sqrt {1+\frac {d x}{c}} \operatorname {AppellF1}\left (\frac {1}{3},\frac {1}{2},\frac {1}{2},\frac {4}{3},-\frac {b x}{a},-\frac {d x}{c}\right )-3 \left (5 b^2 c^2-6 a b c d+5 a^2 d^2\right ) x \sqrt {1+\frac {b x}{a}} \sqrt {1+\frac {d x}{c}} \operatorname {AppellF1}\left (\frac {4}{3},\frac {1}{2},\frac {1}{2},\frac {7}{3},-\frac {b x}{a},-\frac {d x}{c}\right )\right )}{448 b d \sqrt {a+b x} \sqrt {c+d x}} \] Input:

Integrate[x^(1/3)*Sqrt[a + b*x]*Sqrt[c + d*x],x]
 

Output:

(3*x^(1/3)*(8*(a + b*x)*(c + d*x)*(3*b*c + 3*a*d + 8*b*d*x) - 24*a*c*(b*c 
+ a*d)*Sqrt[1 + (b*x)/a]*Sqrt[1 + (d*x)/c]*AppellF1[1/3, 1/2, 1/2, 4/3, -( 
(b*x)/a), -((d*x)/c)] - 3*(5*b^2*c^2 - 6*a*b*c*d + 5*a^2*d^2)*x*Sqrt[1 + ( 
b*x)/a]*Sqrt[1 + (d*x)/c]*AppellF1[4/3, 1/2, 1/2, 7/3, -((b*x)/a), -((d*x) 
/c)]))/(448*b*d*Sqrt[a + b*x]*Sqrt[c + d*x])
 

Rubi [A] (verified)

Time = 0.18 (sec) , antiderivative size = 78, normalized size of antiderivative = 1.00, number of steps used = 3, number of rules used = 3, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.125, Rules used = {152, 152, 150}

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 \sqrt [3]{x} \sqrt {a+b x} \sqrt {c+d x} \, dx\)

\(\Big \downarrow \) 152

\(\displaystyle \frac {\sqrt {a+b x} \int \sqrt [3]{x} \sqrt {\frac {b x}{a}+1} \sqrt {c+d x}dx}{\sqrt {\frac {b x}{a}+1}}\)

\(\Big \downarrow \) 152

\(\displaystyle \frac {\sqrt {a+b x} \sqrt {c+d x} \int \sqrt [3]{x} \sqrt {\frac {b x}{a}+1} \sqrt {\frac {d x}{c}+1}dx}{\sqrt {\frac {b x}{a}+1} \sqrt {\frac {d x}{c}+1}}\)

\(\Big \downarrow \) 150

\(\displaystyle \frac {3 x^{4/3} \sqrt {a+b x} \sqrt {c+d x} \operatorname {AppellF1}\left (\frac {4}{3},-\frac {1}{2},-\frac {1}{2},\frac {7}{3},-\frac {b x}{a},-\frac {d x}{c}\right )}{4 \sqrt {\frac {b x}{a}+1} \sqrt {\frac {d x}{c}+1}}\)

Input:

Int[x^(1/3)*Sqrt[a + b*x]*Sqrt[c + d*x],x]
 

Output:

(3*x^(4/3)*Sqrt[a + b*x]*Sqrt[c + d*x]*AppellF1[4/3, -1/2, -1/2, 7/3, -((b 
*x)/a), -((d*x)/c)])/(4*Sqrt[1 + (b*x)/a]*Sqrt[1 + (d*x)/c])
 

Defintions of rubi rules used

rule 150
Int[((b_.)*(x_))^(m_)*((c_) + (d_.)*(x_))^(n_)*((e_) + (f_.)*(x_))^(p_), x_ 
] :> Simp[c^n*e^p*((b*x)^(m + 1)/(b*(m + 1)))*AppellF1[m + 1, -n, -p, m + 2 
, (-d)*(x/c), (-f)*(x/e)], x] /; FreeQ[{b, c, d, e, f, m, n, p}, x] &&  !In 
tegerQ[m] &&  !IntegerQ[n] && GtQ[c, 0] && (IntegerQ[p] || GtQ[e, 0])
 

rule 152
Int[((b_.)*(x_))^(m_)*((c_) + (d_.)*(x_))^(n_)*((e_) + (f_.)*(x_))^(p_), x_ 
] :> Simp[c^IntPart[n]*((c + d*x)^FracPart[n]/(1 + d*(x/c))^FracPart[n]) 
Int[(b*x)^m*(1 + d*(x/c))^n*(e + f*x)^p, x], x] /; FreeQ[{b, c, d, e, f, m, 
 n, p}, x] &&  !IntegerQ[m] &&  !IntegerQ[n] &&  !GtQ[c, 0]
 
Maple [F]

\[\int x^{\frac {1}{3}} \sqrt {b x +a}\, \sqrt {x d +c}d x\]

Input:

int(x^(1/3)*(b*x+a)^(1/2)*(d*x+c)^(1/2),x)
 

Output:

int(x^(1/3)*(b*x+a)^(1/2)*(d*x+c)^(1/2),x)
 

Fricas [F]

\[ \int \sqrt [3]{x} \sqrt {a+b x} \sqrt {c+d x} \, dx=\int { \sqrt {b x + a} \sqrt {d x + c} x^{\frac {1}{3}} \,d x } \] Input:

integrate(x^(1/3)*(b*x+a)^(1/2)*(d*x+c)^(1/2),x, algorithm="fricas")
 

Output:

integral(sqrt(b*x + a)*sqrt(d*x + c)*x^(1/3), x)
 

Sympy [F]

\[ \int \sqrt [3]{x} \sqrt {a+b x} \sqrt {c+d x} \, dx=\int \sqrt [3]{x} \sqrt {a + b x} \sqrt {c + d x}\, dx \] Input:

integrate(x**(1/3)*(b*x+a)**(1/2)*(d*x+c)**(1/2),x)
 

Output:

Integral(x**(1/3)*sqrt(a + b*x)*sqrt(c + d*x), x)
 

Maxima [F]

\[ \int \sqrt [3]{x} \sqrt {a+b x} \sqrt {c+d x} \, dx=\int { \sqrt {b x + a} \sqrt {d x + c} x^{\frac {1}{3}} \,d x } \] Input:

integrate(x^(1/3)*(b*x+a)^(1/2)*(d*x+c)^(1/2),x, algorithm="maxima")
 

Output:

integrate(sqrt(b*x + a)*sqrt(d*x + c)*x^(1/3), x)
 

Giac [F]

\[ \int \sqrt [3]{x} \sqrt {a+b x} \sqrt {c+d x} \, dx=\int { \sqrt {b x + a} \sqrt {d x + c} x^{\frac {1}{3}} \,d x } \] Input:

integrate(x^(1/3)*(b*x+a)^(1/2)*(d*x+c)^(1/2),x, algorithm="giac")
 

Output:

integrate(sqrt(b*x + a)*sqrt(d*x + c)*x^(1/3), x)
 

Mupad [F(-1)]

Timed out. \[ \int \sqrt [3]{x} \sqrt {a+b x} \sqrt {c+d x} \, dx=\int x^{1/3}\,\sqrt {a+b\,x}\,\sqrt {c+d\,x} \,d x \] Input:

int(x^(1/3)*(a + b*x)^(1/2)*(c + d*x)^(1/2),x)
 

Output:

int(x^(1/3)*(a + b*x)^(1/2)*(c + d*x)^(1/2), x)
 

Reduce [F]

\[ \int \sqrt [3]{x} \sqrt {a+b x} \sqrt {c+d x} \, dx=\frac {-\frac {45 \sqrt {d x +c}\, \sqrt {b x +a}\, a^{2} d^{2}}{112}+\frac {27 \sqrt {d x +c}\, \sqrt {b x +a}\, a b c d}{56}+\frac {9 \sqrt {d x +c}\, \sqrt {b x +a}\, a b \,d^{2} x}{56}-\frac {45 \sqrt {d x +c}\, \sqrt {b x +a}\, b^{2} c^{2}}{112}+\frac {9 \sqrt {d x +c}\, \sqrt {b x +a}\, b^{2} c d x}{56}+\frac {3 \sqrt {d x +c}\, \sqrt {b x +a}\, b^{2} d^{2} x^{2}}{7}-\frac {15 x^{\frac {2}{3}} \left (\int \frac {\sqrt {d x +c}\, \sqrt {b x +a}}{x^{\frac {5}{3}} a c +x^{\frac {8}{3}} a d +x^{\frac {8}{3}} b c +x^{\frac {11}{3}} b d}d x \right ) a^{3} c \,d^{2}}{56}+\frac {9 x^{\frac {2}{3}} \left (\int \frac {\sqrt {d x +c}\, \sqrt {b x +a}}{x^{\frac {5}{3}} a c +x^{\frac {8}{3}} a d +x^{\frac {8}{3}} b c +x^{\frac {11}{3}} b d}d x \right ) a^{2} b \,c^{2} d}{28}-\frac {15 x^{\frac {2}{3}} \left (\int \frac {\sqrt {d x +c}\, \sqrt {b x +a}}{x^{\frac {5}{3}} a c +x^{\frac {8}{3}} a d +x^{\frac {8}{3}} b c +x^{\frac {11}{3}} b d}d x \right ) a \,b^{2} c^{3}}{56}-\frac {15 x^{\frac {2}{3}} \left (\int \frac {\sqrt {d x +c}\, \sqrt {b x +a}}{x^{\frac {2}{3}} a c +x^{\frac {5}{3}} a d +x^{\frac {5}{3}} b c +x^{\frac {8}{3}} b d}d x \right ) a^{3} d^{3}}{224}-\frac {9 x^{\frac {2}{3}} \left (\int \frac {\sqrt {d x +c}\, \sqrt {b x +a}}{x^{\frac {2}{3}} a c +x^{\frac {5}{3}} a d +x^{\frac {5}{3}} b c +x^{\frac {8}{3}} b d}d x \right ) a^{2} b c \,d^{2}}{224}-\frac {9 x^{\frac {2}{3}} \left (\int \frac {\sqrt {d x +c}\, \sqrt {b x +a}}{x^{\frac {2}{3}} a c +x^{\frac {5}{3}} a d +x^{\frac {5}{3}} b c +x^{\frac {8}{3}} b d}d x \right ) a \,b^{2} c^{2} d}{224}-\frac {15 x^{\frac {2}{3}} \left (\int \frac {\sqrt {d x +c}\, \sqrt {b x +a}}{x^{\frac {2}{3}} a c +x^{\frac {5}{3}} a d +x^{\frac {5}{3}} b c +x^{\frac {8}{3}} b d}d x \right ) b^{3} c^{3}}{224}}{x^{\frac {2}{3}} b^{2} d^{2}} \] Input:

int(x^(1/3)*(b*x+a)^(1/2)*(d*x+c)^(1/2),x)
 

Output:

(3*( - 30*sqrt(c + d*x)*sqrt(a + b*x)*a**2*d**2 + 36*sqrt(c + d*x)*sqrt(a 
+ b*x)*a*b*c*d + 12*sqrt(c + d*x)*sqrt(a + b*x)*a*b*d**2*x - 30*sqrt(c + d 
*x)*sqrt(a + b*x)*b**2*c**2 + 12*sqrt(c + d*x)*sqrt(a + b*x)*b**2*c*d*x + 
32*sqrt(c + d*x)*sqrt(a + b*x)*b**2*d**2*x**2 - 20*x**(2/3)*int((sqrt(c + 
d*x)*sqrt(a + b*x))/(x**(2/3)*a*c*x + x**(2/3)*a*d*x**2 + x**(2/3)*b*c*x** 
2 + x**(2/3)*b*d*x**3),x)*a**3*c*d**2 + 24*x**(2/3)*int((sqrt(c + d*x)*sqr 
t(a + b*x))/(x**(2/3)*a*c*x + x**(2/3)*a*d*x**2 + x**(2/3)*b*c*x**2 + x**( 
2/3)*b*d*x**3),x)*a**2*b*c**2*d - 20*x**(2/3)*int((sqrt(c + d*x)*sqrt(a + 
b*x))/(x**(2/3)*a*c*x + x**(2/3)*a*d*x**2 + x**(2/3)*b*c*x**2 + x**(2/3)*b 
*d*x**3),x)*a*b**2*c**3 - 5*x**(2/3)*int((sqrt(c + d*x)*sqrt(a + b*x))/(x* 
*(2/3)*a*c + x**(2/3)*a*d*x + x**(2/3)*b*c*x + x**(2/3)*b*d*x**2),x)*a**3* 
d**3 - 3*x**(2/3)*int((sqrt(c + d*x)*sqrt(a + b*x))/(x**(2/3)*a*c + x**(2/ 
3)*a*d*x + x**(2/3)*b*c*x + x**(2/3)*b*d*x**2),x)*a**2*b*c*d**2 - 3*x**(2/ 
3)*int((sqrt(c + d*x)*sqrt(a + b*x))/(x**(2/3)*a*c + x**(2/3)*a*d*x + x**( 
2/3)*b*c*x + x**(2/3)*b*d*x**2),x)*a*b**2*c**2*d - 5*x**(2/3)*int((sqrt(c 
+ d*x)*sqrt(a + b*x))/(x**(2/3)*a*c + x**(2/3)*a*d*x + x**(2/3)*b*c*x + x* 
*(2/3)*b*d*x**2),x)*b**3*c**3))/(224*x**(2/3)*b**2*d**2)