\(\int (f+g x)^2 \sqrt {d-c^2 d x^2} (a+b \arcsin (c x))^2 \, dx\) [132]

Optimal result

Integrand size = 33, antiderivative size = 723 \[ \int (f+g x)^2 \sqrt {d-c^2 d x^2} (a+b \arcsin (c x))^2 \, dx=\frac {8 b^2 f g \sqrt {d-c^2 d x^2}}{9 c^2}-\frac {1}{4} b^2 f^2 x \sqrt {d-c^2 d x^2}+\frac {b^2 g^2 x \sqrt {d-c^2 d x^2}}{64 c^2}-\frac {1}{32} b^2 g^2 x^3 \sqrt {d-c^2 d x^2}+\frac {4 b^2 f g \left (d-c^2 d x^2\right )^{3/2}}{27 c^2 d}+\frac {b^2 f^2 \sqrt {d-c^2 d x^2} \arcsin (c x)}{4 c \sqrt {1-c^2 x^2}}-\frac {b^2 g^2 \sqrt {d-c^2 d x^2} \arcsin (c x)}{64 c^3 \sqrt {1-c^2 x^2}}+\frac {4 b f g x \sqrt {d-c^2 d x^2} (a+b \arcsin (c x))}{3 c \sqrt {1-c^2 x^2}}-\frac {b c f^2 x^2 \sqrt {d-c^2 d x^2} (a+b \arcsin (c x))}{2 \sqrt {1-c^2 x^2}}+\frac {b g^2 x^2 \sqrt {d-c^2 d x^2} (a+b \arcsin (c x))}{8 c \sqrt {1-c^2 x^2}}-\frac {4 b c f g x^3 \sqrt {d-c^2 d x^2} (a+b \arcsin (c x))}{9 \sqrt {1-c^2 x^2}}-\frac {b c g^2 x^4 \sqrt {d-c^2 d x^2} (a+b \arcsin (c x))}{8 \sqrt {1-c^2 x^2}}+\frac {1}{2} f^2 x \sqrt {d-c^2 d x^2} (a+b \arcsin (c x))^2-\frac {g^2 x \sqrt {d-c^2 d x^2} (a+b \arcsin (c x))^2}{8 c^2}+\frac {1}{4} g^2 x^3 \sqrt {d-c^2 d x^2} (a+b \arcsin (c x))^2-\frac {2 f g \left (d-c^2 d x^2\right )^{3/2} (a+b \arcsin (c x))^2}{3 c^2 d}+\frac {f^2 \sqrt {d-c^2 d x^2} (a+b \arcsin (c x))^3}{6 b c \sqrt {1-c^2 x^2}}+\frac {g^2 \sqrt {d-c^2 d x^2} (a+b \arcsin (c x))^3}{24 b c^3 \sqrt {1-c^2 x^2}} \] Output:

