\(\int \sqrt {a+b x^3} (A+B x^3) \, dx\) [165]

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

Optimal result

Integrand size = 19, antiderivative size = 268 \[ \int \sqrt {a+b x^3} \left (A+B x^3\right ) \, dx=\frac {2 (11 A b-2 a B) x \sqrt {a+b x^3}}{55 b}+\frac {2 B x \left (a+b x^3\right )^{3/2}}{11 b}+\frac {2\ 3^{3/4} \sqrt {2+\sqrt {3}} a (11 A b-2 a B) \left (\sqrt [3]{a}+\sqrt [3]{b} x\right ) \sqrt {\frac {a^{2/3}-\sqrt [3]{a} \sqrt [3]{b} x+b^{2/3} x^2}{\left (\left (1+\sqrt {3}\right ) \sqrt [3]{a}+\sqrt [3]{b} x\right )^2}} \operatorname {EllipticF}\left (\arcsin \left (\frac {\left (1-\sqrt {3}\right ) \sqrt [3]{a}+\sqrt [3]{b} x}{\left (1+\sqrt {3}\right ) \sqrt [3]{a}+\sqrt [3]{b} x}\right ),-7-4 \sqrt {3}\right )}{55 b^{4/3} \sqrt {\frac {\sqrt [3]{a} \left (\sqrt [3]{a}+\sqrt [3]{b} x\right )}{\left (\left (1+\sqrt {3}\right ) \sqrt [3]{a}+\sqrt [3]{b} x\right )^2}} \sqrt {a+b x^3}} \] Output:

2/55*(11*A*b-2*B*a)*x*(b*x^3+a)^(1/2)/b+2/11*B*x*(b*x^3+a)^(3/2)/b+2/55*3^ 
(3/4)*(1/2*6^(1/2)+1/2*2^(1/2))*a*(11*A*b-2*B*a)*(a^(1/3)+b^(1/3)*x)*((a^( 
2/3)-a^(1/3)*b^(1/3)*x+b^(2/3)*x^2)/((1+3^(1/2))*a^(1/3)+b^(1/3)*x)^2)^(1/ 
2)*EllipticF(((1-3^(1/2))*a^(1/3)+b^(1/3)*x)/((1+3^(1/2))*a^(1/3)+b^(1/3)* 
x),I*3^(1/2)+2*I)/b^(4/3)/(a^(1/3)*(a^(1/3)+b^(1/3)*x)/((1+3^(1/2))*a^(1/3 
)+b^(1/3)*x)^2)^(1/2)/(b*x^3+a)^(1/2)
 

Mathematica [C] (verified)

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

Time = 5.76 (sec) , antiderivative size = 75, normalized size of antiderivative = 0.28 \[ \int \sqrt {a+b x^3} \left (A+B x^3\right ) \, dx=\frac {2 x \sqrt {a+b x^3} \left (B \left (a+b x^3\right )+\frac {(11 A b-2 a B) \operatorname {Hypergeometric2F1}\left (-\frac {1}{2},\frac {1}{3},\frac {4}{3},-\frac {b x^3}{a}\right )}{2 \sqrt {1+\frac {b x^3}{a}}}\right )}{11 b} \] Input:

Integrate[Sqrt[a + b*x^3]*(A + B*x^3),x]
 

Output:

(2*x*Sqrt[a + b*x^3]*(B*(a + b*x^3) + ((11*A*b - 2*a*B)*Hypergeometric2F1[ 
-1/2, 1/3, 4/3, -((b*x^3)/a)])/(2*Sqrt[1 + (b*x^3)/a])))/(11*b)
 

Rubi [A] (verified)

Time = 0.49 (sec) , antiderivative size = 264, normalized size of antiderivative = 0.99, number of steps used = 3, number of rules used = 3, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.158, Rules used = {913, 748, 759}

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 {a+b x^3} \left (A+B x^3\right ) \, dx\)

\(\Big \downarrow \) 913

\(\displaystyle \frac {(11 A b-2 a B) \int \sqrt {b x^3+a}dx}{11 b}+\frac {2 B x \left (a+b x^3\right )^{3/2}}{11 b}\)

