import sympy as sym x, n = sym.symbols('x, n') print(sym.integrate(sym.log(x),x)) print(sym.integrate(x**n, x))