The divergence of a vector field \( \mathbf{P} = P_x \hat{i} + P_y \hat{j} + P_z \hat{k} \) is given by:
\[
\nabla \cdot \mathbf{P} = \frac{\partial P_x}{\partial x} + \frac{\partial P_y}{\partial y} + \frac{\partial P_z}{\partial z}
\]
Given \( P_x = x^2 y \) and \( P_y = xy \):
\[
\frac{\partial (x^2 y)}{\partial x} = 2xy, \quad \frac{\partial (xy)}{\partial y} = x
\]
\[
\nabla \cdot \mathbf{P} = 2xy + x
\]