В данной работе рассматривается задача глобальной оптимизации гладкой функции нескольких переменных, заданной на кубоиде. Поиск решения осуществляется  с помощью  вспомогательной функции, полученной путем особого преобразования целевой функции. Вспомогательная функция является функцией одной переменной, ноль которой совпадает со значением глобального минимума целевой функции. Поэтому для решения поставленной задачи применялся метод деления отрезка пополам. Результаты настоящей работы были выявлены на основе большого числа вычислительных экспериментов, проведенных на тестовых функциях с помощью предложенного метода. Эти результаты сформулированы в виде трех теорем и теоретически доказаны.  В первой теореме определены условия, которые указывают на промежуток в котором находится значение глобального минимума. Вторая теорема выражает сходимость итерационной последовательности к значению глобального минимума. В третьей теореме установлена линейная скорость сходимости итерационной процедуры. В качестве примера рассмотрена многоэкстремальная функция Экли двух переменных, определенная в квадрате с центром в начале координат.