def f(t):
    return(1+t**2)

def integrale(a, b, n):
    S = 0
    for i in range (1, n+1): 
        S = S + f(a+i*((b-a)/n))
        I =((b-a)/n)*S
    return (I)
print (integrale (0, 1, 10000))