Поиск максимума в массиве методом касательных и методом Ньютона-Рафсона

%clc;
%clear;
 
x  = 1:length(Z1_var4_max(:,1));
 
u  = zeros(1,100);
 
%f = zeros(1,length(x));
% ОПРЕДЕЛЕНИЕ МИНИМИЗИРУЕМОЙ ФУНКЦИИ 
%alfha = 0.01;
%f = alfha*(x/100-75).^2 + 50;
f = -(Z1_var4_max(:,2)) + 2 ;
 
u(1) = x(1);                                % начальная точка  
deriv = f(1+1)-f(1);
g = deriv*x + f(1);                         % построена касательная
 
u(2)  = (min(g)-f(1))/deriv;                % следующая точка
deriv = f(u(2))-f(u(2)-1);
g1    = deriv*x + f(u(2)) - deriv*u(2);     % построена касательная
h1    = max(g,g1);                          % построение ломаной
minimum = min(h1);                          % минимум ломаной
ind  = find( h1 == minimum );               % поиск индекса
u(3)  = ind(1);                             % следующая точка
 
deriv_L = f(u(3))-f(u(3)-1);                % левая
deriv_R = f(u(3)+1)-f(u(3));                % провая
deriv = min(deriv_L,deriv_R);               % выбрали минимальную производную
g2    = deriv*x + f(u(3)) - deriv*u(3);     % построена касательная
h2    = max(h1,g2);                         % построение ломаной
minimum = min(h2);                          % минимум ломаной
ind  = find( h2 == minimum );               % поиск индекса    
u(4)  = ind(1);
 
deriv_L = f(u(4))-f(u(4)-1);                % левая
deriv_R = f(u(4)+1)-f(u(4));                % провая
% if( (deriv_R >= 0) && (deriv_L <= 0) )      % следующая точка
%     disp ( '_К_О_Н_Е_Ц_'  );
%     disp ( u(4) );
%     flag = 0;
% else
%     flag = 1;
% end
 
u_ = u(4);
h_ = h2;
flag = 1;
i = 0;
n = 0;
while(i <= 100)                               % далее цикл
    disp ( i );
    disp ( u_ );
%   
    deriv = min(deriv_L,deriv_R);
    g_    = deriv*x + f(u_) - deriv*u_;       % построена касательная
    h__   = max(h_,g_);                       % построение ломаной
    minimum = min(h__);
    ind  = find( h__ == minimum );
    u__  = ind(1);
 
     deriv_L = f(u__)-f(u__-1);                % левая
     deriv_R = f(u__+1)-f(u__);                % провая
     
     if( u_ == u__ )
         n = n+1;
     end
     
     u_ = u__;
     h_ = h__;
     
%      deriv = min(deriv_L,deriv_R);
%      g__   = deriv*x + f(u__) - deriv*u__;      % построена касательная
%      h___   = max(h__,g_);                       % построение ломаной
%      minimum = min(h___);
%      ind  = find( h___ == minimum );
%      u___  = ind(1);
% 
%      deriv_L = f(u___)-f(u___-1);                % левая
%      deriv_R = f(u___+1)-f(u___);                % провая
          
    i = i+1;
    
    if( ( (deriv_L <=0) && (deriv_R >= 0)) || n >= 10 )
        i = 101;
        disp ( '_THE_END_'  );
        disp ( u__ );
        
    end
        
end
 
    %plot( [1:length(Z1_var4_max(:,1))], f, [1:length(Z1_var4_max(:,1))], g, [1:length(Z1_var4_max(:,1))], g1, [1:length(Z1_var4_max(:,1))], g2, [1:length(Z1_var4_max(:,1))], h__ );
    %plot( x, f, x, g_, x, h__ ); 
    %plot( [1:length(Z1_var4_max(:,1))], g1 );
    %axis ([ 0 10000 -60 110 ]);
    %grid on;
    
    xx = [1:length(Z1_var4_max(:,1))];
    plot( xx, f(xx), xx, g1(xx), xx, g2(xx), xx, g_(xx) );
    
    %disp(f(u__)  );
    
    % ВЫВОД РЕЗУЛЬТАТОВ
    disp( ' x значение' );
    disp( Z1_var4_max(u__,1)  ); 
    disp( ' y значение' );
    disp( Z1_var4_max(u__,2)  );