\(\int \frac {a c+a d x+b c x^3+b d x^4}{\sqrt {a+b x^3}} \, dx\) [54]

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 = 32, antiderivative size = 525 \[ \int \frac {a c+a d x+b c x^3+b d x^4}{\sqrt {a+b x^3}} \, dx=\frac {6 a d \sqrt {a+b x^3}}{7 b^{2/3} \left (\left (1+\sqrt {3}\right ) \sqrt [3]{a}+\sqrt [3]{b} x\right )}+\frac {2}{35} \left (7 c x+5 d x^2\right ) \sqrt {a+b x^3}-\frac {3 \sqrt [4]{3} \sqrt {2-\sqrt {3}} a^{4/3} d \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}} E\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 )}{7 b^{2/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}}+\frac {2\ 3^{3/4} \sqrt {2+\sqrt {3}} a \left (7 \sqrt [3]{b} c-5 \left (1-\sqrt {3}\right ) \sqrt [3]{a} d\right ) \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 )}{35 b^{2/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:

6/7*a*d*(b*x^3+a)^(1/2)/b^(2/3)/((1+3^(1/2))*a^(1/3)+b^(1/3)*x)+2/35*(5*d* 
x^2+7*c*x)*(b*x^3+a)^(1/2)-3/7*3^(1/4)*(1/2*6^(1/2)-1/2*2^(1/2))*a^(4/3)*d 
*(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)*EllipticE(((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^(2/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)+2/35*3^(3/ 
4)*(1/2*6^(1/2)+1/2*2^(1/2))*a*(7*b^(1/3)*c-5*(1-3^(1/2))*a^(1/3)*d)*(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^(2/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.41 (sec) , antiderivative size = 75, normalized size of antiderivative = 0.14 \[ \int \frac {a c+a d x+b c x^3+b d x^4}{\sqrt {a+b x^3}} \, dx=\frac {x \sqrt {a+b x^3} \left (2 c \operatorname {Hypergeometric2F1}\left (-\frac {1}{2},\frac {1}{3},\frac {4}{3},-\frac {b x^3}{a}\right )+d x \operatorname {Hypergeometric2F1}\left (-\frac {1}{2},\frac {2}{3},\frac {5}{3},-\frac {b x^3}{a}\right )\right )}{2 \sqrt {1+\frac {b x^3}{a}}} \] Input:

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

Output:

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

Rubi [A] (verified)

Time = 1.02 (sec) , antiderivative size = 530, normalized size of antiderivative = 1.01, number of steps used = 6, number of rules used = 6, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.188, Rules used = {2019, 2392, 27, 2417, 759, 2416}

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

\(\Big \downarrow \) 2019

\(\displaystyle \int \sqrt {a+b x^3} (c+d x)dx\)

\(\Big \downarrow \) 2392

\(\displaystyle \frac {3}{2} a \int \frac {2 (7 c+5 d x)}{35 \sqrt {b x^3+a}}dx+\frac {2}{35} \sqrt {a+b x^3} \left (7 c x+5 d x^2\right )\)

\(\Big \downarrow \) 27

\(\displaystyle \frac {3}{35} a \int \frac {7 c+5 d x}{\sqrt {b x^3+a}}dx+\frac {2}{35} \sqrt {a+b x^3} \left (7 c x+5 d x^2\right )\)

\(\Big \downarrow \) 2417

\(\displaystyle \frac {3}{35} a \left (\left (7 c-\frac {5 \left (1-\sqrt {3}\right ) \sqrt [3]{a} d}{\sqrt [3]{b}}\right ) \int \frac {1}{\sqrt {b x^3+a}}dx+\frac {5 d \int \frac {\sqrt [3]{b} x+\left (1-\sqrt {3}\right ) \sqrt [3]{a}}{\sqrt {b x^3+a}}dx}{\sqrt [3]{b}}\right )+\frac {2}{35} \sqrt {a+b x^3} \left (7 c x+5 d x^2\right )\)

\(\Big \downarrow \) 759

\(\displaystyle \frac {3}{35} a \left (\frac {5 d \int \frac {\sqrt [3]{b} x+\left (1-\sqrt {3}\right ) \sqrt [3]{a}}{\sqrt {b x^3+a}}dx}{\sqrt [3]{b}}+\frac {2 \sqrt {2+\sqrt {3}} \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}} \left (7 c-\frac {5 \left (1-\sqrt {3}\right ) \sqrt [3]{a} d}{\sqrt [3]{b}}\right ) \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 )}{\sqrt [4]{3} \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}}\right )+\frac {2}{35} \sqrt {a+b x^3} \left (7 c x+5 d x^2\right )\)

\(\Big \downarrow \) 2416

\(\displaystyle \frac {3}{35} a \left (\frac {2 \sqrt {2+\sqrt {3}} \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}} \left (7 c-\frac {5 \left (1-\sqrt {3}\right ) \sqrt [3]{a} d}{\sqrt [3]{b}}\right ) \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 )}{\sqrt [4]{3} \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 {5 d \left (\frac {2 \sqrt {a+b x^3}}{\sqrt [3]{b} \left (\left (1+\sqrt {3}\right ) \sqrt [3]{a}+\sqrt [3]{b} x\right )}-\frac {\sqrt [4]{3} \sqrt {2-\sqrt {3}} \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}} E\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 )}{\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}}\right )}{\sqrt [3]{b}}\right )+\frac {2}{35} \sqrt {a+b x^3} \left (7 c x+5 d x^2\right )\)

