\(\int \frac {1}{(c e+d e x)^{7/2} \sqrt {1-c^2-2 c d x-d^2 x^2}} \, dx\) [307]

Optimal result

Integrand size = 37, antiderivative size = 159 \[ \int \frac {1}{(c e+d e x)^{7/2} \sqrt {1-c^2-2 c d x-d^2 x^2}} \, dx=-\frac {2 \sqrt {1-c^2-2 c d x-d^2 x^2}}{5 d e (c e+d e x)^{5/2}}-\frac {6 \sqrt {1-c^2-2 c d x-d^2 x^2}}{5 d e^3 \sqrt {c e+d e x}}-\frac {6 E\left (\left .\arcsin \left (\frac {\sqrt {c e+d e x}}{\sqrt {e}}\right )\right |-1\right )}{5 d e^{7/2}}+\frac {6 \operatorname {EllipticF}\left (\arcsin \left (\frac {\sqrt {c e+d e x}}{\sqrt {e}}\right ),-1\right )}{5 d e^{7/2}} \] Output:

-2/5*(-d^2*x^2-2*c*d*x-c^2+1)^(1/2)/d/e/(d*e*x+c*e)^(5/2)-6/5*(-d^2*x^2-2* 
c*d*x-c^2+1)^(1/2)/d/e^3/(d*e*x+c*e)^(1/2)-6/5*EllipticE((d*e*x+c*e)^(1/2) 
/e^(1/2),I)/d/e^(7/2)+6/5*EllipticF((d*e*x+c*e)^(1/2)/e^(1/2),I)/d/e^(7/2)
 

Mathematica [C] (verified)

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

Time = 10.03 (sec) , antiderivative size = 40, normalized size of antiderivative = 0.25 \[ \int \frac {1}{(c e+d e x)^{7/2} \sqrt {1-c^2-2 c d x-d^2 x^2}} \, dx=-\frac {2 (c+d x) \operatorname {Hypergeometric2F1}\left (-\frac {5}{4},\frac {1}{2},-\frac {1}{4},(c+d x)^2\right )}{5 d (e (c+d x))^{7/2}} \] Input:

Integrate[1/((c*e + d*e*x)^(7/2)*Sqrt[1 - c^2 - 2*c*d*x - d^2*x^2]),x]
 

Output:

(-2*(c + d*x)*Hypergeometric2F1[-5/4, 1/2, -1/4, (c + d*x)^2])/(5*d*(e*(c 
+ d*x))^(7/2))
 

Rubi [A] (verified)

Time = 0.37 (sec) , antiderivative size = 163, normalized size of antiderivative = 1.03, number of steps used = 9, number of rules used = 8, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.216, Rules used = {1117, 1117, 1114, 836, 27, 762, 1389, 327}

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 {1}{\sqrt {-c^2-2 c d x-d^2 x^2+1} (c e+d e x)^{7/2}} \, dx\)

\(\Big \downarrow \) 1117

\(\displaystyle \frac {3 \int \frac {1}{(c e+d x e)^{3/2} \sqrt {-c^2-2 d x c-d^2 x^2+1}}dx}{5 e^2}-\frac {2 \sqrt {-c^2-2 c d x-d^2 x^2+1}}{5 d e (c e+d e x)^{5/2}}\)

\(\Big \downarrow \) 1117

\(\displaystyle \frac {3 \left (-\frac {\int \frac {\sqrt {c e+d x e}}{\sqrt {-c^2-2 d x c-d^2 x^2+1}}dx}{e^2}-\frac {2 \sqrt {-c^2-2 c d x-d^2 x^2+1}}{d e \sqrt {c e+d e x}}\right )}{5 e^2}-\frac {2 \sqrt {-c^2-2 c d x-d^2 x^2+1}}{5 d e (c e+d e x)^{5/2}}\)

\(\Big \downarrow \) 1114

\(\displaystyle \frac {3 \left (-\frac {2 \int \frac {c e+d x e}{\sqrt {1-\frac {(c e+d x e)^2}{e^2}}}d\sqrt {c e+d x e}}{d e^3}-\frac {2 \sqrt {-c^2-2 c d x-d^2 x^2+1}}{d e \sqrt {c e+d e x}}\right )}{5 e^2}-\frac {2 \sqrt {-c^2-2 c d x-d^2 x^2+1}}{5 d e (c e+d e x)^{5/2}}\)

\(\Big \downarrow \) 836