\(\Big \downarrow \) 748

\(\displaystyle \frac {(11 A b-2 a B) \left (\frac {3}{5} a \int \frac {1}{\sqrt {b x^3+a}}dx+\frac {2}{5} x \sqrt {a+b x^3}\right )}{11 b}+\frac {2 B x \left (a+b x^3\right )^{3/2}}{11 b}\)

\(\Big \downarrow \) 759

\(\displaystyle \frac {(11 A b-2 a B) \left (\frac {2\ 3^{3/4} \sqrt {2+\sqrt {3}} a \left (\sqrt [3]{a}+\sqrt [3]{b} x\right ) \sqrt {\frac {a^{2/3}-\sqrt [3]{a} \sqrt [3]{b} x+b^{2/3} x^2}{\left (\left (1+\sqrt {3}\right ) \sqrt [3]{a}+\sqrt [3]{b} x\right )^2}} \operatorname {EllipticF}\left (\arcsin \left (\frac {\sqrt [3]{b} x+\left (1-\sqrt {3}\right ) \sqrt [3]{a}}{\sqrt [3]{b} x+\left (1+\sqrt {3}\right ) \sqrt [3]{a}}\right ),-7-4 \sqrt {3}\right )}{5 \sqrt [3]{b} \sqrt {\frac {\sqrt [3]{a} \left (\sqrt [3]{a}+\sqrt [3]{b} x\right )}{\left (\left (1+\sqrt {3}\right ) \sqrt [3]{a}+\sqrt [3]{b} x\right )^2}} \sqrt {a+b x^3}}+\frac {2}{5} x \sqrt {a+b x^3}\right )}{11 b}+\frac {2 B x \left (a+b x^3\right )^{3/2}}{11 b}\)

Input:

Int[Sqrt[a + b*x^3]*(A + B*x^3),x]
 

Output:

(2*B*x*(a + b*x^3)^(3/2))/(11*b) + ((11*A*b - 2*a*B)*((2*x*Sqrt[a + b*x^3] 
)/5 + (2*3^(3/4)*Sqrt[2 + Sqrt[3]]*a*(a^(1/3) + b^(1/3)*x)*Sqrt[(a^(2/3) - 
 a^(1/3)*b^(1/3)*x + b^(2/3)*x^2)/((1 + Sqrt[3])*a^(1/3) + b^(1/3)*x)^2]*E 
llipticF[ArcSin[((1 - Sqrt[3])*a^(1/3) + b^(1/3)*x)/((1 + Sqrt[3])*a^(1/3) 
 + b^(1/3)*x)], -7 - 4*Sqrt[3]])/(5*b^(1/3)*Sqrt[(a^(1/3)*(a^(1/3) + b^(1/ 
3)*x))/((1 + Sqrt[3])*a^(1/3) + b^(1/3)*x)^2]*Sqrt[a + b*x^3])))/(11*b)
 

Defintions of rubi rules used

rule 748
Int[((a_) + (b_.)*(x_)^(n_))^(p_), x_Symbol] :> Simp[x*((a + b*x^n)^p/(n*p 
+ 1)), x] + Simp[a*n*(p/(n*p + 1))   Int[(a + b*x^n)^(p - 1), x], x] /; Fre 
eQ[{a, b}, x] && IGtQ[n, 0] && GtQ[p, 0] && (IntegerQ[2*p] || LtQ[Denominat 
or[p + 1/n], Denominator[p]])
 

rule 759
Int[1/Sqrt[(a_) + (b_.)*(x_)^3], x_Symbol] :> With[{r = Numer[Rt[b/a, 3]], 
s = Denom[Rt[b/a, 3]]}, Simp[2*Sqrt[2 + Sqrt[3]]*(s + r*x)*(Sqrt[(s^2 - r*s 
*x + r^2*x^2)/((1 + Sqrt[3])*s + r*x)^2]/(3^(1/4)*r*Sqrt[a + b*x^3]*Sqrt[s* 
((s + r*x)/((1 + Sqrt[3])*s + r*x)^2)]))*EllipticF[ArcSin[((1 - Sqrt[3])*s 
+ r*x)/((1 + Sqrt[3])*s + r*x)], -7 - 4*Sqrt[3]], x]] /; FreeQ[{a, b}, x] & 
& PosQ[a]
 

