\[ \int \cot (e+f x) \left (a+b \sec ^3(e+f x)\right ) \, dx \]
Optimal antiderivative \[ \frac {\left (a +b \right ) \ln \left (1-\cos \left (f x +e \right )\right )}{2 f}+\frac {\left (a -b \right ) \ln \left (1+\cos \left (f x +e \right )\right )}{2 f}+\frac {b \sec \left (f x +e \right )}{f} \]
command
integrate(cot(f*x+e)*(a+b*sec(f*x+e)^3),x, algorithm="giac")
Giac 1.9.0-11 via sagemath 9.6 output
\[ \frac {{\left (a + b\right )} \log \left (\frac {{\left | -\cos \left (f x + e\right ) + 1 \right |}}{{\left | \cos \left (f x + e\right ) + 1 \right |}}\right ) - 2 \, a \log \left ({\left | -\frac {\cos \left (f x + e\right ) - 1}{\cos \left (f x + e\right ) + 1} + 1 \right |}\right ) + \frac {4 \, b}{\frac {\cos \left (f x + e\right ) - 1}{\cos \left (f x + e\right ) + 1} + 1}}{2 \, f} \]
Giac 1.7.0 via sagemath 9.3 output
\[ \text {Exception raised: NotImplementedError} \]________________________________________________________________________________________