8/9*b^2*f*g*(-c^2*d*x^2+d)^(1/2)/c^2-1/4*b^2*f^2*x*(-c^2*d*x^2+d)^(1/2)+1/ 
64*b^2*g^2*x*(-c^2*d*x^2+d)^(1/2)/c^2-1/32*b^2*g^2*x^3*(-c^2*d*x^2+d)^(1/2 
)+4/27*b^2*f*g*(-c^2*d*x^2+d)^(3/2)/c^2/d+1/4*b^2*f^2*(-c^2*d*x^2+d)^(1/2) 
*arcsin(c*x)/c/(-c^2*x^2+1)^(1/2)-1/64*b^2*g^2*(-c^2*d*x^2+d)^(1/2)*arcsin 
(c*x)/c^3/(-c^2*x^2+1)^(1/2)+4/3*b*f*g*x*(-c^2*d*x^2+d)^(1/2)*(a+b*arcsin( 
c*x))/c/(-c^2*x^2+1)^(1/2)-1/2*b*c*f^2*x^2*(-c^2*d*x^2+d)^(1/2)*(a+b*arcsi 
n(c*x))/(-c^2*x^2+1)^(1/2)+1/8*b*g^2*x^2*(-c^2*d*x^2+d)^(1/2)*(a+b*arcsin( 
c*x))/c/(-c^2*x^2+1)^(1/2)-4/9*b*c*f*g*x^3*(-c^2*d*x^2+d)^(1/2)*(a+b*arcsi 
n(c*x))/(-c^2*x^2+1)^(1/2)-1/8*b*c*g^2*x^4*(-c^2*d*x^2+d)^(1/2)*(a+b*arcsi 
n(c*x))/(-c^2*x^2+1)^(1/2)+1/2*f^2*x*(-c^2*d*x^2+d)^(1/2)*(a+b*arcsin(c*x) 
)^2-1/8*g^2*x*(-c^2*d*x^2+d)^(1/2)*(a+b*arcsin(c*x))^2/c^2+1/4*g^2*x^3*(-c 
^2*d*x^2+d)^(1/2)*(a+b*arcsin(c*x))^2-2/3*f*g*(-c^2*d*x^2+d)^(3/2)*(a+b*ar 
csin(c*x))^2/c^2/d+1/6*f^2*(-c^2*d*x^2+d)^(1/2)*(a+b*arcsin(c*x))^3/b/c/(- 
c^2*x^2+1)^(1/2)+1/24*g^2*(-c^2*d*x^2+d)^(1/2)*(a+b*arcsin(c*x))^3/b/c^3/( 
-c^2*x^2+1)^(1/2)
 

Mathematica [A] (verified)

Time = 0.53 (sec) , antiderivative size = 441, normalized size of antiderivative = 0.61 \[ \int (f+g x)^2 \sqrt {d-c^2 d x^2} (a+b \arcsin (c x))^2 \, dx=\frac {\sqrt {d-c^2 d x^2} \left (\frac {1}{2} f^2 x \sqrt {1-c^2 x^2} (a+b \arcsin (c x))^2+\frac {1}{4} g^2 x^3 \sqrt {1-c^2 x^2} (a+b \arcsin (c x))^2-\frac {2 f g \left (1-c^2 x^2\right )^{3/2} (a+b \arcsin (c x))^2}{3 c^2}+\frac {f^2 (a+b \arcsin (c x))^3}{6 b c}-\frac {4 b f g \left (b \sqrt {1-c^2 x^2} \left (-7+c^2 x^2\right )+3 a c x \left (-3+c^2 x^2\right )+3 b c x \left (-3+c^2 x^2\right ) \arcsin (c x)\right )}{27 c^2}-\frac {b f^2 \left (c x \left (2 a c x+b \sqrt {1-c^2 x^2}\right )+b \left (-1+2 c^2 x^2\right ) \arcsin (c x)\right )}{4 c}-\frac {b g^2 \left (8 a c^4 x^4+b c x \sqrt {1-c^2 x^2} \left (3+2 c^2 x^2\right )+b \left (-3+8 c^4 x^4\right ) \arcsin (c x)\right )}{64 c^3}-\frac {g^2 \left (6 b c x \sqrt {1-c^2 x^2} (a+b \arcsin (c x))^2-2 (a+b \arcsin (c x))^3-3 b^2 \left (c x \left (2 a c x+b \sqrt {1-c^2 x^2}\right )+b \left (-1+2 c^2 x^2\right ) \arcsin (c x)\right )\right )}{48 b c^3}\right )}{\sqrt {1-c^2 x^2}} \] Input:

Integrate[(f + g*x)^2*Sqrt[d - c^2*d*x^2]*(a + b*ArcSin[c*x])^2,x]
 

Output:

