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

Optimal result

Integrand size = 26, antiderivative size = 163 \[ \int \frac {\sqrt {d-c^2 d x^2}}{(a+b \arcsin (c x))^2} \, dx=-\frac {\sqrt {1-c^2 x^2} \sqrt {d-c^2 d x^2}}{b c (a+b \arcsin (c x))}+\frac {\sqrt {d-c^2 d x^2} \operatorname {CosIntegral}\left (\frac {2 (a+b \arcsin (c x))}{b}\right ) \sin \left (\frac {2 a}{b}\right )}{b^2 c \sqrt {1-c^2 x^2}}-\frac {\sqrt {d-c^2 d x^2} \cos \left (\frac {2 a}{b}\right ) \text {Si}\left (\frac {2 (a+b \arcsin (c x))}{b}\right )}{b^2 c \sqrt {1-c^2 x^2}} \] Output:

-(-c^2*x^2+1)^(1/2)*(-c^2*d*x^2+d)^(1/2)/b/c/(a+b*arcsin(c*x))+(-c^2*d*x^2 
+d)^(1/2)*Ci(2*(a+b*arcsin(c*x))/b)*sin(2*a/b)/b^2/c/(-c^2*x^2+1)^(1/2)-(- 
c^2*d*x^2+d)^(1/2)*cos(2*a/b)*Si(2*(a+b*arcsin(c*x))/b)/b^2/c/(-c^2*x^2+1) 
^(1/2)
 

Mathematica [A] (verified)

Time = 0.35 (sec) , antiderivative size = 118, normalized size of antiderivative = 0.72 \[ \int \frac {\sqrt {d-c^2 d x^2}}{(a+b \arcsin (c x))^2} \, dx=\frac {d \sqrt {1-c^2 x^2} \left (b \left (-1+c^2 x^2\right )+(a+b \arcsin (c x)) \operatorname {CosIntegral}\left (2 \left (\frac {a}{b}+\arcsin (c x)\right )\right ) \sin \left (\frac {2 a}{b}\right )-(a+b \arcsin (c x)) \cos \left (\frac {2 a}{b}\right ) \text {Si}\left (2 \left (\frac {a}{b}+\arcsin (c x)\right )\right )\right )}{b^2 c \sqrt {d-c^2 d x^2} (a+b \arcsin (c x))} \] Input:

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

Output:

(d*Sqrt[1 - c^2*x^2]*(b*(-1 + c^2*x^2) + (a + b*ArcSin[c*x])*CosIntegral[2 
*(a/b + ArcSin[c*x])]*Sin[(2*a)/b] - (a + b*ArcSin[c*x])*Cos[(2*a)/b]*SinI 
ntegral[2*(a/b + ArcSin[c*x])]))/(b^2*c*Sqrt[d - c^2*d*x^2]*(a + b*ArcSin[ 
c*x]))
 

Rubi [A] (verified)

Time = 0.83 (sec) , antiderivative size = 131, normalized size of antiderivative = 0.80, number of steps used = 12, number of rules used = 11, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.423, Rules used = {5166, 5146, 25, 4906, 27, 3042, 3784, 25, 3042, 3780, 3783}

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

\(\Big \downarrow \) 5166

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

\(\Big \downarrow \) 5146

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

\(\Big \downarrow \) 25

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

\(\Big \downarrow \) 4906

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

\(\Big \downarrow \) 27

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

\(\Big \downarrow \) 3042

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

\(\Big \downarrow \) 3784

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

\(\Big \downarrow \) 25

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

\(\Big \downarrow \) 3042

\(\displaystyle -\frac {\sqrt {d-c^2 d x^2} \left (\cos \left (\frac {2 a}{b}\right ) \int \frac {\sin \left (\frac {2 (a+b \arcsin (c x))}{b}\right )}{a+b \arcsin (c x)}d(a+b \arcsin (c x))-\sin \left (\frac {2 a}{b}\right ) \int \frac {\sin \left (\frac {2 (a+b \arcsin (c x))}{b}+\frac {\pi }{2}\right )}{a+b \arcsin (c x)}d(a+b \arcsin (c x))\right )}{b^2 c \sqrt {1-c^2 x^2}}-\frac {\sqrt {1-c^2 x^2} \sqrt {d-c^2 d x^2}}{b c (a+b \arcsin (c x))}\)