\(\displaystyle \frac {3 \left (-\frac {2 \left (e \int \frac {c e+d x e+e}{e \sqrt {1-\frac {(c e+d x e)^2}{e^2}}}d\sqrt {c e+d x e}-e \int \frac {1}{\sqrt {1-\frac {(c e+d x e)^2}{e^2}}}d\sqrt {c e+d x e}\right )}{d e^3}-\frac {2 \sqrt {-c^2-2 c d x-d^2 x^2+1}}{d e \sqrt {c e+d e x}}\right )}{5 e^2}-\frac {2 \sqrt {-c^2-2 c d x-d^2 x^2+1}}{5 d e (c e+d e x)^{5/2}}\)

\(\Big \downarrow \) 27

\(\displaystyle \frac {3 \left (-\frac {2 \left (\int \frac {c e+d x e+e}{\sqrt {1-\frac {(c e+d x e)^2}{e^2}}}d\sqrt {c e+d x e}-e \int \frac {1}{\sqrt {1-\frac {(c e+d x e)^2}{e^2}}}d\sqrt {c e+d x e}\right )}{d e^3}-\frac {2 \sqrt {-c^2-2 c d x-d^2 x^2+1}}{d e \sqrt {c e+d e x}}\right )}{5 e^2}-\frac {2 \sqrt {-c^2-2 c d x-d^2 x^2+1}}{5 d e (c e+d e x)^{5/2}}\)

\(\Big \downarrow \) 762

\(\displaystyle \frac {3 \left (-\frac {2 \left (\int \frac {c e+d x e+e}{\sqrt {1-\frac {(c e+d x e)^2}{e^2}}}d\sqrt {c e+d x e}-e^{3/2} \operatorname {EllipticF}\left (\arcsin \left (\frac {\sqrt {c e+d x e}}{\sqrt {e}}\right ),-1\right )\right )}{d e^3}-\frac {2 \sqrt {-c^2-2 c d x-d^2 x^2+1}}{d e \sqrt {c e+d e x}}\right )}{5 e^2}-\frac {2 \sqrt {-c^2-2 c d x-d^2 x^2+1}}{5 d e (c e+d e x)^{5/2}}\)

\(\Big \downarrow \) 1389

\(\displaystyle \frac {3 \left (-\frac {2 \left (e \int \frac {\sqrt {\frac {c e+d x e}{e}+1}}{\sqrt {1-\frac {c e+d x e}{e}}}d\sqrt {c e+d x e}-e^{3/2} \operatorname {EllipticF}\left (\arcsin \left (\frac {\sqrt {c e+d x e}}{\sqrt {e}}\right ),-1\right )\right )}{d e^3}-\frac {2 \sqrt {-c^2-2 c d x-d^2 x^2+1}}{d e \sqrt {c e+d e x}}\right )}{5 e^2}-\frac {2 \sqrt {-c^2-2 c d x-d^2 x^2+1}}{5 d e (c e+d e x)^{5/2}}\)

\(\Big \downarrow \) 327

\(\displaystyle \frac {3 \left (-\frac {2 \left (e^{3/2} E\left (\left .\arcsin \left (\frac {\sqrt {c e+d x e}}{\sqrt {e}}\right )\right |-1\right )-e^{3/2} \operatorname {EllipticF}\left (\arcsin \left (\frac {\sqrt {c e+d x e}}{\sqrt {e}}\right ),-1\right )\right )}{d e^3}-\frac {2 \sqrt {-c^2-2 c d x-d^2 x^2+1}}{d e \sqrt {c e+d e x}}\right )}{5 e^2}-\frac {2 \sqrt {-c^2-2 c d x-d^2 x^2+1}}{5 d e (c e+d e x)^{5/2}}\)

Input:

Int[1/((c*e + d*e*x)^(7/2)*Sqrt[1 - c^2 - 2*c*d*x - d^2*x^2]),x]
 

Output:

(-2*Sqrt[1 - c^2 - 2*c*d*x - d^2*x^2])/(5*d*e*(c*e + d*e*x)^(5/2)) + (3*(( 
-2*Sqrt[1 - c^2 - 2*c*d*x - d^2*x^2])/(d*e*Sqrt[c*e + d*e*x]) - (2*(e^(3/2 
)*EllipticE[ArcSin[Sqrt[c*e + d*e*x]/Sqrt[e]], -1] - e^(3/2)*EllipticF[Arc 
Sin[Sqrt[c*e + d*e*x]/Sqrt[e]], -1]))/(d*e^3)))/(5*e^2)
 

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 327

