Какие методы поиска минимума функции нескольких переменных можно использовать

Есть язык R и функция нескольких переменных:

Holt <- function(data_table, coef)
{
  alpha <- coef[1]
  beta <- coef[2]
  ldt = length(data_table)
  ut = bt = numeric(ldt)
  for(i in (1:(ldt)))
  {
    if(i == 1)
    {
      ut[i] = data_table[i]
      bt[i] = data_table[i+1] - data_table[i]
    }
    else
    {
      ut[i] = alpha * data_table[i-1] + (1 - alpha) * (ut[i-1] + bt[i-1])
      bt[i] = beta * (ut[i] - ut[i-1]) + (1 - beta) * bt[i-1]
    }
  }
  aet = abs(data_table[2:ldt] - ut[2:ldt])
  mape = mean(aet[1:(ldt-1)] / data_table[2:ldt]) * 100
  return(mape)
}

Нужно найти минимальное значение mape при alpha, beta в диапазоне от 0.0001 до 0.9999. Какой алгоритм лучше всего использовать для более точного результата? Если можно с реализацией. Использовал метод покоординатного спуска - оказался неточным.