import numpy as np

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


start = -2.0
stop = 3.0
n = 2000
delta_x = (stop - start)/n

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

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

print('Integralet med foroverintegrasjon =', integral)
print('Eksakt verdi av integralet =', 35/3)