如何计算数值积分中的误差（摆线）

Question

import math
def trap(g, a, b, n):

    h=(b-a)/float(n);
    estimativa = (f(a)+f(b))/2;
    for i in range(n):
        x=a+i*h
        estimativa+=g(x)
    return estimativa*h

def f(x):
    return math.exp(-x**2)

print("approximate value of integral: ",trap(f,0,10,81660))
print("exact value: ",0.8862269254527)

我正在使用飞人方法来计算该积分，有人可以帮助我如何计算一个条件下的近似值吗？错误必须小于10 ^ -7