Локальные максимумы функции

[serge25]

Как подсчитать локальные максимумы для следующей системы?

%Файл–функция

function [dydt] = hodhux(t, y) 
% вводим переменные 
v = y(1); 
n = y(2); 
m = y(3);
h = y(4);
 
a_n = 0.01*(10-v)/(exp((10-v)/10)-1);
b_n = 0.125*exp(-v/80);
a_m = 0.1*(25-v)/(exp((25-v)/10)-1);
b_m = 4*exp(-v/18);
a_h = 0.07*exp(-v/20);
b_h = 1/(exp((30-v)/10)+1);
E_k = -12;
E_na = 120;
E_l = 10.6;
g_k = 36;
g_na = 120;
g_l = 0.3;
c = 1;
I = 6.29107;
n_massiv;
 
dydt = [(0.001*q+ I -g_k*n^4*(v-E_k)-g_na*m^3*h*(v-E_na)-g_l*(v-E_l))/c; 
a_n*(1-n)-b_n*n;
a_m*(1-m)-b_m*m;
a_h*(1-h)-b_h*h];
end

%Программа-сценарий n_massiv, задающая белый шум
for i=0:100
q=rand(1);
i=i+1;
A(1,i)=q;
end

%консольная команда
n_massiv%формируем массив значений белого шума
y0 = [-0.1, 0.2, 0.1, 0.1];
% интервалы интегрирования
tspan = [0, 100];
% решаем систему ДУ
[t_out,y_out] = ode45(@hodhux,tspan,y0);
plot(H,K);set(get(gcf,'CurrentAxes'),'FontName','Times New Roman Cyr','FontSize',14);
grid on; title('График зависимости потенциала от времени');
xlabel(время'); ylabel(значение потенциала');

Добавлено спустя 1 час 52 минуты 30 секунд:

Спасибо за внимание ответ найден)