\(\Big \downarrow \) 3780

\(\displaystyle -\frac {\sqrt {d-c^2 d x^2} \left (\cos \left (\frac {2 a}{b}\right ) \text {Si}\left (\frac {2 (a+b \arcsin (c x))}{b}\right )-\sin \left (\frac {2 a}{b}\right ) \int \frac {\sin \left (\frac {2 (a+b \arcsin (c x))}{b}+\frac {\pi }{2}\right )}{a+b \arcsin (c x)}d(a+b \arcsin (c x))\right )}{b^2 c \sqrt {1-c^2 x^2}}-\frac {\sqrt {1-c^2 x^2} \sqrt {d-c^2 d x^2}}{b c (a+b \arcsin (c x))}\)

\(\Big \downarrow \) 3783

\(\displaystyle -\frac {\sqrt {d-c^2 d x^2} \left (\cos \left (\frac {2 a}{b}\right ) \text {Si}\left (\frac {2 (a+b \arcsin (c x))}{b}\right )-\sin \left (\frac {2 a}{b}\right ) \operatorname {CosIntegral}\left (\frac {2 (a+b \arcsin (c x))}{b}\right )\right )}{b^2 c \sqrt {1-c^2 x^2}}-\frac {\sqrt {1-c^2 x^2} \sqrt {d-c^2 d x^2}}{b c (a+b \arcsin (c x))}\)

Input:

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

Output:

-((Sqrt[1 - c^2*x^2]*Sqrt[d - c^2*d*x^2])/(b*c*(a + b*ArcSin[c*x]))) - (Sq 
rt[d - c^2*d*x^2]*(-(CosIntegral[(2*(a + b*ArcSin[c*x]))/b]*Sin[(2*a)/b]) 
+ Cos[(2*a)/b]*SinIntegral[(2*(a + b*ArcSin[c*x]))/b]))/(b^2*c*Sqrt[1 - c^ 
2*x^2])
 

Defintions of rubi rules used

rule 25

Int[-(Fx_), x_Symbol] :> Simp[Identity[-1]   Int[Fx, x], x]
 

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 3042

Int[u_, x_Symbol] :> Int[DeactivateTrig[u, x], x] /; FunctionOfTrigOfLinear 
Q[u, x]
 

rule 3780

Int[sin[(e_.) + (f_.)*(x_)]/((c_.) + (d_.)*(x_)), x_Symbol] :> Simp[SinInte 
gral[e + f*x]/d, x] /; FreeQ[{c, d, e, f}, x] && EqQ[d*e - c*f, 0]
 

rule 3783

Int[sin[(e_.) + (f_.)*(x_)]/((c_.) + (d_.)*(x_)), x_Symbol] :> Simp[CosInte 
gral[e - Pi/2 + f*x]/d, x] /; FreeQ[{c, d, e, f}, x] && EqQ[d*(e - Pi/2) - 
c*f, 0]
 

rule 3784

Int[sin[(e_.) + (f_.)*(x_)]/((c_.) + (d_.)*(x_)), x_Symbol] :> Simp[Cos[(d* 
e - c*f)/d]   Int[Sin[c*(f/d) + f*x]/(c + d*x), x], x] + Simp[Sin[(d*e - c* 
f)/d]   Int[Cos[c*(f/d) + f*x]/(c + d*x), x], x] /; FreeQ[{c, d, e, f}, x] 
&& NeQ[d*e - c*f, 0]
 

rule 4906

Int[Cos[(a_.) + (b_.)*(x_)]^(p_.)*((c_.) + (d_.)*(x_))^(m_.)*Sin[(a_.) + (b 
_.)*(x_)]^(n_.), x_Symbol] :> Int[ExpandTrigReduce[(c + d*x)^m, Sin[a + b*x 
]^n*Cos[a + b*x]^p, x], x] /; FreeQ[{a, b, c, d, m}, x] && IGtQ[n, 0] && IG 
tQ[p, 0]
 

rule 5146

