ЛР1 / lab1_parabol

function [f_min, x_min] = lab1_parabol(a, b, eps)
a_start = a;
b_start = b;

x1 = input('Введите x1: ');
x2 = input('Введите x2: ');
x3 = input('Введите x3: ');
while (x1 >= x2) || (x3 <= x2) || (myfun(x1) <= myfun(x2)) || (myfun(x3) <= myfun(x2))
x1 = input('Введите x1: ');
x2 = input('Введите x2: ');
x3 = input('Введите x3: ');
end

fx1 = myfun(x1);
fx3 = myfun(x3);
fx2 = myfun(x2);
counter = 3;

a1 = (fx2-fx1)/(x2-x1); % Самая первая итерация по вычислению x_min, x1, x2, x3
a2 = 1/(x3-x2) * ( (fx3-fx1)/(x3-x1) - (fx2-fx1)/(x2-x1) );
x_min = 0.5*(x1+x2-a1/a2);
fx_min=myfun(x_min);
if (x_min<x2)
if (fx_min<fx2)
x3=x2; fx3=fx2; x2=x_min; fx2=fx_min;
else
x1=x_min; fx1=fx_min;
end

else
if (fx_min<fx2)
x1=x2; fx1=fx2; x2=x_min; fx2=fx_min;
else
x3=x_min; fx3=fx_min;
end
end
counter = counter + 1;

delta = 1000; % Т.к. первый раз дельту вычислить невозможно

while abs(delta) > eps
x_min_old = x_min;
a1 = (fx2-fx1)/(x2-x1);
a2 = 1/(x3-x2) * ( (fx3-fx1)/(x3-x1) - (fx2-fx1)/(x2-x1) );
x_min = 0.5*(x1+x2-a1/a2);
fx_min=myfun(x_min);

if (x_min<x2)
if (fx_min<fx2)
x3=x2; fx3=fx2; x2=x_min; fx2=fx_min;
else
x1=x_min; fx1=fx_min;
end

else
if (fx_min<fx2)
x1=x2; fx1=fx2; x2=x_min; fx2=fx_min;
else
x3=x_min; fx3=fx_min;
end
end
counter = counter + 1;
x_min_new = x_min;
delta = abs(x_min_old - x_min_new);

end

f_min = myfun(x_min);
str=sprintf('Минимум функции на отрезке [%d; %d] - %3f, в точке %3f',a_start, b_start, f_min, x_min);
disp(str);
str = sprintf('Количество вычислений функции - %d', counter);
disp(str);
end