ClearAll[x] integrand = Cos[x]^2/Sqrt[Cos[x]^4 + Cos[x]^2 + 1]; res = Integrate[integrand, x]; TeXForm[res]
<< Rubi` ClearAll[x] integrand = Cos[x]^2/Sqrt[Cos[x]^4 + Cos[x]^2 + 1]; res = Int[integrand, x]; TeXForm[res]
restart; integrand:= cos(x)^2/sqrt(cos(x)^4 + cos(x)^2 + 1); res:=int(integrand,x); latex(res)
set output tex off setSimplifyDenomsFlag(true) integrand:= cos(x)^2/sqrt(cos(x)^4 + cos(x)^2 + 1); res:=integrate(integrand,x); latex(res)
integrand : cos(x)^2/sqrt(cos(x)^4 + cos(x)^2 + 1); res : integrate(integrand,x); latex(res)
integrand := cos(x)^2/sqrt(cos(x)^4 + cos(x)^2 + 1); res := integrate(integrand,x); latex(res)
>python Python 3.7.3 (default, Mar 27 2019, 22:11:17) [GCC 7.3.0] :: Anaconda, Inc. on linux from sympy import * x = symbols('x') integrand = integrand = cos(x)**2/sqrt(cos(x)**4 + cos(x)**2 + 1); res = integrate(integrand,x); latex(res)
evalin(symengine,'int(cos(x)^2/sqrt(cos(x)^4 + cos(x)^2 + 1),x)')