Consider function f(x,y,z) = z^2e^{xy}. Find the gradient vector of f at the point (-1,1,3). Compute the directional derivative of f at the point (-1,1,3) in the direction of v = 3i + j - 5k. Give a unit vector that points in the direction of maximum increase for the function f(x,y,z) at P(-1,1,3). Find the maximum rate of change of f at P(-1,1,3).