Int[((a_.) + ArcSin[(c_.)*(x_)]*(b_.))^(n_)*(x_)^(m_.), x_Symbol] :> Simp[1 
/(b*c^(m + 1))   Subst[Int[x^n*Sin[-a/b + x/b]^m*Cos[-a/b + x/b], x], x, a 
+ b*ArcSin[c*x]], x] /; FreeQ[{a, b, c, n}, x] && IGtQ[m, 0]
 

rule 5166

Int[((a_.) + ArcSin[(c_.)*(x_)]*(b_.))^(n_)*((d_) + (e_.)*(x_)^2)^(p_.), x_ 
Symbol] :> Simp[Sqrt[1 - c^2*x^2]*(d + e*x^2)^p*((a + b*ArcSin[c*x])^(n + 1 
)/(b*c*(n + 1))), x] + Simp[c*((2*p + 1)/(b*(n + 1)))*Simp[(d + e*x^2)^p/(1 
 - c^2*x^2)^p]   Int[x*(1 - c^2*x^2)^(p - 1/2)*(a + b*ArcSin[c*x])^(n + 1), 
 x], x] /; FreeQ[{a, b, c, d, e, p}, x] && EqQ[c^2*d + e, 0] && LtQ[n, -1]
 

Maple [C] (verified)

Result contains complex when optimal does not.

Time = 0.16 (sec) , antiderivative size = 281, normalized size of antiderivative = 1.72

method	result	size
default	\(-\frac {\sqrt {-d \left (c^{2} x^{2}-1\right )}\, \left (i c^{2} x^{2}-c x \sqrt {-c^{2} x^{2}+1}-i\right ) \left (2 i \sqrt {-c^{2} x^{2}+1}\, b \,c^{2} x^{2}-2 b \,c^{3} x^{3}+\operatorname {expIntegral}_{1}\left (-2 i \arcsin \left (c x \right )-\frac {2 i a}{b}\right ) {\mathrm e}^{-\frac {i \left (-b \arcsin \left (c x \right )+2 a \right )}{b}} b \arcsin \left (c x \right )-\operatorname {expIntegral}_{1}\left (2 i \arcsin \left (c x \right )+\frac {2 i a}{b}\right ) {\mathrm e}^{\frac {i \left (b \arcsin \left (c x \right )+2 a \right )}{b}} b \arcsin \left (c x \right )-2 i \sqrt {-c^{2} x^{2}+1}\, b +\operatorname {expIntegral}_{1}\left (-2 i \arcsin \left (c x \right )-\frac {2 i a}{b}\right ) {\mathrm e}^{-\frac {i \left (-b \arcsin \left (c x \right )+2 a \right )}{b}} a -\operatorname {expIntegral}_{1}\left (2 i \arcsin \left (c x \right )+\frac {2 i a}{b}\right ) {\mathrm e}^{\frac {i \left (b \arcsin \left (c x \right )+2 a \right )}{b}} a +2 x b c \right )}{2 c \left (c^{2} x^{2}-1\right ) b^{2} \left (a +b \arcsin \left (c x \right )\right )}\)	\(281\)

Input:

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

Output:

-1/2*(-d*(c^2*x^2-1))^(1/2)*(I*c^2*x^2-c*x*(-c^2*x^2+1)^(1/2)-I)*(2*I*(-c^ 
2*x^2+1)^(1/2)*b*c^2*x^2-2*b*c^3*x^3+Ei(1,-2*I*arcsin(c*x)-2*I*a/b)*exp(-I 
*(-b*arcsin(c*x)+2*a)/b)*b*arcsin(c*x)-Ei(1,2*I*arcsin(c*x)+2*I*a/b)*exp(I 
*(b*arcsin(c*x)+2*a)/b)*b*arcsin(c*x)-2*I*(-c^2*x^2+1)^(1/2)*b+Ei(1,-2*I*a 
rcsin(c*x)-2*I*a/b)*exp(-I*(-b*arcsin(c*x)+2*a)/b)*a-Ei(1,2*I*arcsin(c*x)+ 
2*I*a/b)*exp(I*(b*arcsin(c*x)+2*a)/b)*a+2*x*b*c)/c/(c^2*x^2-1)/b^2/(a+b*ar 
csin(c*x))
 

Fricas [F]

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

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

Output:

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

Sympy [F]

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

