How can I calculate using Lagrange multipliers the minimum and maximum of the function $$F(x,y) = \frac{ 2(x^2 +y^2)}{x^2 + y^2 + 1} - {16x^2 + 1}$$
in the domain $D= \Big\{x,y \in \mathbb{R} \ ; \ { x^2 +y^2} \leq 1 \Big\}$.
I am having trouble calculating the partial derivatives of this function and would appreciate any help. Thanks.

I can suggest the following solution: 1) calculate all local min/max at $\mathbb{R}^2$ using standard calculus tools - discard all points outside domain $D$, and then 2) calculate possible min/max using Lagrange multiplier with the constrain $x^2+y^2=1$
Finally, compare the results - the min/max may lie on the boundary line of the domain D.
The partial derivatives are:
$$\frac{\partial f}{\partial x}(x,y)=\frac{2(x^2+y^2)\cdot 2x-4x(x^2+y^2+1)}{(x^2+y^2+1)^2}-32x=\frac{4 x (1-8 (x^2+y^2+1)^2)}{(x^2+y^2+1)^2}$$ $$\frac{\partial f}{\partial y}(x,y)=\frac{4y(x^2+y^2+1)-2(x^2+y^2)\cdot 2y}{(x^2+y^2+1)^2}=\frac{4 y}{\left(x^2+y^2+1\right)^2}$$