ЛР1 / lab1_sechenie

function [f_min, x_min] = lab1_sechenie(a,b,eps)
    tau = (sqrt(5) - 1)/2;
    a_start = a;
    b_start = b;
    
    x1 = a+(1-tau)*(b-a);
    x2 = a + tau*(b-a);
    f1 = myfun(x1);
    f2 = myfun(x2);
    counter = 2;
    eps_n = (b-a)/2;

    while eps_n>eps
        if f1 >= f2
            a = x1;
            x1 = x2;
            f1 = f2;
            x2 = a+b-x2;
            f2 = myfun(x2);
        else
            b = x2;
            x2 = x1;
            f2 = f1;
            x1 = a+b-x1;
            f1 = myfun(x1);
        end
        eps_n = tau*eps_n;
        counter = counter + 1;
    end

    x_min = (a+b)/2;
    f_min = myfun(x_min);
    str=sprintf('Минимум функции на отрезке [%d; %d] - %3f, в точке %3f',a_start, b_start, f_min, x_min);
    disp(str);
    str = sprintf('Количество вычислений функции - %d', counter);
    disp(str);

end