import numpy as np

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


n = 10000  # number of x-values
x1 = -2
x2 = 2
x = np.linspace(x1, x2, n)

x_sol = -np.inf

for i in range(0, n-1):    
    if f(x[i])*f(x[i+1]) <= 0:
        x_sol = (x[i] + x[i+1])/2
        print('Losning:', x_sol)