Поиск максимума функции многих переменных

На примере типовой функции peaks. Разбиваем область определения мелкой сеткой, получаем массив значений функции и в нем ищем максимум. Короче говоря, перебор.

Z = peaks(100);
surf(Z)
Zmax = max(max(Z));
[i,j] = find(Z==Zmax)
hold on
plot3(j,i,Z(i,j),'m.','MarkerSize',24)

Результаты:

i =

    77

j =

    50

Фиолетовая точка на графике показывает расположение найденного максимума.