Вы частные производные правильно нашли

далее Сокращаем на неравную нулю экспоненту и совершенно по-школьному решаем систему. Получаем точку

и окружность

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

.