rule 913
Int[((a_) + (b_.)*(x_)^(n_))^(p_)*((c_) + (d_.)*(x_)^(n_)), x_Symbol] :> Si 
mp[d*x*((a + b*x^n)^(p + 1)/(b*(n*(p + 1) + 1))), x] - Simp[(a*d - b*c*(n*( 
p + 1) + 1))/(b*(n*(p + 1) + 1))   Int[(a + b*x^n)^p, x], x] /; FreeQ[{a, b 
, c, d, n, p}, x] && NeQ[b*c - a*d, 0] && NeQ[n*(p + 1) + 1, 0]
 
Maple [A] (verified)

Time = 0.87 (sec) , antiderivative size = 325, normalized size of antiderivative = 1.21

method result size
risch \(\frac {2 x \left (5 b B \,x^{3}+11 A b +3 B a \right ) \sqrt {b \,x^{3}+a}}{55 b}-\frac {2 i a \left (11 A b -2 B a \right ) \sqrt {3}\, \left (-a \,b^{2}\right )^{\frac {1}{3}} \sqrt {\frac {i \left (x +\frac {\left (-a \,b^{2}\right )^{\frac {1}{3}}}{2 b}-\frac {i \sqrt {3}\, \left (-a \,b^{2}\right )^{\frac {1}{3}}}{2 b}\right ) \sqrt {3}\, b}{\left (-a \,b^{2}\right )^{\frac {1}{3}}}}\, \sqrt {\frac {x -\frac {\left (-a \,b^{2}\right )^{\frac {1}{3}}}{b}}{-\frac {3 \left (-a \,b^{2}\right )^{\frac {1}{3}}}{2 b}+\frac {i \sqrt {3}\, \left (-a \,b^{2}\right )^{\frac {1}{3}}}{2 b}}}\, \sqrt {-\frac {i \left (x +\frac {\left (-a \,b^{2}\right )^{\frac {1}{3}}}{2 b}+\frac {i \sqrt {3}\, \left (-a \,b^{2}\right )^{\frac {1}{3}}}{2 b}\right ) \sqrt {3}\, b}{\left (-a \,b^{2}\right )^{\frac {1}{3}}}}\, \operatorname {EllipticF}\left (\frac {\sqrt {3}\, \sqrt {\frac {i \left (x +\frac {\left (-a \,b^{2}\right )^{\frac {1}{3}}}{2 b}-\frac {i \sqrt {3}\, \left (-a \,b^{2}\right )^{\frac {1}{3}}}{2 b}\right ) \sqrt {3}\, b}{\left (-a \,b^{2}\right )^{\frac {1}{3}}}}}{3}, \sqrt {\frac {i \sqrt {3}\, \left (-a \,b^{2}\right )^{\frac {1}{3}}}{b \left (-\frac {3 \left (-a \,b^{2}\right )^{\frac {1}{3}}}{2 b}+\frac {i \sqrt {3}\, \left (-a \,b^{2}\right )^{\frac {1}{3}}}{2 b}\right )}}\right )}{55 b^{2} \sqrt {b \,x^{3}+a}}\) \(325\)
elliptic \(\frac {2 B \,x^{4} \sqrt {b \,x^{3}+a}}{11}+\frac {2 \left (A b +\frac {3 B a}{11}\right ) x \sqrt {b \,x^{3}+a}}{5 b}-\frac {2 i \left (A a -\frac {2 a \left (A b +\frac {3 B a}{11}\right )}{5 b}\right ) \sqrt {3}\, \left (-a \,b^{2}\right )^{\frac {1}{3}} \sqrt {\frac {i \left (x +\frac {\left (-a \,b^{2}\right )^{\frac {1}{3}}}{2 b}-\frac {i \sqrt {3}\, \left (-a \,b^{2}\right )^{\frac {1}{3}}}{2 b}\right ) \sqrt {3}\, b}{\left (-a \,b^{2}\right )^{\frac {1}{3}}}}\, \sqrt {\frac {x -\frac {\left (-a \,b^{2}\right )^{\frac {1}{3}}}{b}}{-\frac {3 \left (-a \,b^{2}\right )^{\frac {1}{3}}}{2 b}+\frac {i \sqrt {3}\, \left (-a \,b^{2}\right )^{\frac {1}{3}}}{2 b}}}\, \sqrt {-\frac {i \left (x +\frac {\left (-a \,b^{2}\right )^{\frac {1}{3}}}{2 b}+\frac {i \sqrt {3}\, \left (-a \,b^{2}\right )^{\frac {1}{3}}}{2 b}\right ) \sqrt {3}\, b}{\left (-a \,b^{2}\right )^{\frac {1}{3}}}}\, \operatorname {EllipticF}\left (\frac {\sqrt {3}\, \sqrt {\frac {i \left (x +\frac {\left (-a \,b^{2}\right )^{\frac {1}{3}}}{2 b}-\frac {i \sqrt {3}\, \left (-a \,b^{2}\right )^{\frac {1}{3}}}{2 b}\right ) \sqrt {3}\, b}{\left (-a \,b^{2}\right )^{\frac {1}{3}}}}}{3}, \sqrt {\frac {i \sqrt {3}\, \left (-a \,b^{2}\right )^{\frac {1}{3}}}{b \left (-\frac {3 \left (-a \,b^{2}\right )^{\frac {1}{3}}}{2 b}+\frac {i \sqrt {3}\, \left (-a \,b^{2}\right )^{\frac {1}{3}}}{2 b}\right )}}\right )}{3 b \sqrt {b \,x^{3}+a}}\) \(340\)
default \(A \left (\frac {2 x \sqrt {b \,x^{3}+a}}{5}-\frac {2 i a \sqrt {3}\, \left (-a \,b^{2}\right )^{\frac {1}{3}} \sqrt {\frac {i \left (x +\frac {\left (-a \,b^{2}\right )^{\frac {1}{3}}}{2 b}-\frac {i \sqrt {3}\, \left (-a \,b^{2}\right )^{\frac {1}{3}}}{2 b}\right ) \sqrt {3}\, b}{\left (-a \,b^{2}\right )^{\frac {1}{3}}}}\, \sqrt {\frac {x -\frac {\left (-a \,b^{2}\right )^{\frac {1}{3}}}{b}}{-\frac {3 \left (-a \,b^{2}\right )^{\frac {1}{3}}}{2 b}+\frac {i \sqrt {3}\, \left (-a \,b^{2}\right )^{\frac {1}{3}}}{2 b}}}\, \sqrt {-\frac {i \left (x +\frac {\left (-a \,b^{2}\right )^{\frac {1}{3}}}{2 b}+\frac {i \sqrt {3}\, \left (-a \,b^{2}\right )^{\frac {1}{3}}}{2 b}\right ) \sqrt {3}\, b}{\left (-a \,b^{2}\right )^{\frac {1}{3}}}}\, \operatorname {EllipticF}\left (\frac {\sqrt {3}\, \sqrt {\frac {i \left (x +\frac {\left (-a \,b^{2}\right )^{\frac {1}{3}}}{2 b}-\frac {i \sqrt {3}\, \left (-a \,b^{2}\right )^{\frac {1}{3}}}{2 b}\right ) \sqrt {3}\, b}{\left (-a \,b^{2}\right )^{\frac {1}{3}}}}}{3}, \sqrt {\frac {i \sqrt {3}\, \left (-a \,b^{2}\right )^{\frac {1}{3}}}{b \left (-\frac {3 \left (-a \,b^{2}\right )^{\frac {1}{3}}}{2 b}+\frac {i \sqrt {3}\, \left (-a \,b^{2}\right )^{\frac {1}{3}}}{2 b}\right )}}\right )}{5 b \sqrt {b \,x^{3}+a}}\right )+B \left (\frac {2 x^{4} \sqrt {b \,x^{3}+a}}{11}+\frac {6 a x \sqrt {b \,x^{3}+a}}{55 b}+\frac {4 i a^{2} \sqrt {3}\, \left (-a \,b^{2}\right )^{\frac {1}{3}} \sqrt {\frac {i \left (x +\frac {\left (-a \,b^{2}\right )^{\frac {1}{3}}}{2 b}-\frac {i \sqrt {3}\, \left (-a \,b^{2}\right )^{\frac {1}{3}}}{2 b}\right ) \sqrt {3}\, b}{\left (-a \,b^{2}\right )^{\frac {1}{3}}}}\, \sqrt {\frac {x -\frac {\left (-a \,b^{2}\right )^{\frac {1}{3}}}{b}}{-\frac {3 \left (-a \,b^{2}\right )^{\frac {1}{3}}}{2 b}+\frac {i \sqrt {3}\, \left (-a \,b^{2}\right )^{\frac {1}{3}}}{2 b}}}\, \sqrt {-\frac {i \left (x +\frac {\left (-a \,b^{2}\right )^{\frac {1}{3}}}{2 b}+\frac {i \sqrt {3}\, \left (-a \,b^{2}\right )^{\frac {1}{3}}}{2 b}\right ) \sqrt {3}\, b}{\left (-a \,b^{2}\right )^{\frac {1}{3}}}}\, \operatorname {EllipticF}\left (\frac {\sqrt {3}\, \sqrt {\frac {i \left (x +\frac {\left (-a \,b^{2}\right )^{\frac {1}{3}}}{2 b}-\frac {i \sqrt {3}\, \left (-a \,b^{2}\right )^{\frac {1}{3}}}{2 b}\right ) \sqrt {3}\, b}{\left (-a \,b^{2}\right )^{\frac {1}{3}}}}}{3}, \sqrt {\frac {i \sqrt {3}\, \left (-a \,b^{2}\right )^{\frac {1}{3}}}{b \left (-\frac {3 \left (-a \,b^{2}\right )^{\frac {1}{3}}}{2 b}+\frac {i \sqrt {3}\, \left (-a \,b^{2}\right )^{\frac {1}{3}}}{2 b}\right )}}\right )}{55 b^{2} \sqrt {b \,x^{3}+a}}\right )\) \(618\)