Input:

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

Output:

(2*(7*c*x + 5*d*x^2)*Sqrt[a + b*x^3])/35 + (3*a*((5*d*((2*Sqrt[a + b*x^3]) 
/(b^(1/3)*((1 + Sqrt[3])*a^(1/3) + b^(1/3)*x)) - (3^(1/4)*Sqrt[2 - Sqrt[3] 
]*a^(1/3)*(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]*EllipticE[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]])/(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])))/b^(1/3) + (2*Sqrt[2 + Sqrt[3]]*(7*c - 
(5*(1 - Sqrt[3])*a^(1/3)*d)/b^(1/3))*(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]])/(3^(1/4)*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])))/35
 

Defintions of rubi rules used

rule 27
Int[(a_)*(Fx_), x_Symbol] :> Simp[a   Int[Fx, x], x] /; FreeQ[a, x] &&  !Ma 
tchQ[Fx, (b_)*(Gx_) /; FreeQ[b, x]]
 

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 2019
Int[(u_.)*(Px_)^(p_.)*(Qx_)^(q_.), x_Symbol] :> Int[u*PolynomialQuotient[Px 
, Qx, x]^p*Qx^(p + q), x] /; FreeQ[q, x] && PolyQ[Px, x] && PolyQ[Qx, x] && 
 EqQ[PolynomialRemainder[Px, Qx, x], 0] && IntegerQ[p] && LtQ[p*q, 0]
 

rule 2392
Int[(Pq_)*((a_) + (b_.)*(x_)^(n_.))^(p_), x_Symbol] :> Module[{q = Expon[Pq 
, x], i}, Simp[(a + b*x^n)^p*Sum[Coeff[Pq, x, i]*(x^(i + 1)/(n*p + i + 1)), 
 {i, 0, q}], x] + Simp[a*n*p   Int[(a + b*x^n)^(p - 1)*Sum[Coeff[Pq, x, i]* 
(x^i/(n*p + i + 1)), {i, 0, q}], x], x]] /; FreeQ[{a, b}, x] && PolyQ[Pq, x 
] && IGtQ[(n - 1)/2, 0] && GtQ[p, 0]
 

