\(\int \frac {\sqrt {a-a \sin (e+f x)} \sqrt {c+c \sin (e+f x)}}{x} \, dx\) [168]

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

Optimal result

Integrand size = 33, antiderivative size = 86 \[ \int \frac {\sqrt {a-a \sin (e+f x)} \sqrt {c+c \sin (e+f x)}}{x} \, dx=\cos (e) \operatorname {CosIntegral}(f x) \sec (e+f x) \sqrt {a-a \sin (e+f x)} \sqrt {c+c \sin (e+f x)}-\sec (e+f x) \sin (e) \sqrt {a-a \sin (e+f x)} \sqrt {c+c \sin (e+f x)} \text {Si}(f x) \]

[Out]

Ci(f*x)*cos(e)*sec(f*x+e)*(a-a*sin(f*x+e))^(1/2)*(c+c*sin(f*x+e))^(1/2)-sec(f*x+e)*Si(f*x)*sin(e)*(a-a*sin(f*x
+e))^(1/2)*(c+c*sin(f*x+e))^(1/2)

Rubi [A] (verified)

Time = 0.22 (sec) , antiderivative size = 86, normalized size of antiderivative = 1.00, number of steps used = 4, number of rules used = 4, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.121, Rules used = {4700, 3384, 3380, 3383} \[ \int \frac {\sqrt {a-a \sin (e+f x)} \sqrt {c+c \sin (e+f x)}}{x} \, dx=\cos (e) \operatorname {CosIntegral}(f x) \sec (e+f x) \sqrt {a-a \sin (e+f x)} \sqrt {c \sin (e+f x)+c}-\sin (e) \text {Si}(f x) \sec (e+f x) \sqrt {a-a \sin (e+f x)} \sqrt {c \sin (e+f x)+c} \]

[In]

Int[(Sqrt[a - a*Sin[e + f*x]]*Sqrt[c + c*Sin[e + f*x]])/x,x]

[Out]

Cos[e]*CosIntegral[f*x]*Sec[e + f*x]*Sqrt[a - a*Sin[e + f*x]]*Sqrt[c + c*Sin[e + f*x]] - Sec[e + f*x]*Sin[e]*S
qrt[a - a*Sin[e + f*x]]*Sqrt[c + c*Sin[e + f*x]]*SinIntegral[f*x]

Rule 3380

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

Rule 3383

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

Rule 3384

