এসসিপি nquad

from scipy import integrate

result, error = integrate.nquad(
    lambda x, y: x ** 2 + y ** 2,
    [[0, 1], [0, 2]]
)

print(result)
print(error)