(Sqrt[d - c^2*d*x^2]*((f^2*x*Sqrt[1 - c^2*x^2]*(a + b*ArcSin[c*x])^2)/2 + 
(g^2*x^3*Sqrt[1 - c^2*x^2]*(a + b*ArcSin[c*x])^2)/4 - (2*f*g*(1 - c^2*x^2) 
^(3/2)*(a + b*ArcSin[c*x])^2)/(3*c^2) + (f^2*(a + b*ArcSin[c*x])^3)/(6*b*c 
) - (4*b*f*g*(b*Sqrt[1 - c^2*x^2]*(-7 + c^2*x^2) + 3*a*c*x*(-3 + c^2*x^2) 
+ 3*b*c*x*(-3 + c^2*x^2)*ArcSin[c*x]))/(27*c^2) - (b*f^2*(c*x*(2*a*c*x + b 
*Sqrt[1 - c^2*x^2]) + b*(-1 + 2*c^2*x^2)*ArcSin[c*x]))/(4*c) - (b*g^2*(8*a 
*c^4*x^4 + b*c*x*Sqrt[1 - c^2*x^2]*(3 + 2*c^2*x^2) + b*(-3 + 8*c^4*x^4)*Ar 
cSin[c*x]))/(64*c^3) - (g^2*(6*b*c*x*Sqrt[1 - c^2*x^2]*(a + b*ArcSin[c*x]) 
^2 - 2*(a + b*ArcSin[c*x])^3 - 3*b^2*(c*x*(2*a*c*x + b*Sqrt[1 - c^2*x^2]) 
+ b*(-1 + 2*c^2*x^2)*ArcSin[c*x])))/(48*b*c^3)))/Sqrt[1 - c^2*x^2]
 

Rubi [A] (verified)

Time = 1.25 (sec) , antiderivative size = 477, normalized size of antiderivative = 0.66, number of steps used = 3, number of rules used = 3, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.091, Rules used = {5276, 5262, 2009}

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 {d-c^2 d x^2} (f+g x)^2 (a+b \arcsin (c x))^2 \, dx\)

\(\Big \downarrow \) 5276

\(\displaystyle \frac {\sqrt {d-c^2 d x^2} \int (f+g x)^2 \sqrt {1-c^2 x^2} (a+b \arcsin (c x))^2dx}{\sqrt {1-c^2 x^2}}\)

\(\Big \downarrow \) 5262

\(\displaystyle \frac {\sqrt {d-c^2 d x^2} \int \left (f^2 \sqrt {1-c^2 x^2} (a+b \arcsin (c x))^2+g^2 x^2 \sqrt {1-c^2 x^2} (a+b \arcsin (c x))^2+2 f g x \sqrt {1-c^2 x^2} (a+b \arcsin (c x))^2\right )dx}{\sqrt {1-c^2 x^2}}\)

\(\Big \downarrow \) 2009

\(\displaystyle \frac {\sqrt {d-c^2 d x^2} \left (\frac {g^2 (a+b \arcsin (c x))^3}{24 b c^3}+\frac {1}{2} f^2 x \sqrt {1-c^2 x^2} (a+b \arcsin (c x))^2-\frac {2 f g \left (1-c^2 x^2\right )^{3/2} (a+b \arcsin (c x))^2}{3 c^2}-\frac {g^2 x \sqrt {1-c^2 x^2} (a+b \arcsin (c x))^2}{8 c^2}+\frac {1}{4} g^2 x^3 \sqrt {1-c^2 x^2} (a+b \arcsin (c x))^2-\frac {1}{2} b c f^2 x^2 (a+b \arcsin (c x))+\frac {f^2 (a+b \arcsin (c x))^3}{6 b c}-\frac {4}{9} b c f g x^3 (a+b \arcsin (c x))+\frac {4 b f g x (a+b \arcsin (c x))}{3 c}-\frac {1}{8} b c g^2 x^4 (a+b \arcsin (c x))+\frac {b g^2 x^2 (a+b \arcsin (c x))}{8 c}-\frac {b^2 g^2 \arcsin (c x)}{64 c^3}+\frac {b^2 f^2 \arcsin (c x)}{4 c}-\frac {1}{4} b^2 f^2 x \sqrt {1-c^2 x^2}+\frac {4 b^2 f g \left (1-c^2 x^2\right )^{3/2}}{27 c^2}+\frac {8 b^2 f g \sqrt {1-c^2 x^2}}{9 c^2}+\frac {b^2 g^2 x \sqrt {1-c^2 x^2}}{64 c^2}-\frac {1}{32} b^2 g^2 x^3 \sqrt {1-c^2 x^2}\right )}{\sqrt {1-c^2 x^2}}\)

