import scipy
from scipy import integrate

value, error = integrate.quad(lambda x: x**2, 0, 1)

print(f"SciPy {scipy.__version__}")
print(f"integral={value:.6f}")
print(f"estimated_error={error:.2e}")