4. (10%) Find the local maxima of the function: $f(x, y) = \left(\frac{1}{3}x^3 + xy^2 - 4xy + 1\right)$