rule 2416
Int[((c_) + (d_.)*(x_))/Sqrt[(a_) + (b_.)*(x_)^3], x_Symbol] :> With[{r = N 
umer[Simplify[(1 - Sqrt[3])*(d/c)]], s = Denom[Simplify[(1 - Sqrt[3])*(d/c) 
]]}, Simp[2*d*s^3*(Sqrt[a + b*x^3]/(a*r^2*((1 + Sqrt[3])*s + r*x))), x] - S 
imp[3^(1/4)*Sqrt[2 - Sqrt[3]]*d*s*(s + r*x)*(Sqrt[(s^2 - r*s*x + r^2*x^2)/( 
(1 + Sqrt[3])*s + r*x)^2]/(r^2*Sqrt[a + b*x^3]*Sqrt[s*((s + r*x)/((1 + Sqrt 
[3])*s + r*x)^2)]))*EllipticE[ArcSin[((1 - Sqrt[3])*s + r*x)/((1 + Sqrt[3]) 
*s + r*x)], -7 - 4*Sqrt[3]], x]] /; FreeQ[{a, b, c, d}, x] && PosQ[a] && Eq 
Q[b*c^3 - 2*(5 - 3*Sqrt[3])*a*d^3, 0]
 

rule 2417
Int[((c_) + (d_.)*(x_))/Sqrt[(a_) + (b_.)*(x_)^3], x_Symbol] :> With[{r = N 
umer[Rt[b/a, 3]], s = Denom[Rt[b/a, 3]]}, Simp[(c*r - (1 - Sqrt[3])*d*s)/r 
  Int[1/Sqrt[a + b*x^3], x], x] + Simp[d/r   Int[((1 - Sqrt[3])*s + r*x)/Sq 
rt[a + b*x^3], x], x]] /; FreeQ[{a, b, c, d}, x] && PosQ[a] && NeQ[b*c^3 - 
2*(5 - 3*Sqrt[3])*a*d^3, 0]
 
Maple [A] (verified)

Time = 0.63 (sec) , antiderivative size = 744, normalized size of antiderivative = 1.42

method result size
risch \(\frac {2 x \left (5 d x +7 c \right ) \sqrt {b \,x^{3}+a}}{35}+\frac {3 a \left (-\frac {14 i c \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}}-\frac {10 i d \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}}}}\, \left (\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 ) \operatorname {EllipticE}\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 )+\frac {\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 )}{b}\right )}{3 b \sqrt {b \,x^{3}+a}}\right )}{35}\) \(744\)
elliptic \(\frac {2 d \,x^{2} \sqrt {b \,x^{3}+a}}{7}+\frac {2 c x \sqrt {b \,x^{3}+a}}{5}-\frac {2 i a c \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}}-\frac {2 i a d \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}}}}\, \left (\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 ) \operatorname {EllipticE}\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 )+\frac {\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 )}{b}\right )}{7 b \sqrt {b \,x^{3}+a}}\) \(750\)
default \(\text {Expression too large to display}\) \(1480\)

Input:

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

Output:

2/35*x*(5*d*x+7*c)*(b*x^3+a)^(1/2)+3/35*a*(-14/3*I*c*3^(1/2)/b*(-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)*Ellipti 
cF(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))-10/3*I*d*3^(1/2)/b*(-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)*(( 
-3/2/b*(-a*b^2)^(1/3)+1/2*I*3^(1/2)/b*(-a*b^2)^(1/3))*EllipticE(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))+1/b*(-a*b^2)^(1/3)*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.10 (sec) , antiderivative size = 69, normalized size of antiderivative = 0.13 \[ \int \frac {a c+a d x+b c x^3+b d x^4}{\sqrt {a+b x^3}} \, dx=\frac {2 \, {\left (21 \, a \sqrt {b} c {\rm weierstrassPInverse}\left (0, -\frac {4 \, a}{b}, x\right ) - 15 \, a \sqrt {b} d {\rm weierstrassZeta}\left (0, -\frac {4 \, a}{b}, {\rm weierstrassPInverse}\left (0, -\frac {4 \, a}{b}, x\right )\right ) + {\left (5 \, b d x^{2} + 7 \, b c x\right )} \sqrt {b x^{3} + a}\right )}}{35 \, b} \] Input:

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