Input:

Int[(f + g*x)^2*Sqrt[d - c^2*d*x^2]*(a + b*ArcSin[c*x])^2,x]
 

Output:

(Sqrt[d - c^2*d*x^2]*((8*b^2*f*g*Sqrt[1 - c^2*x^2])/(9*c^2) - (b^2*f^2*x*S 
qrt[1 - c^2*x^2])/4 + (b^2*g^2*x*Sqrt[1 - c^2*x^2])/(64*c^2) - (b^2*g^2*x^ 
3*Sqrt[1 - c^2*x^2])/32 + (4*b^2*f*g*(1 - c^2*x^2)^(3/2))/(27*c^2) + (b^2* 
f^2*ArcSin[c*x])/(4*c) - (b^2*g^2*ArcSin[c*x])/(64*c^3) + (4*b*f*g*x*(a + 
b*ArcSin[c*x]))/(3*c) - (b*c*f^2*x^2*(a + b*ArcSin[c*x]))/2 + (b*g^2*x^2*( 
a + b*ArcSin[c*x]))/(8*c) - (4*b*c*f*g*x^3*(a + b*ArcSin[c*x]))/9 - (b*c*g 
^2*x^4*(a + b*ArcSin[c*x]))/8 + (f^2*x*Sqrt[1 - c^2*x^2]*(a + b*ArcSin[c*x 
])^2)/2 - (g^2*x*Sqrt[1 - c^2*x^2]*(a + b*ArcSin[c*x])^2)/(8*c^2) + (g^2*x 
^3*Sqrt[1 - c^2*x^2]*(a + b*ArcSin[c*x])^2)/4 - (2*f*g*(1 - c^2*x^2)^(3/2) 
*(a + b*ArcSin[c*x])^2)/(3*c^2) + (f^2*(a + b*ArcSin[c*x])^3)/(6*b*c) + (g 
^2*(a + b*ArcSin[c*x])^3)/(24*b*c^3)))/Sqrt[1 - c^2*x^2]
 

Defintions of rubi rules used

rule 2009

Int[u_, x_Symbol] :> Simp[IntSum[u, x], x] /; SumQ[u]
 

rule 5262

Int[((a_.) + ArcSin[(c_.)*(x_)]*(b_.))^(n_.)*((f_) + (g_.)*(x_))^(m_.)*((d_ 
) + (e_.)*(x_)^2)^(p_), x_Symbol] :> Int[ExpandIntegrand[(d + e*x^2)^p*(a + 
 b*ArcSin[c*x])^n, (f + g*x)^m, x], x] /; FreeQ[{a, b, c, d, e, f, g}, x] & 
& EqQ[c^2*d + e, 0] && IGtQ[m, 0] && IntegerQ[p + 1/2] && GtQ[d, 0] && IGtQ 
[n, 0] && (m == 1 || p > 0 || (n == 1 && p > -1) || (m == 2 && p < -2))
 

rule 5276

Int[((a_.) + ArcSin[(c_.)*(x_)]*(b_.))^(n_.)*((f_) + (g_.)*(x_))^(m_.)*((d_ 
) + (e_.)*(x_)^2)^(p_), x_Symbol] :> Simp[Simp[(d + e*x^2)^p/(1 - c^2*x^2)^ 
p]   Int[(f + g*x)^m*(1 - c^2*x^2)^p*(a + b*ArcSin[c*x])^n, x], x] /; FreeQ 
[{a, b, c, d, e, f, g, n}, x] && EqQ[c^2*d + e, 0] && IntegerQ[m] && Intege 
rQ[p - 1/2] &&  !GtQ[d, 0]
 

Maple [C] (verified)

Result contains complex when optimal does not.

Time = 0.78 (sec) , antiderivative size = 1852, normalized size of antiderivative = 2.56

