The gradient at point \( P(3, -1) \) is \( \nabla f(3, -1) = \left( 0.4615, -0.3077 \right) \), the unit vector is \( \mathbf{u} \approx \left( 0.9487, 0.3162 \right) \), and the directional derivative is approximately \( 0.3406 \).
Thus, the final answers are:
\[
\boxed{\nabla f(3, -1) = \left( 0.4615, -0.3077 \right)}
\]
\[
\boxed{\mathbf{u} \approx \left( 0.9487, 0.3162 \right)}
\]
\[
\boxed{D_{\mathbf{u}} f \approx 0.3406}
\]