import sympy as sym x, a = sym.symbols('x, a') print(sym.integrate(x**2, (x, -1, 1))) print(sym.integrate(sym.exp(-a*x**2),(x,-sym.oo, sym.oo)))