method	result	size
default	\(\text {Expression too large to display}\)	\(1852\)
parts	\(\text {Expression too large to display}\)	\(1852\)

Input:

int((g*x+f)^2*(-c^2*d*x^2+d)^(1/2)*(a+b*arcsin(c*x))^2,x,method=_RETURNVER 
BOSE)
 

Output:

a^2*(f^2*(1/2*x*(-c^2*d*x^2+d)^(1/2)+1/2*d/(c^2*d)^(1/2)*arctan((c^2*d)^(1 
/2)*x/(-c^2*d*x^2+d)^(1/2)))+g^2*(-1/4*x*(-c^2*d*x^2+d)^(3/2)/c^2/d+1/4/c^ 
2*(1/2*x*(-c^2*d*x^2+d)^(1/2)+1/2*d/(c^2*d)^(1/2)*arctan((c^2*d)^(1/2)*x/( 
-c^2*d*x^2+d)^(1/2))))-2/3*f*g/c^2/d*(-c^2*d*x^2+d)^(3/2))+b^2*(-1/24*(-d* 
(c^2*x^2-1))^(1/2)*(-c^2*x^2+1)^(1/2)/c^3/(c^2*x^2-1)*arcsin(c*x)^3*(4*c^2 
*f^2+g^2)+1/512*(-d*(c^2*x^2-1))^(1/2)*(-8*I*(-c^2*x^2+1)^(1/2)*x^4*c^4+8* 
c^5*x^5+8*I*(-c^2*x^2+1)^(1/2)*x^2*c^2-12*c^3*x^3-I*(-c^2*x^2+1)^(1/2)+4*c 
*x)*g^2*(4*I*arcsin(c*x)+8*arcsin(c*x)^2-1)/c^3/(c^2*x^2-1)+1/108*(-d*(c^2 
*x^2-1))^(1/2)*(4*c^4*x^4-5*c^2*x^2-4*I*(-c^2*x^2+1)^(1/2)*x^3*c^3+3*I*(-c 
^2*x^2+1)^(1/2)*x*c+1)*f*g*(6*I*arcsin(c*x)+9*arcsin(c*x)^2-2)/c^2/(c^2*x^ 
2-1)+1/16*(-d*(c^2*x^2-1))^(1/2)*(-2*I*(-c^2*x^2+1)^(1/2)*x^2*c^2+2*c^3*x^ 
3+I*(-c^2*x^2+1)^(1/2)-2*c*x)*f^2*(2*I*arcsin(c*x)+2*arcsin(c*x)^2-1)/c/(c 
^2*x^2-1)-1/4*(-d*(c^2*x^2-1))^(1/2)*(c^2*x^2-I*c*x*(-c^2*x^2+1)^(1/2)-1)* 
f*g*(arcsin(c*x)^2-2+2*I*arcsin(c*x))/c^2/(c^2*x^2-1)-1/4*(-d*(c^2*x^2-1)) 
^(1/2)*(I*(-c^2*x^2+1)^(1/2)*c*x+c^2*x^2-1)*f*g*(arcsin(c*x)^2-2-2*I*arcsi 
n(c*x))/c^2/(c^2*x^2-1)+1/16*(-d*(c^2*x^2-1))^(1/2)*(2*I*(-c^2*x^2+1)^(1/2 
)*x^2*c^2+2*c^3*x^3-I*(-c^2*x^2+1)^(1/2)-2*c*x)*f^2*(2*arcsin(c*x)^2-1-2*I 
*arcsin(c*x))/c/(c^2*x^2-1)+1/108*(-d*(c^2*x^2-1))^(1/2)*(4*I*c^3*x^3*(-c^ 
2*x^2+1)^(1/2)+4*c^4*x^4-3*I*(-c^2*x^2+1)^(1/2)*x*c-5*c^2*x^2+1)*f*g*(-6*I 
*arcsin(c*x)+9*arcsin(c*x)^2-2)/c^2/(c^2*x^2-1)+1/512*(-d*(c^2*x^2-1))^...
 

Fricas [F]

\[ \int (f+g x)^2 \sqrt {d-c^2 d x^2} (a+b \arcsin (c x))^2 \, dx=\int { \sqrt {-c^{2} d x^{2} + d} {\left (g x + f\right )}^{2} {\left (b \arcsin \left (c x\right ) + a\right )}^{2} \,d x } \] Input:

integrate((g*x+f)^2*(-c^2*d*x^2+d)^(1/2)*(a+b*arcsin(c*x))^2,x, algorithm= 
"fricas")
 

Output:

integral((a^2*g^2*x^2 + 2*a^2*f*g*x + a^2*f^2 + (b^2*g^2*x^2 + 2*b^2*f*g*x 
 + b^2*f^2)*arcsin(c*x)^2 + 2*(a*b*g^2*x^2 + 2*a*b*f*g*x + a*b*f^2)*arcsin 
(c*x))*sqrt(-c^2*d*x^2 + d), x)
 

Sympy [F]

\[ \int (f+g x)^2 \sqrt {d-c^2 d x^2} (a+b \arcsin (c x))^2 \, dx=\int \sqrt {- d \left (c x - 1\right ) \left (c x + 1\right )} \left (a + b \operatorname {asin}{\left (c x \right )}\right )^{2} \left (f + g x\right )^{2}\, dx \] Input:

integrate((g*x+f)**2*(-c**2*d*x**2+d)**(1/2)*(a+b*asin(c*x))**2,x)
 

Output:

Integral(sqrt(-d*(c*x - 1)*(c*x + 1))*(a + b*asin(c*x))**2*(f + g*x)**2, x 
)
 

Maxima [F]

\[ \int (f+g x)^2 \sqrt {d-c^2 d x^2} (a+b \arcsin (c x))^2 \, dx=\int { \sqrt {-c^{2} d x^{2} + d} {\left (g x + f\right )}^{2} {\left (b \arcsin \left (c x\right ) + a\right )}^{2} \,d x } \] Input:

integrate((g*x+f)^2*(-c^2*d*x^2+d)^(1/2)*(a+b*arcsin(c*x))^2,x, algorithm= 
"maxima")
 

Output:

1/2*(sqrt(-c^2*d*x^2 + d)*x + sqrt(d)*arcsin(c*x)/c)*a^2*f^2 + 1/8*a^2*g^2 
*(sqrt(-c^2*d*x^2 + d)*x/c^2 - 2*(-c^2*d*x^2 + d)^(3/2)*x/(c^2*d) + sqrt(d 
)*arcsin(c*x)/c^3) - 2/3*(-c^2*d*x^2 + d)^(3/2)*a^2*f*g/(c^2*d) + sqrt(d)* 
integrate(((b^2*g^2*x^2 + 2*b^2*f*g*x + b^2*f^2)*arctan2(c*x, sqrt(c*x + 1 
)*sqrt(-c*x + 1))^2 + 2*(a*b*g^2*x^2 + 2*a*b*f*g*x + a*b*f^2)*arctan2(c*x, 
 sqrt(c*x + 1)*sqrt(-c*x + 1)))*sqrt(c*x + 1)*sqrt(-c*x + 1), x)
 

Giac [F(-2)]

Exception generated. \[ \int (f+g x)^2 \sqrt {d-c^2 d x^2} (a+b \arcsin (c x))^2 \, dx=\text {Exception raised: RuntimeError} \] Input:

integrate((g*x+f)^2*(-c^2*d*x^2+d)^(1/2)*(a+b*arcsin(c*x))^2,x, algorithm= 
"giac")
 

Output:

Exception raised: RuntimeError >> an error occurred running a Giac command 
:INPUT:sage2OUTPUT:sym2poly/r2sym(const gen & e,const index_m & i,const ve 
cteur & l) Error: Bad Argument Value
 

Mupad [F(-1)]

Timed out. \[ \int (f+g x)^2 \sqrt {d-c^2 d x^2} (a+b \arcsin (c x))^2 \, dx=\int {\left (f+g\,x\right )}^2\,{\left (a+b\,\mathrm {asin}\left (c\,x\right )\right )}^2\,\sqrt {d-c^2\,d\,x^2} \,d x \] Input:

int((f + g*x)^2*(a + b*asin(c*x))^2*(d - c^2*d*x^2)^(1/2),x)
 

Output:

int((f + g*x)^2*(a + b*asin(c*x))^2*(d - c^2*d*x^2)^(1/2), x)
 

Reduce [F]

\[ \int (f+g x)^2 \sqrt {d-c^2 d x^2} (a+b \arcsin (c x))^2 \, dx=\frac {\sqrt {d}\, \left (12 \mathit {asin} \left (c x \right ) a^{2} c^{2} f^{2}+3 \mathit {asin} \left (c x \right ) a^{2} g^{2}+12 \sqrt {-c^{2} x^{2}+1}\, a^{2} c^{3} f^{2} x +16 \sqrt {-c^{2} x^{2}+1}\, a^{2} c^{3} f g \,x^{2}+6 \sqrt {-c^{2} x^{2}+1}\, a^{2} c^{3} g^{2} x^{3}-16 \sqrt {-c^{2} x^{2}+1}\, a^{2} c f g -3 \sqrt {-c^{2} x^{2}+1}\, a^{2} c \,g^{2} x +48 \left (\int \sqrt {-c^{2} x^{2}+1}\, \mathit {asin} \left (c x \right ) x^{2}d x \right ) a b \,c^{3} g^{2}+96 \left (\int \sqrt {-c^{2} x^{2}+1}\, \mathit {asin} \left (c x \right ) x d x \right ) a b \,c^{3} f g +48 \left (\int \sqrt {-c^{2} x^{2}+1}\, \mathit {asin} \left (c x \right )d x \right ) a b \,c^{3} f^{2}+24 \left (\int \sqrt {-c^{2} x^{2}+1}\, \mathit {asin} \left (c x \right )^{2} x^{2}d x \right ) b^{2} c^{3} g^{2}+48 \left (\int \sqrt {-c^{2} x^{2}+1}\, \mathit {asin} \left (c x \right )^{2} x d x \right ) b^{2} c^{3} f g +24 \left (\int \sqrt {-c^{2} x^{2}+1}\, \mathit {asin} \left (c x \right )^{2}d x \right ) b^{2} c^{3} f^{2}+16 a^{2} c f g \right )}{24 c^{3}} \] Input:

int((g*x+f)^2*(-c^2*d*x^2+d)^(1/2)*(a+b*asin(c*x))^2,x)
 

Output:

(sqrt(d)*(12*asin(c*x)*a**2*c**2*f**2 + 3*asin(c*x)*a**2*g**2 + 12*sqrt( - 
 c**2*x**2 + 1)*a**2*c**3*f**2*x + 16*sqrt( - c**2*x**2 + 1)*a**2*c**3*f*g 
*x**2 + 6*sqrt( - c**2*x**2 + 1)*a**2*c**3*g**2*x**3 - 16*sqrt( - c**2*x** 
2 + 1)*a**2*c*f*g - 3*sqrt( - c**2*x**2 + 1)*a**2*c*g**2*x + 48*int(sqrt( 
- c**2*x**2 + 1)*asin(c*x)*x**2,x)*a*b*c**3*g**2 + 96*int(sqrt( - c**2*x** 
2 + 1)*asin(c*x)*x,x)*a*b*c**3*f*g + 48*int(sqrt( - c**2*x**2 + 1)*asin(c* 
x),x)*a*b*c**3*f**2 + 24*int(sqrt( - c**2*x**2 + 1)*asin(c*x)**2*x**2,x)*b 
**2*c**3*g**2 + 48*int(sqrt( - c**2*x**2 + 1)*asin(c*x)**2*x,x)*b**2*c**3* 
f*g + 24*int(sqrt( - c**2*x**2 + 1)*asin(c*x)**2,x)*b**2*c**3*f**2 + 16*a* 
*2*c*f*g))/(24*c**3)