integrate((-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, x)
 

Maxima [F]

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

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

Output:

(c^2*x^2 - 2*(b^2*c^2*arctan2(c*x, sqrt(c*x + 1)*sqrt(-c*x + 1)) + a*b*c^2 
)*integrate(x/(b^2*arctan2(c*x, sqrt(c*x + 1)*sqrt(-c*x + 1)) + a*b), x) - 
 1)*sqrt(d)/(b^2*c*arctan2(c*x, sqrt(c*x + 1)*sqrt(-c*x + 1)) + a*b*c)
 

Giac [B] (verification not implemented)

Leaf count of result is larger than twice the leaf count of optimal. 323 vs. \(2 (151) = 302\).

Time = 0.69 (sec) , antiderivative size = 323, normalized size of antiderivative = 1.98 \[ \int \frac {\sqrt {d-c^2 d x^2}}{(a+b \arcsin (c x))^2} \, dx={\left (\frac {2 \, b \arcsin \left (c x\right ) \cos \left (\frac {a}{b}\right ) \operatorname {Ci}\left (\frac {2 \, a}{b} + 2 \, \arcsin \left (c x\right )\right ) \sin \left (\frac {a}{b}\right )}{b^{3} c^{2} \arcsin \left (c x\right ) + a b^{2} c^{2}} - \frac {2 \, b \arcsin \left (c x\right ) \cos \left (\frac {a}{b}\right )^{2} \operatorname {Si}\left (\frac {2 \, a}{b} + 2 \, \arcsin \left (c x\right )\right )}{b^{3} c^{2} \arcsin \left (c x\right ) + a b^{2} c^{2}} + \frac {2 \, a \cos \left (\frac {a}{b}\right ) \operatorname {Ci}\left (\frac {2 \, a}{b} + 2 \, \arcsin \left (c x\right )\right ) \sin \left (\frac {a}{b}\right )}{b^{3} c^{2} \arcsin \left (c x\right ) + a b^{2} c^{2}} - \frac {2 \, a \cos \left (\frac {a}{b}\right )^{2} \operatorname {Si}\left (\frac {2 \, a}{b} + 2 \, \arcsin \left (c x\right )\right )}{b^{3} c^{2} \arcsin \left (c x\right ) + a b^{2} c^{2}} + \frac {b \arcsin \left (c x\right ) \operatorname {Si}\left (\frac {2 \, a}{b} + 2 \, \arcsin \left (c x\right )\right )}{b^{3} c^{2} \arcsin \left (c x\right ) + a b^{2} c^{2}} + \frac {{\left (c^{2} x^{2} - 1\right )} b}{b^{3} c^{2} \arcsin \left (c x\right ) + a b^{2} c^{2}} + \frac {a \operatorname {Si}\left (\frac {2 \, a}{b} + 2 \, \arcsin \left (c x\right )\right )}{b^{3} c^{2} \arcsin \left (c x\right ) + a b^{2} c^{2}}\right )} c \sqrt {d} \] Input:

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

Output:

(2*b*arcsin(c*x)*cos(a/b)*cos_integral(2*a/b + 2*arcsin(c*x))*sin(a/b)/(b^ 
3*c^2*arcsin(c*x) + a*b^2*c^2) - 2*b*arcsin(c*x)*cos(a/b)^2*sin_integral(2 
*a/b + 2*arcsin(c*x))/(b^3*c^2*arcsin(c*x) + a*b^2*c^2) + 2*a*cos(a/b)*cos 
_integral(2*a/b + 2*arcsin(c*x))*sin(a/b)/(b^3*c^2*arcsin(c*x) + a*b^2*c^2 
) - 2*a*cos(a/b)^2*sin_integral(2*a/b + 2*arcsin(c*x))/(b^3*c^2*arcsin(c*x 
) + a*b^2*c^2) + b*arcsin(c*x)*sin_integral(2*a/b + 2*arcsin(c*x))/(b^3*c^ 
2*arcsin(c*x) + a*b^2*c^2) + (c^2*x^2 - 1)*b/(b^3*c^2*arcsin(c*x) + a*b^2* 
c^2) + a*sin_integral(2*a/b + 2*arcsin(c*x))/(b^3*c^2*arcsin(c*x) + a*b^2* 
c^2))*c*sqrt(d)
 

Mupad [F(-1)]

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

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

Output:

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

Reduce [F]

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

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

Output:

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