For Easom function (or any other functions): $f(x,y)=\cos x*\cos y*\exp\left(-(〖x-π)〗^2-(y〖-π)〗^2 \right)$, defined in a range $0\le x,y\le 10$, the plot will look as follows: enter image description here
Is it possible to analytically compute maximum functional value within a circle of radius '$R$' (using double integral), which would be a subset of the original region.