Metoda dihotomica

function dihotomic_method(f, interval, epsilon)
    delta = 0.02
    a = first(interval)
    b = last(interval)
    while (abs(b-a) > epsilon)
        lambda = (a+b)/2 - delta
        miu = (a+b)/2 + delta
        if (f(lambda) < f(miu))
            b = miu
        else
            a = lambda
        end
    end
    x_optim = (b+a)/2
end

f(x) = x*(x-1.5)
interval = range(0, stop=1)
epsilon = 0.1
x_optim = dihotomic_method(f, interval, epsilon)