def f(x):
    return 8.0 - x**2

x0 = 1.0
x1 = 2.0
for aantal_rh in range (100,100000,100):
    oppervlakte = 0
    h = (x1-x0)/aantal_rh
    for i in range(aantal_rh + 1):
        oppervlakte += f(x0 + (i * h))
    oppervlakte -= (f(x0) + f(x1)) / 2
    print(h * oppervlakte)
