БДЗ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)