Int[Sqrt[(a_) + (b_.)*(x_)^2]/Sqrt[(c_) + (d_.)*(x_)^2], x_Symbol] :> Simp[ 
(Sqrt[a]/(Sqrt[c]*Rt[-d/c, 2]))*EllipticE[ArcSin[Rt[-d/c, 2]*x], b*(c/(a*d) 
)], x] /; FreeQ[{a, b, c, d}, x] && NegQ[d/c] && GtQ[c, 0] && GtQ[a, 0]
 

rule 762

Int[1/Sqrt[(a_) + (b_.)*(x_)^4], x_Symbol] :> Simp[(1/(Sqrt[a]*Rt[-b/a, 4]) 
)*EllipticF[ArcSin[Rt[-b/a, 4]*x], -1], x] /; FreeQ[{a, b}, x] && NegQ[b/a] 
 && GtQ[a, 0]
 

rule 836

Int[(x_)^2/Sqrt[(a_) + (b_.)*(x_)^4], x_Symbol] :> With[{q = Rt[-b/a, 2]}, 
Simp[-q^(-1)   Int[1/Sqrt[a + b*x^4], x], x] + Simp[1/q   Int[(1 + q*x^2)/S 
qrt[a + b*x^4], x], x]] /; FreeQ[{a, b}, x] && NegQ[b/a]
 

rule 1114

Int[Sqrt[(d_) + (e_.)*(x_)]/Sqrt[(a_.) + (b_.)*(x_) + (c_.)*(x_)^2], x_Symb 
ol] :> Simp[(4/e)*Sqrt[-c/(b^2 - 4*a*c)]   Subst[Int[x^2/Sqrt[Simp[1 - b^2* 
(x^4/(d^2*(b^2 - 4*a*c))), x]], x], x, Sqrt[d + e*x]], x] /; FreeQ[{a, b, c 
, d, e}, x] && EqQ[2*c*d - b*e, 0] && LtQ[c/(b^2 - 4*a*c), 0]
 

rule 1117

Int[((d_) + (e_.)*(x_))^(m_)*((a_.) + (b_.)*(x_) + (c_.)*(x_)^2)^(p_.), x_S 
ymbol] :> Simp[-2*b*d*(d + e*x)^(m + 1)*((a + b*x + c*x^2)^(p + 1)/(d^2*(m 
+ 1)*(b^2 - 4*a*c))), x] + Simp[b^2*((m + 2*p + 3)/(d^2*(m + 1)*(b^2 - 4*a* 
c)))   Int[(d + e*x)^(m + 2)*(a + b*x + c*x^2)^p, x], x] /; FreeQ[{a, b, c, 
 d, e, p}, x] && EqQ[2*c*d - b*e, 0] && NeQ[m + 2*p + 3, 0] && LtQ[m, -1] & 
& (IntegerQ[2*p] || (IntegerQ[m] && RationalQ[p]) || IntegerQ[(m + 2*p + 3) 
/2])
 

rule 1389

Int[((d_) + (e_.)*(x_)^2)/Sqrt[(a_) + (c_.)*(x_)^4], x_Symbol] :> Simp[d/Sq 
rt[a]   Int[Sqrt[1 + e*(x^2/d)]/Sqrt[1 - e*(x^2/d)], x], x] /; FreeQ[{a, c, 
 d, e}, x] && EqQ[c*d^2 + a*e^2, 0] && NegQ[c/a] && GtQ[a, 0]
 

Maple [A] (verified)

Time = 4.83 (sec) , antiderivative size = 263, normalized size of antiderivative = 1.65

