БДЗ1 / dichotomy

clc,clear
syms x
fun = x.^5 - 2*sin(x);
a = 0;
b = 1;
e = 0.001;
delta = e/2;
x_min = a;
f_min = subs(fun,x_min);
A = [];
B = [];

while abs(b-a)/2 > e
x1 = (b+a-delta)/2;
x2 = (b+a+delta)/2;

if subs(fun,x1) > subs(fun,x2)
a = x1;
else
b = x2;
end

A(end+1) = a;
B(end+1) = b;

end

A((end-2):end)
B((end-2):end)

x_min = (a+b)/2
f_min = subs(fun,x_min)