БДЗ1 / golden_ratio

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

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

if subs(fun,x1) > subs(fun,x2)
a = x1;
x2 = a+b-x2;
else
b = x2;
x2 = x1;
x1 = a+b-x1;
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)