I think I know this one but the solution given is different than what I get so maybe I've missed something important.
From $f(x,y)=x^2y^2+x^2y+2y^2-4y$ we have $f'_x=2xy^2+2xy$ and $f'_y=2x^2y+x^2+4y-4$. Setting both equations equal to $0$ we have $(1)$ $2xy^2+2xy=2xy(y+1)=0$ and $(2)$ $2x^2y+x^2+4y-4=0$.
The points that satisfy $(1)$ is $(0,y)$, $(x,0)$ and $(x,-1)$ where $x,y\in\mathbb{R}$.
Putting $(0,y)$ in $(2)$ gives $4y-4=0$ and it follows that $(0,1)$ is one stationary point.
Putting $(x,0)$ in $(2)$ gives $x^2-4=0\iff x=2\text{ or } x=-2$ so the stationary points are $(2,0)$ and $(-2,0)$.
Similarly with $(x,-1)$ we get $-x^2-8=0$ which does not have any real solutions.
So the stationary points I got are $(0,1)$, $(2,0)$ and $(-2,0)$ but the solution that was given to me is only $(0,1)$. I know the other two points satisfy $(1)$ and $(2)$ but I'm trying to make sure I haven't missed something.