I should analyze whether the function
$$f(x,y) = x^2y + x y^2 \text{ where } x,y > 0$$
is (quasi-) concave or convex.
Thus, as usual, I set up the Hessian as
$$ D^2f(x,y) = \left( \begin{array}{cc} 2y & 2x + 2y \\ 2x + 2y & 2x \\ \end{array} \right) $$
which, given the constraints $x,y > 0$, is indefinite.
So I need the bordered Hessian:
$$ H= \left( \begin{array}{ccc} 0 & 2xy + y^2 & x^2 + 2xy \\ 2xy + y^2 & 2y & 2x + 2y \\ x^2 + 2xy & 2x+2y & 2x \end{array} \right) $$
and check its determinant to find out about its properties. I end up with
$$ \det|H| = 2x(2x^3 + x^3 y - 2x^2 y + 4x^2 y^2 + 4 x y^3 + 2 y^4)$$
which, considering $x,y > 0$, looks quite positive to me, so $f$ would be quasiconcave. But how can I prove that $\det|H| > 0$? Or is there any easier approach?
Thanks!