Int[sin[(e_.) + (f_.)*(x_)]/((c_.) + (d_.)*(x_)), x_Symbol] :> Dist[Cos[(d*e - c*f)/d], Int[Sin[c*(f/d) + f*x]
/(c + d*x), x], x] + Dist[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 4700

Int[((g_.) + (h_.)*(x_))^(p_.)*((a_) + (b_.)*Sin[(e_.) + (f_.)*(x_)])^(m_)*((c_) + (d_.)*Sin[(e_.) + (f_.)*(x_
)])^(n_), x_Symbol] :> Dist[a^IntPart[m]*c^IntPart[m]*(a + b*Sin[e + f*x])^FracPart[m]*((c + d*Sin[e + f*x])^F
racPart[m]/Cos[e + f*x]^(2*FracPart[m])), Int[(g + h*x)^p*Cos[e + f*x]^(2*m)*(c + d*Sin[e + f*x])^(n - m), x],
 x] /; FreeQ[{a, b, c, d, e, f, g, h}, x] && EqQ[b*c + a*d, 0] && EqQ[a^2 - b^2, 0] && IntegerQ[p] && IntegerQ
[2*m] && IGeQ[n - m, 0]

Rubi steps \begin{align*} \text {integral}& = \left (\sec (e+f x) \sqrt {a-a \sin (e+f x)} \sqrt {c+c \sin (e+f x)}\right ) \int \frac {\cos (e+f x)}{x} \, dx \\ & = \left (\cos (e) \sec (e+f x) \sqrt {a-a \sin (e+f x)} \sqrt {c+c \sin (e+f x)}\right ) \int \frac {\cos (f x)}{x} \, dx-\left (\sec (e+f x) \sin (e) \sqrt {a-a \sin (e+f x)} \sqrt {c+c \sin (e+f x)}\right ) \int \frac {\sin (f x)}{x} \, dx \\ & = \cos (e) \operatorname {CosIntegral}(f x) \sec (e+f x) \sqrt {a-a \sin (e+f x)} \sqrt {c+c \sin (e+f x)}-\sec (e+f x) \sin (e) \sqrt {a-a \sin (e+f x)} \sqrt {c+c \sin (e+f x)} \text {Si}(f x) \\ \end{align*}

Mathematica [A] (verified)

Time = 0.76 (sec) , antiderivative size = 52, normalized size of antiderivative = 0.60 \[ \int \frac {\sqrt {a-a \sin (e+f x)} \sqrt {c+c \sin (e+f x)}}{x} \, dx=\sec (e+f x) \sqrt {c (1+\sin (e+f x))} \sqrt {a-a \sin (e+f x)} (\cos (e) \operatorname {CosIntegral}(f x)-\sin (e) \text {Si}(f x)) \]

[In]

Integrate[(Sqrt[a - a*Sin[e + f*x]]*Sqrt[c + c*Sin[e + f*x]])/x,x]

[Out]

Sec[e + f*x]*Sqrt[c*(1 + Sin[e + f*x])]*Sqrt[a - a*Sin[e + f*x]]*(Cos[e]*CosIntegral[f*x] - Sin[e]*SinIntegral
[f*x])

Maple [F]

\[\int \frac {\sqrt {a -\sin \left (f x +e \right ) a}\, \sqrt {c +c \sin \left (f x +e \right )}}{x}d x\]

[In]

int((a-sin(f*x+e)*a)^(1/2)*(c+c*sin(f*x+e))^(1/2)/x,x)

[Out]

int((a-sin(f*x+e)*a)^(1/2)*(c+c*sin(f*x+e))^(1/2)/x,x)

Fricas [F(-2)]

Exception generated. \[ \int \frac {\sqrt {a-a \sin (e+f x)} \sqrt {c+c \sin (e+f x)}}{x} \, dx=\text {Exception raised: TypeError} \]

[In]

integrate((a-a*sin(f*x+e))^(1/2)*(c+c*sin(f*x+e))^(1/2)/x,x, algorithm="fricas")

[Out]

Exception raised: TypeError >>  Error detected within library code:   integrate: implementation incomplete (ha
s polynomial part)

Sympy [F]

\[ \int \frac {\sqrt {a-a \sin (e+f x)} \sqrt {c+c \sin (e+f x)}}{x} \, dx=\int \frac {\sqrt {c \left (\sin {\left (e + f x \right )} + 1\right )} \sqrt {- a \left (\sin {\left (e + f x \right )} - 1\right )}}{x}\, dx \]

[In]

integrate((a-a*sin(f*x+e))**(1/2)*(c+c*sin(f*x+e))**(1/2)/x,x)

[Out]

Integral(sqrt(c*(sin(e + f*x) + 1))*sqrt(-a*(sin(e + f*x) - 1))/x, x)

Maxima [F]

\[ \int \frac {\sqrt {a-a \sin (e+f x)} \sqrt {c+c \sin (e+f x)}}{x} \, dx=\int { \frac {\sqrt {-a \sin \left (f x + e\right ) + a} \sqrt {c \sin \left (f x + e\right ) + c}}{x} \,d x } \]

[In]

integrate((a-a*sin(f*x+e))^(1/2)*(c+c*sin(f*x+e))^(1/2)/x,x, algorithm="maxima")

[Out]

integrate(sqrt(-a*sin(f*x + e) + a)*sqrt(c*sin(f*x + e) + c)/x, x)

Giac [C] (verification not implemented)

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

Time = 0.32 (sec) , antiderivative size = 274, normalized size of antiderivative = 3.19 \[ \int \frac {\sqrt {a-a \sin (e+f x)} \sqrt {c+c \sin (e+f x)}}{x} \, dx=\frac {{\left (\Re \left ( \operatorname {Ci}\left (f x\right ) \right ) \mathrm {sgn}\left (\cos \left (-\frac {1}{4} \, \pi + \frac {1}{2} \, f x + \frac {1}{2} \, e\right )\right ) \mathrm {sgn}\left (\sin \left (-\frac {1}{4} \, \pi + \frac {1}{2} \, f x + \frac {1}{2} \, e\right )\right ) \tan \left (\frac {1}{2} \, e\right )^{2} + \Re \left ( \operatorname {Ci}\left (-f x\right ) \right ) \mathrm {sgn}\left (\cos \left (-\frac {1}{4} \, \pi + \frac {1}{2} \, f x + \frac {1}{2} \, e\right )\right ) \mathrm {sgn}\left (\sin \left (-\frac {1}{4} \, \pi + \frac {1}{2} \, f x + \frac {1}{2} \, e\right )\right ) \tan \left (\frac {1}{2} \, e\right )^{2} + 2 \, \Im \left ( \operatorname {Ci}\left (f x\right ) \right ) \mathrm {sgn}\left (\cos \left (-\frac {1}{4} \, \pi + \frac {1}{2} \, f x + \frac {1}{2} \, e\right )\right ) \mathrm {sgn}\left (\sin \left (-\frac {1}{4} \, \pi + \frac {1}{2} \, f x + \frac {1}{2} \, e\right )\right ) \tan \left (\frac {1}{2} \, e\right ) - 2 \, \Im \left ( \operatorname {Ci}\left (-f x\right ) \right ) \mathrm {sgn}\left (\cos \left (-\frac {1}{4} \, \pi + \frac {1}{2} \, f x + \frac {1}{2} \, e\right )\right ) \mathrm {sgn}\left (\sin \left (-\frac {1}{4} \, \pi + \frac {1}{2} \, f x + \frac {1}{2} \, e\right )\right ) \tan \left (\frac {1}{2} \, e\right ) + 4 \, \mathrm {sgn}\left (\cos \left (-\frac {1}{4} \, \pi + \frac {1}{2} \, f x + \frac {1}{2} \, e\right )\right ) \mathrm {sgn}\left (\sin \left (-\frac {1}{4} \, \pi + \frac {1}{2} \, f x + \frac {1}{2} \, e\right )\right ) \operatorname {Si}\left (f x\right ) \tan \left (\frac {1}{2} \, e\right ) - \Re \left ( \operatorname {Ci}\left (f x\right ) \right ) \mathrm {sgn}\left (\cos \left (-\frac {1}{4} \, \pi + \frac {1}{2} \, f x + \frac {1}{2} \, e\right )\right ) \mathrm {sgn}\left (\sin \left (-\frac {1}{4} \, \pi + \frac {1}{2} \, f x + \frac {1}{2} \, e\right )\right ) - \Re \left ( \operatorname {Ci}\left (-f x\right ) \right ) \mathrm {sgn}\left (\cos \left (-\frac {1}{4} \, \pi + \frac {1}{2} \, f x + \frac {1}{2} \, e\right )\right ) \mathrm {sgn}\left (\sin \left (-\frac {1}{4} \, \pi + \frac {1}{2} \, f x + \frac {1}{2} \, e\right )\right )\right )} \sqrt {a} \sqrt {c}}{2 \, {\left (\tan \left (\frac {1}{2} \, e\right )^{2} + 1\right )}} \]

[In]

integrate((a-a*sin(f*x+e))^(1/2)*(c+c*sin(f*x+e))^(1/2)/x,x, algorithm="giac")

[Out]

1/2*(real_part(cos_integral(f*x))*sgn(cos(-1/4*pi + 1/2*f*x + 1/2*e))*sgn(sin(-1/4*pi + 1/2*f*x + 1/2*e))*tan(
1/2*e)^2 + real_part(cos_integral(-f*x))*sgn(cos(-1/4*pi + 1/2*f*x + 1/2*e))*sgn(sin(-1/4*pi + 1/2*f*x + 1/2*e
))*tan(1/2*e)^2 + 2*imag_part(cos_integral(f*x))*sgn(cos(-1/4*pi + 1/2*f*x + 1/2*e))*sgn(sin(-1/4*pi + 1/2*f*x
 + 1/2*e))*tan(1/2*e) - 2*imag_part(cos_integral(-f*x))*sgn(cos(-1/4*pi + 1/2*f*x + 1/2*e))*sgn(sin(-1/4*pi +
1/2*f*x + 1/2*e))*tan(1/2*e) + 4*sgn(cos(-1/4*pi + 1/2*f*x + 1/2*e))*sgn(sin(-1/4*pi + 1/2*f*x + 1/2*e))*sin_i
ntegral(f*x)*tan(1/2*e) - real_part(cos_integral(f*x))*sgn(cos(-1/4*pi + 1/2*f*x + 1/2*e))*sgn(sin(-1/4*pi + 1
/2*f*x + 1/2*e)) - real_part(cos_integral(-f*x))*sgn(cos(-1/4*pi + 1/2*f*x + 1/2*e))*sgn(sin(-1/4*pi + 1/2*f*x
 + 1/2*e)))*sqrt(a)*sqrt(c)/(tan(1/2*e)^2 + 1)

Mupad [F(-1)]

Timed out. \[ \int \frac {\sqrt {a-a \sin (e+f x)} \sqrt {c+c \sin (e+f x)}}{x} \, dx=\int \frac {\sqrt {a-a\,\sin \left (e+f\,x\right )}\,\sqrt {c+c\,\sin \left (e+f\,x\right )}}{x} \,d x \]

[In]

int(((a - a*sin(e + f*x))^(1/2)*(c + c*sin(e + f*x))^(1/2))/x,x)

[Out]

int(((a - a*sin(e + f*x))^(1/2)*(c + c*sin(e + f*x))^(1/2))/x, x)