import sympy as sp from IPython.display import display a,x = sp.symbols('a x',real=True) eq=[x**a,sp.sin(x),sp.cos(x), 1/sp.cos(x)**2,1/sp.sin(x)**2,sp.E**x,a**x] for f in eq: display(sp.Eq(sp.Integral(f, x), sp.simplify(sp.Integral(f, x).doit())))