БДЗ2 / midpoint

clc,clear
syms x
fun = x.^5 - 2*sin(x);
a = 0;
b = 1;
e = 0.001;
der = diff(fun);
xmin = (a+b)/2;
A = [];
B = [];

while abs(subs(der,(a+b)/2))>=e

der_value = subs(der,(a+b)/2);
xmin = (a+b)/2;

if der_value >0
b = xmin;
else
a = xmin;
end
A(end+1) = xmin;
B(end+1) = subs(fun,xmin);

end

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

xmin
ymin = subs(fun,xmin)