Минимизация унимодальной функции методом дихотомии

clear all, clc
S = 2;
R = -3;
% начальные границы диапазона
a0 = 1;
b0 = 2;
% S = input('Введите S: ');
% R = input('Введите R: ');
% a0 = input('Введите а: ')
% b0 = input('Введите b: ')
% функция
f = @(x) x - x.^2-S/2 - 1 - R;
 
fi = (1+sqrt(5))/2; % константа, пропорция золотого сечения
eps = 0.001; % точность
a = a0;
b = b0;
% начальные условия
x1 = b - (b-a)/fi;
x2 = a + (b-a)/fi;
y1 = f(x1);
y2 = f(x2);
 
while (b-a)>eps
    if abs(y1) > abs(y2)
        a = x1;
        x1 = b - (b-a)/fi;
        x2 = a + (b-a)/fi;
        y1 = f(x1);
        y2 = f(x2);
    else
        b = x2;
        x1 = b - (b-a)/fi;
        x2 = a + (b-a)/fi;
        y1 = f(x1);
        y2 = f(x2);
    end
end
x = (a+b)/2 % корень
 
% рисуем график
t = a0:(b0-a0)/299:b0;
y = f(t);
plot(t,y,'b',x,0,'or')

clear all; clc;
% Основные условия
f=@(x) cos(2*x)+x-5; %это искомая функция
a=1;  % находим точку минимума нижняя граница оси У
b=-4; % находим точку минимума верхняя граница оси Х
e=0.0001; % точность до 10-4
 
% Осноной цикл 
while abs(b-a)>e
    c=(b+a)/2;
    if sign(f(c))==sign(f(a))
        a=c;
    else
        b=c;
    end
end
disp(['min y = ' num2str(f(c),6)]);
disp(['Ответ x=' num2str(c,3)]);
 
fplot(f)
grid on
hold on
plot(c, f(c),'or')
 
% disp('для проверки точный ответ =');
% % disp(solve('cos(2*x)+x-5'));
% 
% syms x;
% disp(solve(f(x)));