Find the directional derivative of $\displaystyle f(x,y,z) = 2xy+z^{2}$ at the point $(1,-3,3)$ in the direction of the maximum rate of change of $f$.

$f_{\boldsymbol{u}} \, (1,-3,3) = D_{\boldsymbol{u}} \, f(1,-3,3) =$