Using Lagrange multiplier , find minimum value of $$f(x,y)=xy(x^2 + y^2) +4$$ , Given that $g(x,y)=x^2 + y^2 +xy -1 = 0 , $for all values of $ x,y \in \mathbb R$.
My attempt
So i formed a function $$F(x,y , \lambda ) = f(x,y) + {\lambda}g(x,y) $$ , where $\lambda$ is the Lagrange multiplier .
I then solved the partial derivatives w.r.t $x,y$ and $\lambda$.
$${{\partial {F} } \over {\partial x}} = 3x^2 y + y^3 + 2x{ \lambda} + y{\lambda} = 0 \space ....(1)$$ $${{\partial {F} } \over {\partial y}} = 3y^2 x + x^3 + 2y{ \lambda} + x{\lambda} = 0 \space ....(2)$$
$${{\partial {F} } \over {\partial \lambda}} = x^2 + y^2 +xy -1 = 0 \space ....(3).$$
Solving equations (1) ,(2) and (3) , I got $$xy = \frac12$$
Therefore , putting value of $xy$ in $x^2 + y^2 + xy -1 = 0$ I got $$x^2 + y^2 = \frac12$$.
And hence ,putting these value in $f(x,y)$ I got the value $$f_{min} (x,y) =\frac{17}{4}$$
However , this is not the correct answer , the correct answer is 2 . Please help me out and tell where am I going wrong in the method . I know a trigonometric solution to this question , but I want a solution using Lagrange multipliers . Any help is appreciated .

The equations obtained are correct, maybe you did some mistake deriving the solutions.
Notably, multiplying the first equation by $y$ and the second one by $x$ we obtain
and subtracting we obtain
$$y^4-x^4+y^2{\lambda}-x^2{\lambda}=0 \iff (y^2-x^2)(x^2+y^2+\lambda)=0$$
$$\iff y^2-x^2=0 \quad \lor \quad x^2+y^2+\lambda=0$$
which leads to solutions for $x^2=y^2 \iff x=\pm y$, which can easily be obtained from the constraint, indeed the case $\lambda=-(x^2+y^2)$ leads to no solution.