Метод наискорейшего спуска

clc
clear all
%---- Задаём функцию
y=@(x)211*x(1).^2-420*x(1).*x(2)+211*x(2).^2-192*x(1)+50*x(2)-25;
%------
dyx1=@(x)422*x(1) - 420*x(2) - 192;
dyx2=@(x)422*x(2) - 420*x(1) + 50;
%------
E=0.0001;
% Метод наискорейшего спуска
a=1;
x=[0,0];
k=0;
while (sqrt(dyx1(x)^2+dyx2(x)^2)>E)
    f=@(a)y(x-a*[dyx1(x),dyx2(x)]);
    local_min=dihotomy(f,0,3,E);
    a=local_min(1);
    x=x-a*[dyx1(x),dyx2(x)];
    k=k+1;
end
 
%График
[x1,x2] = meshgrid(-200:1:200);
z=211*x1.^2+-420*x1.*x2+211*x2.^2-192*x1+50*x2-25;
surf(x1,x2,z)
shading interp;
text(x(1), x(2), 'Минимум', 'color', 'r'); 
%Ответ
k
x
y(x)

function [ ret ] = dihotomy(cb,a,b,eps)
    d2=eps/10;
    n=0;
    x_min=a;
    while (b-a)>eps
        x_min = (b+a)/2;
        fp=cb(x_min+d2);
        fm=cb(x_min-d2);
        if fp<fm
            a=x_min-d2;
        else
            b=x_min+d2;
        end
        n=n+1;
    end
    ret=[x_min,n];
end

a=local_min(1);

while (sqrt(dyx1(x)^2+dyx2(x)^2+dyx3(x))>E)
    f=@(a)y(x-a*[dyx1(x),dyx2(x),dyx3(x)]);
    local_min=dihotomy(f,0,3,E);
    a=local_min(1);
    x=x-a*[dyx1(x),dyx2(x),dyx3(x)];
    k=k+1;
end