from scipy.integrate import quad
result = quad(my_function, a, b)
"""
quad 1-Dim function
dblquad, tplquad, nquad 2, 3, N-Dim respectively
my_function is a Python-def function
a special function from scipy (check documentation)
or a lambda function
a, b: integration interval
"""