Input:

int((b*x^3+a)^(1/2)*(B*x^3+A),x,method=_RETURNVERBOSE)
 

Output:

2/55*x*(5*B*b*x^3+11*A*b+3*B*a)/b*(b*x^3+a)^(1/2)-2/55*I*a*(11*A*b-2*B*a)/ 
b^2*3^(1/2)*(-a*b^2)^(1/3)*(I*(x+1/2/b*(-a*b^2)^(1/3)-1/2*I*3^(1/2)/b*(-a* 
b^2)^(1/3))*3^(1/2)*b/(-a*b^2)^(1/3))^(1/2)*((x-1/b*(-a*b^2)^(1/3))/(-3/2/ 
b*(-a*b^2)^(1/3)+1/2*I*3^(1/2)/b*(-a*b^2)^(1/3)))^(1/2)*(-I*(x+1/2/b*(-a*b 
^2)^(1/3)+1/2*I*3^(1/2)/b*(-a*b^2)^(1/3))*3^(1/2)*b/(-a*b^2)^(1/3))^(1/2)/ 
(b*x^3+a)^(1/2)*EllipticF(1/3*3^(1/2)*(I*(x+1/2/b*(-a*b^2)^(1/3)-1/2*I*3^( 
1/2)/b*(-a*b^2)^(1/3))*3^(1/2)*b/(-a*b^2)^(1/3))^(1/2),(I*3^(1/2)/b*(-a*b^ 
2)^(1/3)/(-3/2/b*(-a*b^2)^(1/3)+1/2*I*3^(1/2)/b*(-a*b^2)^(1/3)))^(1/2))
 