method	result	size
default	\(\frac {\left (3 \operatorname {EllipticE}\left (\frac {\sqrt {-2 d x -2 c +2}}{2}, \sqrt {2}\right ) d^{2} x^{2} \sqrt {-2 d x -2 c +2}\, \sqrt {d x +c}\, \sqrt {2 d x +2 c +2}+6 d^{4} x^{4}+6 \operatorname {EllipticE}\left (\frac {\sqrt {-2 d x -2 c +2}}{2}, \sqrt {2}\right ) c d x \sqrt {-2 d x -2 c +2}\, \sqrt {d x +c}\, \sqrt {2 d x +2 c +2}+24 c \,d^{3} x^{3}+3 \operatorname {EllipticE}\left (\frac {\sqrt {-2 d x -2 c +2}}{2}, \sqrt {2}\right ) c^{2} \sqrt {-2 d x -2 c +2}\, \sqrt {d x +c}\, \sqrt {2 d x +2 c +2}+36 d^{2} c^{2} x^{2}+24 c^{3} d x +6 c^{4}-4 d^{2} x^{2}-8 c d x -4 c^{2}-2\right ) \sqrt {e \left (d x +c \right )}}{5 \left (d x +c \right )^{3} \sqrt {-d^{2} x^{2}-2 c d x -c^{2}+1}\, e^{4} d}\)	\(263\)
elliptic	\(\frac {\sqrt {-e \left (d x +c \right ) \left (d^{2} x^{2}+2 c d x +c^{2}-1\right )}\, \left (-\frac {2 \sqrt {-d^{3} e \,x^{3}-3 c \,d^{2} e \,x^{2}-3 c^{2} d e x -c^{3} e +d e x +c e}}{5 d^{4} e^{4} \left (x +\frac {c}{d}\right )^{3}}-\frac {6 \left (-d^{3} e \,x^{2}-2 c \,d^{2} e x -c^{2} d e +d e \right )}{5 d^{2} e^{4} \sqrt {\left (x +\frac {c}{d}\right ) \left (-d^{3} e \,x^{2}-2 c \,d^{2} e x -c^{2} d e +d e \right )}}-\frac {6 c \left (-\frac {c +1}{d}+\frac {c -1}{d}\right ) \sqrt {\frac {x +\frac {c -1}{d}}{-\frac {c +1}{d}+\frac {c -1}{d}}}\, \sqrt {\frac {x +\frac {c}{d}}{-\frac {c -1}{d}+\frac {c}{d}}}\, \sqrt {\frac {x +\frac {c +1}{d}}{-\frac {c -1}{d}+\frac {c +1}{d}}}\, \operatorname {EllipticF}\left (\sqrt {\frac {x +\frac {c -1}{d}}{-\frac {c +1}{d}+\frac {c -1}{d}}}, \sqrt {\frac {-\frac {c -1}{d}+\frac {c +1}{d}}{-\frac {c -1}{d}+\frac {c}{d}}}\right )}{5 e^{3} \sqrt {-d^{3} e \,x^{3}-3 c \,d^{2} e \,x^{2}-3 c^{2} d e x -c^{3} e +d e x +c e}}-\frac {6 d \left (-\frac {c +1}{d}+\frac {c -1}{d}\right ) \sqrt {\frac {x +\frac {c -1}{d}}{-\frac {c +1}{d}+\frac {c -1}{d}}}\, \sqrt {\frac {x +\frac {c}{d}}{-\frac {c -1}{d}+\frac {c}{d}}}\, \sqrt {\frac {x +\frac {c +1}{d}}{-\frac {c -1}{d}+\frac {c +1}{d}}}\, \left (\left (-\frac {c -1}{d}+\frac {c}{d}\right ) \operatorname {EllipticE}\left (\sqrt {\frac {x +\frac {c -1}{d}}{-\frac {c +1}{d}+\frac {c -1}{d}}}, \sqrt {\frac {-\frac {c -1}{d}+\frac {c +1}{d}}{-\frac {c -1}{d}+\frac {c}{d}}}\right )-\frac {c \operatorname {EllipticF}\left (\sqrt {\frac {x +\frac {c -1}{d}}{-\frac {c +1}{d}+\frac {c -1}{d}}}, \sqrt {\frac {-\frac {c -1}{d}+\frac {c +1}{d}}{-\frac {c -1}{d}+\frac {c}{d}}}\right )}{d}\right )}{5 e^{3} \sqrt {-d^{3} e \,x^{3}-3 c \,d^{2} e \,x^{2}-3 c^{2} d e x -c^{3} e +d e x +c e}}\right )}{\sqrt {e \left (d x +c \right )}\, \sqrt {-d^{2} x^{2}-2 c d x -c^{2}+1}}\)	\(717\)

Input:

int(1/(d*e*x+c*e)^(7/2)/(-d^2*x^2-2*c*d*x-c^2+1)^(1/2),x,method=_RETURNVER 
BOSE)
 

Output:

