import numpy as np

def f(x):
    return 1/x


start = 1.0
stop = 3.0
n = 10
delta_x = (stop - start)/n

x = np.linspace(start, stop, n+1)


integral = 0.5*(f(start) + f(stop))*delta_x

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

eksakt_verdi = np.log(3)    
print('integral = ', integral)
print('integral - eksakt_verdi = ',integral-eksakt_verdi)