Fricas [A] (verification not implemented)

Time = 0.09 (sec) , antiderivative size = 67, normalized size of antiderivative = 0.25 \[ \int \sqrt {a+b x^3} \left (A+B x^3\right ) \, dx=-\frac {2 \, {\left (3 \, {\left (2 \, B a^{2} - 11 \, A a b\right )} \sqrt {b} {\rm weierstrassPInverse}\left (0, -\frac {4 \, a}{b}, x\right ) - {\left (5 \, B b^{2} x^{4} + {\left (3 \, B a b + 11 \, A b^{2}\right )} x\right )} \sqrt {b x^{3} + a}\right )}}{55 \, b^{2}} \] Input:

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

Output:

-2/55*(3*(2*B*a^2 - 11*A*a*b)*sqrt(b)*weierstrassPInverse(0, -4*a/b, x) - 
(5*B*b^2*x^4 + (3*B*a*b + 11*A*b^2)*x)*sqrt(b*x^3 + a))/b^2
 

Sympy [A] (verification not implemented)

Time = 1.12 (sec) , antiderivative size = 82, normalized size of antiderivative = 0.31 \[ \int \sqrt {a+b x^3} \left (A+B x^3\right ) \, dx=\frac {A \sqrt {a} x \Gamma \left (\frac {1}{3}\right ) {{}_{2}F_{1}\left (\begin {matrix} - \frac {1}{2}, \frac {1}{3} \\ \frac {4}{3} \end {matrix}\middle | {\frac {b x^{3} e^{i \pi }}{a}} \right )}}{3 \Gamma \left (\frac {4}{3}\right )} + \frac {B \sqrt {a} x^{4} \Gamma \left (\frac {4}{3}\right ) {{}_{2}F_{1}\left (\begin {matrix} - \frac {1}{2}, \frac {4}{3} \\ \frac {7}{3} \end {matrix}\middle | {\frac {b x^{3} e^{i \pi }}{a}} \right )}}{3 \Gamma \left (\frac {7}{3}\right )} \] Input:

integrate((b*x**3+a)**(1/2)*(B*x**3+A),x)
 

Output:

A*sqrt(a)*x*gamma(1/3)*hyper((-1/2, 1/3), (4/3,), b*x**3*exp_polar(I*pi)/a 
)/(3*gamma(4/3)) + B*sqrt(a)*x**4*gamma(4/3)*hyper((-1/2, 4/3), (7/3,), b* 
x**3*exp_polar(I*pi)/a)/(3*gamma(7/3))
 

Maxima [F]

\[ \int \sqrt {a+b x^3} \left (A+B x^3\right ) \, dx=\int { {\left (B x^{3} + A\right )} \sqrt {b x^{3} + a} \,d x } \] Input:

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

Output:

integrate((B*x^3 + A)*sqrt(b*x^3 + a), x)
 

Giac [F]

\[ \int \sqrt {a+b x^3} \left (A+B x^3\right ) \, dx=\int { {\left (B x^{3} + A\right )} \sqrt {b x^{3} + a} \,d x } \] Input:

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

Output:

integrate((B*x^3 + A)*sqrt(b*x^3 + a), x)
 

Mupad [F(-1)]

Timed out. \[ \int \sqrt {a+b x^3} \left (A+B x^3\right ) \, dx=\int \left (B\,x^3+A\right )\,\sqrt {b\,x^3+a} \,d x \] Input:

int((A + B*x^3)*(a + b*x^3)^(1/2),x)
 

Output:

int((A + B*x^3)*(a + b*x^3)^(1/2), x)
 

Reduce [F]

\[ \int \sqrt {a+b x^3} \left (A+B x^3\right ) \, dx=\frac {28 \sqrt {b \,x^{3}+a}\, a x}{55}+\frac {2 \sqrt {b \,x^{3}+a}\, b \,x^{4}}{11}+\frac {27 \left (\int \frac {\sqrt {b \,x^{3}+a}}{b \,x^{3}+a}d x \right ) a^{2}}{55} \] Input:

int((b*x^3+a)^(1/2)*(B*x^3+A),x)
                                                                                    
                                                                                    
 

Output:

(28*sqrt(a + b*x**3)*a*x + 10*sqrt(a + b*x**3)*b*x**4 + 27*int(sqrt(a + b* 
x**3)/(a + b*x**3),x)*a**2)/55