1/5*(3*EllipticE(1/2*(-2*d*x-2*c+2)^(1/2),2^(1/2))*d^2*x^2*(-2*d*x-2*c+2)^ 
(1/2)*(d*x+c)^(1/2)*(2*d*x+2*c+2)^(1/2)+6*d^4*x^4+6*EllipticE(1/2*(-2*d*x- 
2*c+2)^(1/2),2^(1/2))*c*d*x*(-2*d*x-2*c+2)^(1/2)*(d*x+c)^(1/2)*(2*d*x+2*c+ 
2)^(1/2)+24*c*d^3*x^3+3*EllipticE(1/2*(-2*d*x-2*c+2)^(1/2),2^(1/2))*c^2*(- 
2*d*x-2*c+2)^(1/2)*(d*x+c)^(1/2)*(2*d*x+2*c+2)^(1/2)+36*d^2*c^2*x^2+24*c^3 
*d*x+6*c^4-4*d^2*x^2-8*c*d*x-4*c^2-2)/(d*x+c)^3/(-d^2*x^2-2*c*d*x-c^2+1)^( 
1/2)*(e*(d*x+c))^(1/2)/e^4/d
 

Fricas [A] (verification not implemented)

Time = 0.11 (sec) , antiderivative size = 168, normalized size of antiderivative = 1.06 \[ \int \frac {1}{(c e+d e x)^{7/2} \sqrt {1-c^2-2 c d x-d^2 x^2}} \, dx=-\frac {2 \, {\left (3 \, {\left (d^{3} x^{3} + 3 \, c d^{2} x^{2} + 3 \, c^{2} d x + c^{3}\right )} \sqrt {-d^{3} e} {\rm weierstrassZeta}\left (\frac {4}{d^{2}}, 0, {\rm weierstrassPInverse}\left (\frac {4}{d^{2}}, 0, \frac {d x + c}{d}\right )\right ) + {\left (3 \, d^{3} x^{2} + 6 \, c d^{2} x + {\left (3 \, c^{2} + 1\right )} d\right )} \sqrt {-d^{2} x^{2} - 2 \, c d x - c^{2} + 1} \sqrt {d e x + c e}\right )}}{5 \, {\left (d^{5} e^{4} x^{3} + 3 \, c d^{4} e^{4} x^{2} + 3 \, c^{2} d^{3} e^{4} x + c^{3} d^{2} e^{4}\right )}} \] Input:

integrate(1/(d*e*x+c*e)^(7/2)/(-d^2*x^2-2*c*d*x-c^2+1)^(1/2),x, algorithm= 
"fricas")
 

Output:

-2/5*(3*(d^3*x^3 + 3*c*d^2*x^2 + 3*c^2*d*x + c^3)*sqrt(-d^3*e)*weierstrass 
Zeta(4/d^2, 0, weierstrassPInverse(4/d^2, 0, (d*x + c)/d)) + (3*d^3*x^2 + 
6*c*d^2*x + (3*c^2 + 1)*d)*sqrt(-d^2*x^2 - 2*c*d*x - c^2 + 1)*sqrt(d*e*x + 
 c*e))/(d^5*e^4*x^3 + 3*c*d^4*e^4*x^2 + 3*c^2*d^3*e^4*x + c^3*d^2*e^4)
 

Sympy [F]

\[ \int \frac {1}{(c e+d e x)^{7/2} \sqrt {1-c^2-2 c d x-d^2 x^2}} \, dx=\int \frac {1}{\left (e \left (c + d x\right )\right )^{\frac {7}{2}} \sqrt {- \left (c + d x - 1\right ) \left (c + d x + 1\right )}}\, dx \] Input:

integrate(1/(d*e*x+c*e)**(7/2)/(-d**2*x**2-2*c*d*x-c**2+1)**(1/2),x)
 

Output:

Integral(1/((e*(c + d*x))**(7/2)*sqrt(-(c + d*x - 1)*(c + d*x + 1))), x)
 

Maxima [F]

\[ \int \frac {1}{(c e+d e x)^{7/2} \sqrt {1-c^2-2 c d x-d^2 x^2}} \, dx=\int { \frac {1}{\sqrt {-d^{2} x^{2} - 2 \, c d x - c^{2} + 1} {\left (d e x + c e\right )}^{\frac {7}{2}}} \,d x } \] Input:

integrate(1/(d*e*x+c*e)^(7/2)/(-d^2*x^2-2*c*d*x-c^2+1)^(1/2),x, algorithm= 
"maxima")
 

Output:

integrate(1/(sqrt(-d^2*x^2 - 2*c*d*x - c^2 + 1)*(d*e*x + c*e)^(7/2)), x)
 

Giac [F]

\[ \int \frac {1}{(c e+d e x)^{7/2} \sqrt {1-c^2-2 c d x-d^2 x^2}} \, dx=\int { \frac {1}{\sqrt {-d^{2} x^{2} - 2 \, c d x - c^{2} + 1} {\left (d e x + c e\right )}^{\frac {7}{2}}} \,d x } \] Input:

integrate(1/(d*e*x+c*e)^(7/2)/(-d^2*x^2-2*c*d*x-c^2+1)^(1/2),x, algorithm= 
"giac")
 

Output:

integrate(1/(sqrt(-d^2*x^2 - 2*c*d*x - c^2 + 1)*(d*e*x + c*e)^(7/2)), x)
 

Mupad [F(-1)]

Timed out. \[ \int \frac {1}{(c e+d e x)^{7/2} \sqrt {1-c^2-2 c d x-d^2 x^2}} \, dx=\int \frac {1}{{\left (c\,e+d\,e\,x\right )}^{7/2}\,\sqrt {-c^2-2\,c\,d\,x-d^2\,x^2+1}} \,d x \] Input:

int(1/((c*e + d*e*x)^(7/2)*(1 - d^2*x^2 - 2*c*d*x - c^2)^(1/2)),x)
 

Output:

int(1/((c*e + d*e*x)^(7/2)*(1 - d^2*x^2 - 2*c*d*x - c^2)^(1/2)), x)
 

Reduce [F]

\[ \int \frac {1}{(c e+d e x)^{7/2} \sqrt {1-c^2-2 c d x-d^2 x^2}} \, dx =\text {Too large to display} \] Input:

int(1/(d*e*x+c*e)^(7/2)/(-d^2*x^2-2*c*d*x-c^2+1)^(1/2),x)
 

Output:

(sqrt(e)*(int((sqrt(c + d*x)*sqrt( - c**2 - 2*c*d*x - d**2*x**2 + 1)*x**3) 
/(c**10 + 6*c**9*d*x + 15*c**8*d**2*x**2 - c**8 + 20*c**7*d**3*x**3 - 4*c* 
*7*d*x + 15*c**6*d**4*x**4 - 6*c**6*d**2*x**2 - c**6 + 6*c**5*d**5*x**5 - 
4*c**5*d**3*x**3 - 6*c**5*d*x + c**4*d**6*x**6 - c**4*d**4*x**4 - 15*c**4* 
d**2*x**2 + c**4 - 20*c**3*d**3*x**3 + 4*c**3*d*x - 15*c**2*d**4*x**4 + 6* 
c**2*d**2*x**2 - 6*c*d**5*x**5 + 4*c*d**3*x**3 - d**6*x**6 + d**4*x**4),x) 
*c**2*d**3 + int((sqrt(c + d*x)*sqrt( - c**2 - 2*c*d*x - d**2*x**2 + 1)*x* 
*3)/(c**10 + 6*c**9*d*x + 15*c**8*d**2*x**2 - c**8 + 20*c**7*d**3*x**3 - 4 
*c**7*d*x + 15*c**6*d**4*x**4 - 6*c**6*d**2*x**2 - c**6 + 6*c**5*d**5*x**5 
 - 4*c**5*d**3*x**3 - 6*c**5*d*x + c**4*d**6*x**6 - c**4*d**4*x**4 - 15*c* 
*4*d**2*x**2 + c**4 - 20*c**3*d**3*x**3 + 4*c**3*d*x - 15*c**2*d**4*x**4 + 
 6*c**2*d**2*x**2 - 6*c*d**5*x**5 + 4*c*d**3*x**3 - d**6*x**6 + d**4*x**4) 
,x)*d**3 - int((sqrt(c + d*x)*sqrt( - c**2 - 2*c*d*x - d**2*x**2 + 1)*x**3 
)/(c**8 + 6*c**7*d*x + 15*c**6*d**2*x**2 - 2*c**6 + 20*c**5*d**3*x**3 - 10 
*c**5*d*x + 15*c**4*d**4*x**4 - 21*c**4*d**2*x**2 + c**4 + 6*c**3*d**5*x** 
5 - 24*c**3*d**3*x**3 + 4*c**3*d*x + c**2*d**6*x**6 - 16*c**2*d**4*x**4 + 
6*c**2*d**2*x**2 - 6*c*d**5*x**5 + 4*c*d**3*x**3 - d**6*x**6 + d**4*x**4), 
x)*d**3 - 2*int((sqrt(c + d*x)*sqrt( - c**2 - 2*c*d*x - d**2*x**2 + 1))/(c 
**10 + 6*c**9*d*x + 15*c**8*d**2*x**2 - c**8 + 20*c**7*d**3*x**3 - 4*c**7* 
d*x + 15*c**6*d**4*x**4 - 6*c**6*d**2*x**2 - c**6 + 6*c**5*d**5*x**5 - ...