$$\begin{align}&f_x=e^{xy}\left(1+xy-y^2\right)=0\iff y^2-xy=1\\{}\\&f_y=e^{xy}\left(-1+x^2-xy\right)=0\iff x^2-xy=1\end{align}$$
Resumiendo las dos ecuaciones anteriores obtenemos $\;(x-y)(x+y)=0\implies x=\pm y\;$ y tenemos muchos puntos críticos.
Ahora:
$$\begin{align}&f_{xx}=e^{xy}\left(2y+xy^2-y^3\right)\\{}\\&f_{xy}=e^{xy}\left(2x+x^2y-xy^2\right)\\{}\\&f_{yy}=e^{xy}\left(-2x+x^3-x^2y\right)\end{align}$$
Evalúa ahora el hessiano en los puntos críticos y verás que hay muchos puntos máx.-mín. en la plaza. Por último, tendrás que compararlo con los valores de la función en el perímetro del cuadrado.
Descargo de responsabilidad: Compruebe cuidadosamente los cálculos anteriores, que sólo vienen a servirle de guía. Puede haber errores.