Output:

2/35*(21*a*sqrt(b)*c*weierstrassPInverse(0, -4*a/b, x) - 15*a*sqrt(b)*d*we 
ierstrassZeta(0, -4*a/b, weierstrassPInverse(0, -4*a/b, x)) + (5*b*d*x^2 + 
 7*b*c*x)*sqrt(b*x^3 + a))/b
 

Sympy [A] (verification not implemented)

Time = 2.23 (sec) , antiderivative size = 163, normalized size of antiderivative = 0.31 \[ \int \frac {a c+a d x+b c x^3+b d x^4}{\sqrt {a+b x^3}} \, dx=\frac {\sqrt {a} c x \Gamma \left (\frac {1}{3}\right ) {{}_{2}F_{1}\left (\begin {matrix} \frac {1}{3}, \frac {1}{2} \\ \frac {4}{3} \end {matrix}\middle | {\frac {b x^{3} e^{i \pi }}{a}} \right )}}{3 \Gamma \left (\frac {4}{3}\right )} + \frac {\sqrt {a} d x^{2} \Gamma \left (\frac {2}{3}\right ) {{}_{2}F_{1}\left (\begin {matrix} \frac {1}{2}, \frac {2}{3} \\ \frac {5}{3} \end {matrix}\middle | {\frac {b x^{3} e^{i \pi }}{a}} \right )}}{3 \Gamma \left (\frac {5}{3}\right )} + \frac {b c 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 \sqrt {a} \Gamma \left (\frac {7}{3}\right )} + \frac {b d x^{5} \Gamma \left (\frac {5}{3}\right ) {{}_{2}F_{1}\left (\begin {matrix} \frac {1}{2}, \frac {5}{3} \\ \frac {8}{3} \end {matrix}\middle | {\frac {b x^{3} e^{i \pi }}{a}} \right )}}{3 \sqrt {a} \Gamma \left (\frac {8}{3}\right )} \] Input:

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

Output:

sqrt(a)*c*x*gamma(1/3)*hyper((1/3, 1/2), (4/3,), b*x**3*exp_polar(I*pi)/a) 
/(3*gamma(4/3)) + sqrt(a)*d*x**2*gamma(2/3)*hyper((1/2, 2/3), (5/3,), b*x* 
*3*exp_polar(I*pi)/a)/(3*gamma(5/3)) + b*c*x**4*gamma(4/3)*hyper((1/2, 4/3 
), (7/3,), b*x**3*exp_polar(I*pi)/a)/(3*sqrt(a)*gamma(7/3)) + b*d*x**5*gam 
ma(5/3)*hyper((1/2, 5/3), (8/3,), b*x**3*exp_polar(I*pi)/a)/(3*sqrt(a)*gam 
ma(8/3))
 

Maxima [F]

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

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

Output:

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

Giac [F]

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

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

Output:

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

Mupad [F(-1)]

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

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

Output:

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

Reduce [F]

\[ \int \frac {a c+a d x+b c x^3+b d x^4}{\sqrt {a+b x^3}} \, dx=\frac {2 \sqrt {b \,x^{3}+a}\, c x}{5}+\frac {2 \sqrt {b \,x^{3}+a}\, d \,x^{2}}{7}+\frac {3 \left (\int \frac {\sqrt {b \,x^{3}+a}}{b \,x^{3}+a}d x \right ) a c}{5}+\frac {3 \left (\int \frac {\sqrt {b \,x^{3}+a}\, x}{b \,x^{3}+a}d x \right ) a d}{7} \] Input:

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

Output:

(14*sqrt(a + b*x**3)*c*x + 10*sqrt(a + b*x**3)*d*x**2 + 21*int(sqrt(a + b* 
x**3)/(a + b*x**3),x)*a*c + 15*int((sqrt(a + b*x**3)*x)/(a + b*x**3),x)*a* 
d)/35