3.140 \(\int \tan (c-b x) \tan (a+b x) \, dx\)
Optimal. Leaf size=34 \[ -\frac {\cot (a+c) \log (\cos (c-b x))}{b}+\frac {\cot (a+c) \log (\cos (a+b x))}{b}+x \]
[Out]
x-cot(a+c)*ln(cos(b*x-c))/b+cot(a+c)*ln(cos(b*x+a))/b
________________________________________________________________________________________
Rubi [A] time = 0.07, antiderivative size = 34, normalized size of antiderivative = 1.00,
number of steps used = 4, number of rules used = 3, integrand size = 14, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.214, Rules used =
{4612, 4610, 3475} \[ -\frac {\cot (a+c) \log (\cos (c-b x))}{b}+\frac {\cot (a+c) \log (\cos (a+b x))}{b}+x \]
Antiderivative was successfully verified.
[In]
Int[Tan[c - b*x]*Tan[a + b*x],x]
[Out]
x - (Cot[a + c]*Log[Cos[c - b*x]])/b + (Cot[a + c]*Log[Cos[a + b*x]])/b
Rule 3475
Int[tan[(c_.) + (d_.)*(x_)], x_Symbol] :> -Simp[Log[RemoveContent[Cos[c + d*x], x]]/d, x] /; FreeQ[{c, d}, x]
Rule 4610
Int[Sec[(a_.) + (b_.)*(x_)]*Sec[(c_) + (d_.)*(x_)], x_Symbol] :> -Dist[Csc[(b*c - a*d)/d], Int[Tan[a + b*x], x
], x] + Dist[Csc[(b*c - a*d)/b], Int[Tan[c + d*x], x], x] /; FreeQ[{a, b, c, d}, x] && EqQ[b^2 - d^2, 0] && Ne
Q[b*c - a*d, 0]
Rule 4612
Int[Tan[(a_.) + (b_.)*(x_)]*Tan[(c_) + (d_.)*(x_)], x_Symbol] :> -Simp[(b*x)/d, x] + Dist[(b*Cos[(b*c - a*d)/d
])/d, Int[Sec[a + b*x]*Sec[c + d*x], x], x] /; FreeQ[{a, b, c, d}, x] && EqQ[b^2 - d^2, 0] && NeQ[b*c - a*d, 0
]
Rubi steps
\begin {align*} \int \tan (c-b x) \tan (a+b x) \, dx &=x-\cos (a+c) \int \sec (c-b x) \sec (a+b x) \, dx\\ &=x-\cot (a+c) \int \tan (c-b x) \, dx-\cot (a+c) \int \tan (a+b x) \, dx\\ &=x-\frac {\cot (a+c) \log (\cos (c-b x))}{b}+\frac {\cot (a+c) \log (\cos (a+b x))}{b}\\ \end {align*}
________________________________________________________________________________________
Mathematica [A] time = 0.53, size = 28, normalized size = 0.82 \[ \frac {\cot (a+c) (\log (\cos (a+b x))-\log (\cos (c-b x)))}{b}+x \]
Antiderivative was successfully verified.
[In]
Integrate[Tan[c - b*x]*Tan[a + b*x],x]
[Out]
x + (Cot[a + c]*(-Log[Cos[c - b*x]] + Log[Cos[a + b*x]]))/b
________________________________________________________________________________________
fricas [B] time = 2.01, size = 145, normalized size = 4.26 \[ \frac {2 \, b x \sin \left (2 \, a + 2 \, c\right ) - {\left (\cos \left (2 \, a + 2 \, c\right ) + 1\right )} \log \left (-\frac {{\left (\cos \left (2 \, a + 2 \, c\right ) - 1\right )} \tan \left (b x + a\right )^{2} - 2 \, \sin \left (2 \, a + 2 \, c\right ) \tan \left (b x + a\right ) - \cos \left (2 \, a + 2 \, c\right ) - 1}{{\left (\cos \left (2 \, a + 2 \, c\right ) + 1\right )} \tan \left (b x + a\right )^{2} + \cos \left (2 \, a + 2 \, c\right ) + 1}\right ) + {\left (\cos \left (2 \, a + 2 \, c\right ) + 1\right )} \log \left (\frac {1}{\tan \left (b x + a\right )^{2} + 1}\right )}{2 \, b \sin \left (2 \, a + 2 \, c\right )} \]
Verification of antiderivative is not currently implemented for this CAS.
[In]
integrate(-tan(b*x-c)*tan(b*x+a),x, algorithm="fricas")
[Out]
1/2*(2*b*x*sin(2*a + 2*c) - (cos(2*a + 2*c) + 1)*log(-((cos(2*a + 2*c) - 1)*tan(b*x + a)^2 - 2*sin(2*a + 2*c)*
tan(b*x + a) - cos(2*a + 2*c) - 1)/((cos(2*a + 2*c) + 1)*tan(b*x + a)^2 + cos(2*a + 2*c) + 1)) + (cos(2*a + 2*
c) + 1)*log(1/(tan(b*x + a)^2 + 1)))/(b*sin(2*a + 2*c))
________________________________________________________________________________________
giac [F(-2)] time = 0.00, size = 0, normalized size = 0.00 \[ \text {Exception raised: NotImplementedError} \]
Verification of antiderivative is not currently implemented for this CAS.
[In]
integrate(-tan(b*x-c)*tan(b*x+a),x, algorithm="giac")
[Out]
Exception raised: NotImplementedError >> Unable to parse Giac output: Unable to check sign: (2*pi/x/2)>(-2*pi/
x/2)-2/b*(-1/2*b*x+(-2*tan(a/2)^2*tan(c/2)^3+2*tan(a/2)^2*tan(c/2)+8*tan(a/2)*tan(c/2)^2+2*tan(c/2)^3-2*tan(c/
2))/(8*tan(a/2)^2*tan(c/2)^2+8*tan(a/2)*tan(c/2)^3-8*tan(a/2)*tan(c/2)-8*tan(c/2)^2)*ln(abs(2*tan(b*x)*tan(c/2
)-tan(c/2)^2+1))+(-2*tan(a/2)^3*tan(c/2)^2+2*tan(a/2)^3+8*tan(a/2)^2*tan(c/2)+2*tan(a/2)*tan(c/2)^2-2*tan(a/2)
)/(-8*tan(a/2)^3*tan(c/2)-8*tan(a/2)^2*tan(c/2)^2+8*tan(a/2)^2+8*tan(a/2)*tan(c/2))*ln(abs(2*tan(b*x)*tan(a/2)
+tan(a/2)^2-1)))
________________________________________________________________________________________
maple [C] time = 0.14, size = 145, normalized size = 4.26 \[ x +\frac {i \ln \left (1+{\mathrm e}^{2 i \left (b x +a \right )}\right ) {\mathrm e}^{2 i \left (a +c \right )}}{b \left ({\mathrm e}^{2 i \left (a +c \right )}-1\right )}+\frac {i \ln \left (1+{\mathrm e}^{2 i \left (b x +a \right )}\right )}{b \left ({\mathrm e}^{2 i \left (a +c \right )}-1\right )}-\frac {i \ln \left ({\mathrm e}^{2 i \left (a +c \right )}+{\mathrm e}^{2 i \left (b x +a \right )}\right ) {\mathrm e}^{2 i \left (a +c \right )}}{b \left ({\mathrm e}^{2 i \left (a +c \right )}-1\right )}-\frac {i \ln \left ({\mathrm e}^{2 i \left (a +c \right )}+{\mathrm e}^{2 i \left (b x +a \right )}\right )}{b \left ({\mathrm e}^{2 i \left (a +c \right )}-1\right )} \]
Verification of antiderivative is not currently implemented for this CAS.
[In]
int(-tan(b*x-c)*tan(b*x+a),x)
[Out]
x+I/b/(exp(2*I*(a+c))-1)*ln(1+exp(2*I*(b*x+a)))*exp(2*I*(a+c))+I/b/(exp(2*I*(a+c))-1)*ln(1+exp(2*I*(b*x+a)))-I
/b/(exp(2*I*(a+c))-1)*ln(exp(2*I*(a+c))+exp(2*I*(b*x+a)))*exp(2*I*(a+c))-I/b/(exp(2*I*(a+c))-1)*ln(exp(2*I*(a+
c))+exp(2*I*(b*x+a)))
________________________________________________________________________________________
maxima [B] time = 0.35, size = 290, normalized size = 8.53 \[ \frac {{\left (b \cos \left (2 \, a + 2 \, c\right )^{2} + b \sin \left (2 \, a + 2 \, c\right )^{2} - 2 \, b \cos \left (2 \, a + 2 \, c\right ) + b\right )} x - {\left (\cos \left (2 \, a + 2 \, c\right )^{2} + \sin \left (2 \, a + 2 \, c\right )^{2} - 1\right )} \arctan \left (\sin \left (2 \, b x\right ) - \sin \left (2 \, a\right ), \cos \left (2 \, b x\right ) + \cos \left (2 \, a\right )\right ) + {\left (\cos \left (2 \, a + 2 \, c\right )^{2} + \sin \left (2 \, a + 2 \, c\right )^{2} - 1\right )} \arctan \left (\sin \left (2 \, b x\right ) + \sin \left (2 \, c\right ), \cos \left (2 \, b x\right ) + \cos \left (2 \, c\right )\right ) + \log \left (\cos \left (2 \, b x\right )^{2} + 2 \, \cos \left (2 \, b x\right ) \cos \left (2 \, a\right ) + \cos \left (2 \, a\right )^{2} + \sin \left (2 \, b x\right )^{2} - 2 \, \sin \left (2 \, b x\right ) \sin \left (2 \, a\right ) + \sin \left (2 \, a\right )^{2}\right ) \sin \left (2 \, a + 2 \, c\right ) - \log \left (\cos \left (2 \, b x\right )^{2} + 2 \, \cos \left (2 \, b x\right ) \cos \left (2 \, c\right ) + \cos \left (2 \, c\right )^{2} + \sin \left (2 \, b x\right )^{2} + 2 \, \sin \left (2 \, b x\right ) \sin \left (2 \, c\right ) + \sin \left (2 \, c\right )^{2}\right ) \sin \left (2 \, a + 2 \, c\right )}{b \cos \left (2 \, a + 2 \, c\right )^{2} + b \sin \left (2 \, a + 2 \, c\right )^{2} - 2 \, b \cos \left (2 \, a + 2 \, c\right ) + b} \]
Verification of antiderivative is not currently implemented for this CAS.
[In]
integrate(-tan(b*x-c)*tan(b*x+a),x, algorithm="maxima")
[Out]
((b*cos(2*a + 2*c)^2 + b*sin(2*a + 2*c)^2 - 2*b*cos(2*a + 2*c) + b)*x - (cos(2*a + 2*c)^2 + sin(2*a + 2*c)^2 -
1)*arctan2(sin(2*b*x) - sin(2*a), cos(2*b*x) + cos(2*a)) + (cos(2*a + 2*c)^2 + sin(2*a + 2*c)^2 - 1)*arctan2(
sin(2*b*x) + sin(2*c), cos(2*b*x) + cos(2*c)) + log(cos(2*b*x)^2 + 2*cos(2*b*x)*cos(2*a) + cos(2*a)^2 + sin(2*
b*x)^2 - 2*sin(2*b*x)*sin(2*a) + sin(2*a)^2)*sin(2*a + 2*c) - log(cos(2*b*x)^2 + 2*cos(2*b*x)*cos(2*c) + cos(2
*c)^2 + sin(2*b*x)^2 + 2*sin(2*b*x)*sin(2*c) + sin(2*c)^2)*sin(2*a + 2*c))/(b*cos(2*a + 2*c)^2 + b*sin(2*a + 2
*c)^2 - 2*b*cos(2*a + 2*c) + b)
________________________________________________________________________________________
mupad [B] time = 5.00, size = 196, normalized size = 5.76 \[ \frac {\frac {x}{2}+x\,\left ({\sin \left (a+c\right )}^2-\frac {1}{2}\right )}{{\sin \left (a+c\right )}^2}+\frac {\frac {\sin \left (2\,a+2\,c\right )\,\ln \left ({\sin \left (2\,a+2\,c\right )}^2\,2{}\mathrm {i}-{\sin \left (a+b\,x\right )}^2\,2{}\mathrm {i}+{\sin \left (3\,a+2\,c+b\,x\right )}^2\,2{}\mathrm {i}+\sin \left (4\,a+4\,c\right )+\sin \left (6\,a+4\,c+2\,b\,x\right )-\sin \left (2\,a+2\,b\,x\right )\right )}{2}-\frac {\sin \left (2\,a+2\,c\right )\,\ln \left ({\sin \left (2\,a+c+b\,x\right )}^2\,2{}\mathrm {i}+{\sin \left (2\,a+2\,c\right )}^2\,2{}\mathrm {i}-{\sin \left (c-b\,x\right )}^2\,2{}\mathrm {i}+\sin \left (4\,a+4\,c\right )+\sin \left (4\,a+2\,c+2\,b\,x\right )+\sin \left (2\,c-2\,b\,x\right )\right )}{2}}{b\,{\sin \left (a+c\right )}^2} \]
Verification of antiderivative is not currently implemented for this CAS.
[In]
int(tan(a + b*x)*tan(c - b*x),x)
[Out]
(x/2 + x*(sin(a + c)^2 - 1/2))/sin(a + c)^2 + ((sin(2*a + 2*c)*log(sin(4*a + 4*c) + sin(6*a + 4*c + 2*b*x) - s
in(2*a + 2*b*x) + sin(2*a + 2*c)^2*2i - sin(a + b*x)^2*2i + sin(3*a + 2*c + b*x)^2*2i))/2 - (sin(2*a + 2*c)*lo
g(sin(4*a + 4*c) + sin(4*a + 2*c + 2*b*x) + sin(2*c - 2*b*x) + sin(2*a + c + b*x)^2*2i + sin(2*a + 2*c)^2*2i -
sin(c - b*x)^2*2i))/2)/(b*sin(a + c)^2)
________________________________________________________________________________________
sympy [B] time = 8.52, size = 7720, normalized size = 227.06 \[ \text {result too large to display} \]
Verification of antiderivative is not currently implemented for this CAS.
[In]
integrate(-tan(b*x-c)*tan(b*x+a),x)
[Out]
Piecewise((0, Eq(a, 0) & Eq(b, 0) & Eq(c, 0)), (2*b*x*tan(a)/(2*b*tan(a)**2 + 2*b) + 2*log(tan(b*x) - 1/tan(a)
)/(2*b*tan(a)**2 + 2*b) - log(tan(b*x)**2 + 1)/(2*b*tan(a)**2 + 2*b), Eq(c, 0)), (-4*b*x*tan(c)**2*tan(b*x)/(2
*b*tan(c)**5*tan(b*x) + 2*b*tan(c)**4 + 4*b*tan(c)**3*tan(b*x) + 4*b*tan(c)**2 + 2*b*tan(c)*tan(b*x) + 2*b) -
4*b*x*tan(c)/(2*b*tan(c)**5*tan(b*x) + 2*b*tan(c)**4 + 4*b*tan(c)**3*tan(b*x) + 4*b*tan(c)**2 + 2*b*tan(c)*tan
(b*x) + 2*b) - 2*log(tan(b*x) + 1/tan(c))*tan(c)**3*tan(b*x)/(2*b*tan(c)**5*tan(b*x) + 2*b*tan(c)**4 + 4*b*tan
(c)**3*tan(b*x) + 4*b*tan(c)**2 + 2*b*tan(c)*tan(b*x) + 2*b) - 2*log(tan(b*x) + 1/tan(c))*tan(c)**2/(2*b*tan(c
)**5*tan(b*x) + 2*b*tan(c)**4 + 4*b*tan(c)**3*tan(b*x) + 4*b*tan(c)**2 + 2*b*tan(c)*tan(b*x) + 2*b) + 2*log(ta
n(b*x) + 1/tan(c))*tan(c)*tan(b*x)/(2*b*tan(c)**5*tan(b*x) + 2*b*tan(c)**4 + 4*b*tan(c)**3*tan(b*x) + 4*b*tan(
c)**2 + 2*b*tan(c)*tan(b*x) + 2*b) + 2*log(tan(b*x) + 1/tan(c))/(2*b*tan(c)**5*tan(b*x) + 2*b*tan(c)**4 + 4*b*
tan(c)**3*tan(b*x) + 4*b*tan(c)**2 + 2*b*tan(c)*tan(b*x) + 2*b) + log(tan(b*x)**2 + 1)*tan(c)**3*tan(b*x)/(2*b
*tan(c)**5*tan(b*x) + 2*b*tan(c)**4 + 4*b*tan(c)**3*tan(b*x) + 4*b*tan(c)**2 + 2*b*tan(c)*tan(b*x) + 2*b) + lo
g(tan(b*x)**2 + 1)*tan(c)**2/(2*b*tan(c)**5*tan(b*x) + 2*b*tan(c)**4 + 4*b*tan(c)**3*tan(b*x) + 4*b*tan(c)**2
+ 2*b*tan(c)*tan(b*x) + 2*b) - log(tan(b*x)**2 + 1)*tan(c)*tan(b*x)/(2*b*tan(c)**5*tan(b*x) + 2*b*tan(c)**4 +
4*b*tan(c)**3*tan(b*x) + 4*b*tan(c)**2 + 2*b*tan(c)*tan(b*x) + 2*b) - log(tan(b*x)**2 + 1)/(2*b*tan(c)**5*tan(
b*x) + 2*b*tan(c)**4 + 4*b*tan(c)**3*tan(b*x) + 4*b*tan(c)**2 + 2*b*tan(c)*tan(b*x) + 2*b) - 2*tan(c)**2/(2*b*
tan(c)**5*tan(b*x) + 2*b*tan(c)**4 + 4*b*tan(c)**3*tan(b*x) + 4*b*tan(c)**2 + 2*b*tan(c)*tan(b*x) + 2*b) - 2/(
2*b*tan(c)**5*tan(b*x) + 2*b*tan(c)**4 + 4*b*tan(c)**3*tan(b*x) + 4*b*tan(c)**2 + 2*b*tan(c)*tan(b*x) + 2*b),
Eq(a, -atan(tan(c)) - pi*floor((c - pi/2)/pi) - pi*floor(c/pi - 1/2))), (0, Eq(b, 0)), (-2*b*x*tan(c)/(2*b*tan
(c)**2 + 2*b) + 2*log(tan(b*x) + 1/tan(c))/(2*b*tan(c)**2 + 2*b) - log(tan(b*x)**2 + 1)/(2*b*tan(c)**2 + 2*b),
Eq(a, 0)), (2*b*x*tan(a)**2/(2*b*tan(a)**3*tan(c)**2 + 2*b*tan(a)**3 + 2*b*tan(a)**2*tan(c)**3 + 2*b*tan(a)**
2*tan(c) + 2*b*tan(a)*tan(c)**2 + 2*b*tan(a) + 2*b*tan(c)**3 + 2*b*tan(c)) - 2*b*x*tan(c)**2/(2*b*tan(a)**3*ta
n(c)**2 + 2*b*tan(a)**3 + 2*b*tan(a)**2*tan(c)**3 + 2*b*tan(a)**2*tan(c) + 2*b*tan(a)*tan(c)**2 + 2*b*tan(a) +
2*b*tan(c)**3 + 2*b*tan(c)) + 2*log(tan(b*x) - 1/tan(a))*tan(a)*tan(c)**2/(2*b*tan(a)**3*tan(c)**2 + 2*b*tan(
a)**3 + 2*b*tan(a)**2*tan(c)**3 + 2*b*tan(a)**2*tan(c) + 2*b*tan(a)*tan(c)**2 + 2*b*tan(a) + 2*b*tan(c)**3 + 2
*b*tan(c)) + 2*log(tan(b*x) - 1/tan(a))*tan(a)/(2*b*tan(a)**3*tan(c)**2 + 2*b*tan(a)**3 + 2*b*tan(a)**2*tan(c)
**3 + 2*b*tan(a)**2*tan(c) + 2*b*tan(a)*tan(c)**2 + 2*b*tan(a) + 2*b*tan(c)**3 + 2*b*tan(c)) + 2*log(tan(b*x)
+ 1/tan(c))*tan(a)**2*tan(c)/(2*b*tan(a)**3*tan(c)**2 + 2*b*tan(a)**3 + 2*b*tan(a)**2*tan(c)**3 + 2*b*tan(a)**
2*tan(c) + 2*b*tan(a)*tan(c)**2 + 2*b*tan(a) + 2*b*tan(c)**3 + 2*b*tan(c)) + 2*log(tan(b*x) + 1/tan(c))*tan(c)
/(2*b*tan(a)**3*tan(c)**2 + 2*b*tan(a)**3 + 2*b*tan(a)**2*tan(c)**3 + 2*b*tan(a)**2*tan(c) + 2*b*tan(a)*tan(c)
**2 + 2*b*tan(a) + 2*b*tan(c)**3 + 2*b*tan(c)) - log(tan(b*x)**2 + 1)*tan(a)**2*tan(c)/(2*b*tan(a)**3*tan(c)**
2 + 2*b*tan(a)**3 + 2*b*tan(a)**2*tan(c)**3 + 2*b*tan(a)**2*tan(c) + 2*b*tan(a)*tan(c)**2 + 2*b*tan(a) + 2*b*t
an(c)**3 + 2*b*tan(c)) - log(tan(b*x)**2 + 1)*tan(a)*tan(c)**2/(2*b*tan(a)**3*tan(c)**2 + 2*b*tan(a)**3 + 2*b*
tan(a)**2*tan(c)**3 + 2*b*tan(a)**2*tan(c) + 2*b*tan(a)*tan(c)**2 + 2*b*tan(a) + 2*b*tan(c)**3 + 2*b*tan(c)) -
log(tan(b*x)**2 + 1)*tan(a)/(2*b*tan(a)**3*tan(c)**2 + 2*b*tan(a)**3 + 2*b*tan(a)**2*tan(c)**3 + 2*b*tan(a)**
2*tan(c) + 2*b*tan(a)*tan(c)**2 + 2*b*tan(a) + 2*b*tan(c)**3 + 2*b*tan(c)) - log(tan(b*x)**2 + 1)*tan(c)/(2*b*
tan(a)**3*tan(c)**2 + 2*b*tan(a)**3 + 2*b*tan(a)**2*tan(c)**3 + 2*b*tan(a)**2*tan(c) + 2*b*tan(a)*tan(c)**2 +
2*b*tan(a) + 2*b*tan(c)**3 + 2*b*tan(c)), True))*tan(a) - Piecewise((0, Eq(a, 0) & Eq(b, 0) & Eq(c, 0)), (2*b*
x*tan(a)/(2*b*tan(a)**2 + 2*b) + 2*log(tan(b*x) - 1/tan(a))/(2*b*tan(a)**2 + 2*b) - log(tan(b*x)**2 + 1)/(2*b*
tan(a)**2 + 2*b), Eq(c, 0)), (-4*b*x*tan(c)**2*tan(b*x)/(2*b*tan(c)**5*tan(b*x) + 2*b*tan(c)**4 + 4*b*tan(c)**
3*tan(b*x) + 4*b*tan(c)**2 + 2*b*tan(c)*tan(b*x) + 2*b) - 4*b*x*tan(c)/(2*b*tan(c)**5*tan(b*x) + 2*b*tan(c)**4
+ 4*b*tan(c)**3*tan(b*x) + 4*b*tan(c)**2 + 2*b*tan(c)*tan(b*x) + 2*b) - 2*log(tan(b*x) + 1/tan(c))*tan(c)**3*
tan(b*x)/(2*b*tan(c)**5*tan(b*x) + 2*b*tan(c)**4 + 4*b*tan(c)**3*tan(b*x) + 4*b*tan(c)**2 + 2*b*tan(c)*tan(b*x
) + 2*b) - 2*log(tan(b*x) + 1/tan(c))*tan(c)**2/(2*b*tan(c)**5*tan(b*x) + 2*b*tan(c)**4 + 4*b*tan(c)**3*tan(b*
x) + 4*b*tan(c)**2 + 2*b*tan(c)*tan(b*x) + 2*b) + 2*log(tan(b*x) + 1/tan(c))*tan(c)*tan(b*x)/(2*b*tan(c)**5*ta
n(b*x) + 2*b*tan(c)**4 + 4*b*tan(c)**3*tan(b*x) + 4*b*tan(c)**2 + 2*b*tan(c)*tan(b*x) + 2*b) + 2*log(tan(b*x)
+ 1/tan(c))/(2*b*tan(c)**5*tan(b*x) + 2*b*tan(c)**4 + 4*b*tan(c)**3*tan(b*x) + 4*b*tan(c)**2 + 2*b*tan(c)*tan(
b*x) + 2*b) + log(tan(b*x)**2 + 1)*tan(c)**3*tan(b*x)/(2*b*tan(c)**5*tan(b*x) + 2*b*tan(c)**4 + 4*b*tan(c)**3*
tan(b*x) + 4*b*tan(c)**2 + 2*b*tan(c)*tan(b*x) + 2*b) + log(tan(b*x)**2 + 1)*tan(c)**2/(2*b*tan(c)**5*tan(b*x)
+ 2*b*tan(c)**4 + 4*b*tan(c)**3*tan(b*x) + 4*b*tan(c)**2 + 2*b*tan(c)*tan(b*x) + 2*b) - log(tan(b*x)**2 + 1)*
tan(c)*tan(b*x)/(2*b*tan(c)**5*tan(b*x) + 2*b*tan(c)**4 + 4*b*tan(c)**3*tan(b*x) + 4*b*tan(c)**2 + 2*b*tan(c)*
tan(b*x) + 2*b) - log(tan(b*x)**2 + 1)/(2*b*tan(c)**5*tan(b*x) + 2*b*tan(c)**4 + 4*b*tan(c)**3*tan(b*x) + 4*b*
tan(c)**2 + 2*b*tan(c)*tan(b*x) + 2*b) - 2*tan(c)**2/(2*b*tan(c)**5*tan(b*x) + 2*b*tan(c)**4 + 4*b*tan(c)**3*t
an(b*x) + 4*b*tan(c)**2 + 2*b*tan(c)*tan(b*x) + 2*b) - 2/(2*b*tan(c)**5*tan(b*x) + 2*b*tan(c)**4 + 4*b*tan(c)*
*3*tan(b*x) + 4*b*tan(c)**2 + 2*b*tan(c)*tan(b*x) + 2*b), Eq(a, -atan(tan(c)) - pi*floor((c - pi/2)/pi) - pi*f
loor(c/pi - 1/2))), (0, Eq(b, 0)), (-2*b*x*tan(c)/(2*b*tan(c)**2 + 2*b) + 2*log(tan(b*x) + 1/tan(c))/(2*b*tan(
c)**2 + 2*b) - log(tan(b*x)**2 + 1)/(2*b*tan(c)**2 + 2*b), Eq(a, 0)), (2*b*x*tan(a)**2/(2*b*tan(a)**3*tan(c)**
2 + 2*b*tan(a)**3 + 2*b*tan(a)**2*tan(c)**3 + 2*b*tan(a)**2*tan(c) + 2*b*tan(a)*tan(c)**2 + 2*b*tan(a) + 2*b*t
an(c)**3 + 2*b*tan(c)) - 2*b*x*tan(c)**2/(2*b*tan(a)**3*tan(c)**2 + 2*b*tan(a)**3 + 2*b*tan(a)**2*tan(c)**3 +
2*b*tan(a)**2*tan(c) + 2*b*tan(a)*tan(c)**2 + 2*b*tan(a) + 2*b*tan(c)**3 + 2*b*tan(c)) + 2*log(tan(b*x) - 1/ta
n(a))*tan(a)*tan(c)**2/(2*b*tan(a)**3*tan(c)**2 + 2*b*tan(a)**3 + 2*b*tan(a)**2*tan(c)**3 + 2*b*tan(a)**2*tan(
c) + 2*b*tan(a)*tan(c)**2 + 2*b*tan(a) + 2*b*tan(c)**3 + 2*b*tan(c)) + 2*log(tan(b*x) - 1/tan(a))*tan(a)/(2*b*
tan(a)**3*tan(c)**2 + 2*b*tan(a)**3 + 2*b*tan(a)**2*tan(c)**3 + 2*b*tan(a)**2*tan(c) + 2*b*tan(a)*tan(c)**2 +
2*b*tan(a) + 2*b*tan(c)**3 + 2*b*tan(c)) + 2*log(tan(b*x) + 1/tan(c))*tan(a)**2*tan(c)/(2*b*tan(a)**3*tan(c)**
2 + 2*b*tan(a)**3 + 2*b*tan(a)**2*tan(c)**3 + 2*b*tan(a)**2*tan(c) + 2*b*tan(a)*tan(c)**2 + 2*b*tan(a) + 2*b*t
an(c)**3 + 2*b*tan(c)) + 2*log(tan(b*x) + 1/tan(c))*tan(c)/(2*b*tan(a)**3*tan(c)**2 + 2*b*tan(a)**3 + 2*b*tan(
a)**2*tan(c)**3 + 2*b*tan(a)**2*tan(c) + 2*b*tan(a)*tan(c)**2 + 2*b*tan(a) + 2*b*tan(c)**3 + 2*b*tan(c)) - log
(tan(b*x)**2 + 1)*tan(a)**2*tan(c)/(2*b*tan(a)**3*tan(c)**2 + 2*b*tan(a)**3 + 2*b*tan(a)**2*tan(c)**3 + 2*b*ta
n(a)**2*tan(c) + 2*b*tan(a)*tan(c)**2 + 2*b*tan(a) + 2*b*tan(c)**3 + 2*b*tan(c)) - log(tan(b*x)**2 + 1)*tan(a)
*tan(c)**2/(2*b*tan(a)**3*tan(c)**2 + 2*b*tan(a)**3 + 2*b*tan(a)**2*tan(c)**3 + 2*b*tan(a)**2*tan(c) + 2*b*tan
(a)*tan(c)**2 + 2*b*tan(a) + 2*b*tan(c)**3 + 2*b*tan(c)) - log(tan(b*x)**2 + 1)*tan(a)/(2*b*tan(a)**3*tan(c)**
2 + 2*b*tan(a)**3 + 2*b*tan(a)**2*tan(c)**3 + 2*b*tan(a)**2*tan(c) + 2*b*tan(a)*tan(c)**2 + 2*b*tan(a) + 2*b*t
an(c)**3 + 2*b*tan(c)) - log(tan(b*x)**2 + 1)*tan(c)/(2*b*tan(a)**3*tan(c)**2 + 2*b*tan(a)**3 + 2*b*tan(a)**2*
tan(c)**3 + 2*b*tan(a)**2*tan(c) + 2*b*tan(a)*tan(c)**2 + 2*b*tan(a) + 2*b*tan(c)**3 + 2*b*tan(c)), True))*tan
(c) + Piecewise((0, Eq(a, 0) & Eq(b, 0) & Eq(c, 0)), (2*b*x*tan(a)/(2*b*tan(a)**3 + 2*b*tan(a)) + 2*log(tan(b*
x) - 1/tan(a))/(2*b*tan(a)**3 + 2*b*tan(a)) + log(tan(b*x)**2 + 1)*tan(a)**2/(2*b*tan(a)**3 + 2*b*tan(a)), Eq(
c, 0)), (-b*x*tan(c)**4*tan(b*x)/(b*tan(c)**6*tan(b*x) + b*tan(c)**5 + 2*b*tan(c)**4*tan(b*x) + 2*b*tan(c)**3
+ b*tan(c)**2*tan(b*x) + b*tan(c)) - b*x*tan(c)**3/(b*tan(c)**6*tan(b*x) + b*tan(c)**5 + 2*b*tan(c)**4*tan(b*x
) + 2*b*tan(c)**3 + b*tan(c)**2*tan(b*x) + b*tan(c)) + b*x*tan(c)**2*tan(b*x)/(b*tan(c)**6*tan(b*x) + b*tan(c)
**5 + 2*b*tan(c)**4*tan(b*x) + 2*b*tan(c)**3 + b*tan(c)**2*tan(b*x) + b*tan(c)) + b*x*tan(c)/(b*tan(c)**6*tan(
b*x) + b*tan(c)**5 + 2*b*tan(c)**4*tan(b*x) + 2*b*tan(c)**3 + b*tan(c)**2*tan(b*x) + b*tan(c)) + 2*log(tan(b*x
) + 1/tan(c))*tan(c)**3*tan(b*x)/(b*tan(c)**6*tan(b*x) + b*tan(c)**5 + 2*b*tan(c)**4*tan(b*x) + 2*b*tan(c)**3
+ b*tan(c)**2*tan(b*x) + b*tan(c)) + 2*log(tan(b*x) + 1/tan(c))*tan(c)**2/(b*tan(c)**6*tan(b*x) + b*tan(c)**5
+ 2*b*tan(c)**4*tan(b*x) + 2*b*tan(c)**3 + b*tan(c)**2*tan(b*x) + b*tan(c)) - log(tan(b*x)**2 + 1)*tan(c)**3*t
an(b*x)/(b*tan(c)**6*tan(b*x) + b*tan(c)**5 + 2*b*tan(c)**4*tan(b*x) + 2*b*tan(c)**3 + b*tan(c)**2*tan(b*x) +
b*tan(c)) - log(tan(b*x)**2 + 1)*tan(c)**2/(b*tan(c)**6*tan(b*x) + b*tan(c)**5 + 2*b*tan(c)**4*tan(b*x) + 2*b*
tan(c)**3 + b*tan(c)**2*tan(b*x) + b*tan(c)) + tan(c)**2/(b*tan(c)**6*tan(b*x) + b*tan(c)**5 + 2*b*tan(c)**4*t
an(b*x) + 2*b*tan(c)**3 + b*tan(c)**2*tan(b*x) + b*tan(c)) + 1/(b*tan(c)**6*tan(b*x) + b*tan(c)**5 + 2*b*tan(c
)**4*tan(b*x) + 2*b*tan(c)**3 + b*tan(c)**2*tan(b*x) + b*tan(c)), Eq(a, -atan(tan(c)) - pi*floor((c - pi/2)/pi
) - pi*floor(c/pi - 1/2))), (0, Eq(b, 0)), (2*b*x*tan(c)/(2*b*tan(c)**3 + 2*b*tan(c)) - 2*log(tan(b*x) + 1/tan
(c))/(2*b*tan(c)**3 + 2*b*tan(c)) - log(tan(b*x)**2 + 1)*tan(c)**2/(2*b*tan(c)**3 + 2*b*tan(c)), Eq(a, 0)), (2
*b*x*tan(a)**2*tan(c)/(2*b*tan(a)**3*tan(c)**2 + 2*b*tan(a)**3 + 2*b*tan(a)**2*tan(c)**3 + 2*b*tan(a)**2*tan(c
) + 2*b*tan(a)*tan(c)**2 + 2*b*tan(a) + 2*b*tan(c)**3 + 2*b*tan(c)) + 2*b*x*tan(a)*tan(c)**2/(2*b*tan(a)**3*ta
n(c)**2 + 2*b*tan(a)**3 + 2*b*tan(a)**2*tan(c)**3 + 2*b*tan(a)**2*tan(c) + 2*b*tan(a)*tan(c)**2 + 2*b*tan(a) +
2*b*tan(c)**3 + 2*b*tan(c)) + 2*b*x*tan(a)/(2*b*tan(a)**3*tan(c)**2 + 2*b*tan(a)**3 + 2*b*tan(a)**2*tan(c)**3
+ 2*b*tan(a)**2*tan(c) + 2*b*tan(a)*tan(c)**2 + 2*b*tan(a) + 2*b*tan(c)**3 + 2*b*tan(c)) + 2*b*x*tan(c)/(2*b*
tan(a)**3*tan(c)**2 + 2*b*tan(a)**3 + 2*b*tan(a)**2*tan(c)**3 + 2*b*tan(a)**2*tan(c) + 2*b*tan(a)*tan(c)**2 +
2*b*tan(a) + 2*b*tan(c)**3 + 2*b*tan(c)) + 2*log(tan(b*x) - 1/tan(a))*tan(c)**2/(2*b*tan(a)**3*tan(c)**2 + 2*b
*tan(a)**3 + 2*b*tan(a)**2*tan(c)**3 + 2*b*tan(a)**2*tan(c) + 2*b*tan(a)*tan(c)**2 + 2*b*tan(a) + 2*b*tan(c)**
3 + 2*b*tan(c)) + 2*log(tan(b*x) - 1/tan(a))/(2*b*tan(a)**3*tan(c)**2 + 2*b*tan(a)**3 + 2*b*tan(a)**2*tan(c)**
3 + 2*b*tan(a)**2*tan(c) + 2*b*tan(a)*tan(c)**2 + 2*b*tan(a) + 2*b*tan(c)**3 + 2*b*tan(c)) - 2*log(tan(b*x) +
1/tan(c))*tan(a)**2/(2*b*tan(a)**3*tan(c)**2 + 2*b*tan(a)**3 + 2*b*tan(a)**2*tan(c)**3 + 2*b*tan(a)**2*tan(c)
+ 2*b*tan(a)*tan(c)**2 + 2*b*tan(a) + 2*b*tan(c)**3 + 2*b*tan(c)) - 2*log(tan(b*x) + 1/tan(c))/(2*b*tan(a)**3*
tan(c)**2 + 2*b*tan(a)**3 + 2*b*tan(a)**2*tan(c)**3 + 2*b*tan(a)**2*tan(c) + 2*b*tan(a)*tan(c)**2 + 2*b*tan(a)
+ 2*b*tan(c)**3 + 2*b*tan(c)) + log(tan(b*x)**2 + 1)*tan(a)**2/(2*b*tan(a)**3*tan(c)**2 + 2*b*tan(a)**3 + 2*b
*tan(a)**2*tan(c)**3 + 2*b*tan(a)**2*tan(c) + 2*b*tan(a)*tan(c)**2 + 2*b*tan(a) + 2*b*tan(c)**3 + 2*b*tan(c))
- log(tan(b*x)**2 + 1)*tan(c)**2/(2*b*tan(a)**3*tan(c)**2 + 2*b*tan(a)**3 + 2*b*tan(a)**2*tan(c)**3 + 2*b*tan(
a)**2*tan(c) + 2*b*tan(a)*tan(c)**2 + 2*b*tan(a) + 2*b*tan(c)**3 + 2*b*tan(c)), True)) - Piecewise((-x, Eq(a,
0) & Eq(b, 0) & Eq(c, 0)), (-2*b*x/(2*b*tan(a)**2 + 2*b) + 2*log(tan(b*x) - 1/tan(a))*tan(a)/(2*b*tan(a)**2 +
2*b) - log(tan(b*x)**2 + 1)*tan(a)/(2*b*tan(a)**2 + 2*b), Eq(c, 0)), (b*x*tan(c)**3*tan(b*x)/(b*tan(c)**5*tan(
b*x) + b*tan(c)**4 + 2*b*tan(c)**3*tan(b*x) + 2*b*tan(c)**2 + b*tan(c)*tan(b*x) + b) + b*x*tan(c)**2/(b*tan(c)
**5*tan(b*x) + b*tan(c)**4 + 2*b*tan(c)**3*tan(b*x) + 2*b*tan(c)**2 + b*tan(c)*tan(b*x) + b) - b*x*tan(c)*tan(
b*x)/(b*tan(c)**5*tan(b*x) + b*tan(c)**4 + 2*b*tan(c)**3*tan(b*x) + 2*b*tan(c)**2 + b*tan(c)*tan(b*x) + b) - b
*x/(b*tan(c)**5*tan(b*x) + b*tan(c)**4 + 2*b*tan(c)**3*tan(b*x) + 2*b*tan(c)**2 + b*tan(c)*tan(b*x) + b) - 2*l
og(tan(b*x) + 1/tan(c))*tan(c)**2*tan(b*x)/(b*tan(c)**5*tan(b*x) + b*tan(c)**4 + 2*b*tan(c)**3*tan(b*x) + 2*b*
tan(c)**2 + b*tan(c)*tan(b*x) + b) - 2*log(tan(b*x) + 1/tan(c))*tan(c)/(b*tan(c)**5*tan(b*x) + b*tan(c)**4 + 2
*b*tan(c)**3*tan(b*x) + 2*b*tan(c)**2 + b*tan(c)*tan(b*x) + b) + log(tan(b*x)**2 + 1)*tan(c)**2*tan(b*x)/(b*ta
n(c)**5*tan(b*x) + b*tan(c)**4 + 2*b*tan(c)**3*tan(b*x) + 2*b*tan(c)**2 + b*tan(c)*tan(b*x) + b) + log(tan(b*x
)**2 + 1)*tan(c)/(b*tan(c)**5*tan(b*x) + b*tan(c)**4 + 2*b*tan(c)**3*tan(b*x) + 2*b*tan(c)**2 + b*tan(c)*tan(b
*x) + b) + tan(c)**3/(b*tan(c)**5*tan(b*x) + b*tan(c)**4 + 2*b*tan(c)**3*tan(b*x) + 2*b*tan(c)**2 + b*tan(c)*t
an(b*x) + b) + tan(c)/(b*tan(c)**5*tan(b*x) + b*tan(c)**4 + 2*b*tan(c)**3*tan(b*x) + 2*b*tan(c)**2 + b*tan(c)*
tan(b*x) + b), Eq(a, -atan(tan(c)) - pi*floor((c - pi/2)/pi) - pi*floor(c/pi - 1/2))), (-x, Eq(b, 0)), (-2*b*x
/(2*b*tan(c)**2 + 2*b) - 2*log(tan(b*x) + 1/tan(c))*tan(c)/(2*b*tan(c)**2 + 2*b) + log(tan(b*x)**2 + 1)*tan(c)
/(2*b*tan(c)**2 + 2*b), Eq(a, 0)), (-2*b*x*tan(a)**2*tan(c)/(2*b*tan(a)**3*tan(c)**2 + 2*b*tan(a)**3 + 2*b*tan
(a)**2*tan(c)**3 + 2*b*tan(a)**2*tan(c) + 2*b*tan(a)*tan(c)**2 + 2*b*tan(a) + 2*b*tan(c)**3 + 2*b*tan(c)) - 2*
b*x*tan(a)*tan(c)**2/(2*b*tan(a)**3*tan(c)**2 + 2*b*tan(a)**3 + 2*b*tan(a)**2*tan(c)**3 + 2*b*tan(a)**2*tan(c)
+ 2*b*tan(a)*tan(c)**2 + 2*b*tan(a) + 2*b*tan(c)**3 + 2*b*tan(c)) - 2*b*x*tan(a)/(2*b*tan(a)**3*tan(c)**2 + 2
*b*tan(a)**3 + 2*b*tan(a)**2*tan(c)**3 + 2*b*tan(a)**2*tan(c) + 2*b*tan(a)*tan(c)**2 + 2*b*tan(a) + 2*b*tan(c)
**3 + 2*b*tan(c)) - 2*b*x*tan(c)/(2*b*tan(a)**3*tan(c)**2 + 2*b*tan(a)**3 + 2*b*tan(a)**2*tan(c)**3 + 2*b*tan(
a)**2*tan(c) + 2*b*tan(a)*tan(c)**2 + 2*b*tan(a) + 2*b*tan(c)**3 + 2*b*tan(c)) + 2*log(tan(b*x) - 1/tan(a))*ta
n(a)**2*tan(c)**2/(2*b*tan(a)**3*tan(c)**2 + 2*b*tan(a)**3 + 2*b*tan(a)**2*tan(c)**3 + 2*b*tan(a)**2*tan(c) +
2*b*tan(a)*tan(c)**2 + 2*b*tan(a) + 2*b*tan(c)**3 + 2*b*tan(c)) + 2*log(tan(b*x) - 1/tan(a))*tan(a)**2/(2*b*ta
n(a)**3*tan(c)**2 + 2*b*tan(a)**3 + 2*b*tan(a)**2*tan(c)**3 + 2*b*tan(a)**2*tan(c) + 2*b*tan(a)*tan(c)**2 + 2*
b*tan(a) + 2*b*tan(c)**3 + 2*b*tan(c)) - 2*log(tan(b*x) + 1/tan(c))*tan(a)**2*tan(c)**2/(2*b*tan(a)**3*tan(c)*
*2 + 2*b*tan(a)**3 + 2*b*tan(a)**2*tan(c)**3 + 2*b*tan(a)**2*tan(c) + 2*b*tan(a)*tan(c)**2 + 2*b*tan(a) + 2*b*
tan(c)**3 + 2*b*tan(c)) - 2*log(tan(b*x) + 1/tan(c))*tan(c)**2/(2*b*tan(a)**3*tan(c)**2 + 2*b*tan(a)**3 + 2*b*
tan(a)**2*tan(c)**3 + 2*b*tan(a)**2*tan(c) + 2*b*tan(a)*tan(c)**2 + 2*b*tan(a) + 2*b*tan(c)**3 + 2*b*tan(c)) -
log(tan(b*x)**2 + 1)*tan(a)**2/(2*b*tan(a)**3*tan(c)**2 + 2*b*tan(a)**3 + 2*b*tan(a)**2*tan(c)**3 + 2*b*tan(a
)**2*tan(c) + 2*b*tan(a)*tan(c)**2 + 2*b*tan(a) + 2*b*tan(c)**3 + 2*b*tan(c)) + log(tan(b*x)**2 + 1)*tan(c)**2
/(2*b*tan(a)**3*tan(c)**2 + 2*b*tan(a)**3 + 2*b*tan(a)**2*tan(c)**3 + 2*b*tan(a)**2*tan(c) + 2*b*tan(a)*tan(c)
**2 + 2*b*tan(a) + 2*b*tan(c)**3 + 2*b*tan(c)), True))*tan(a)*tan(c)
________________________________________________________________________________________