import numpy as np

def f(x):
    return np.sqrt(1 + (4*x)**2) 


start = -0.5
stop = 2
n = 100
delta_x = (stop - start)/n

x = np.linspace(start + delta_x/2, stop - delta_x/2, n)


kurvelengde = 0.0
for i in range(0, n, 1):
    kurvelengde = kurvelengde + f(x[i])*delta_x


print('kurvelengde =', f'{kurvelengde:.3f}')