= 1. Let f(x, y, z) = xyz + x +y +z + 1. Find the gradient vf and divergence div(vf), and then calculate curl(vf) at point (1, 1, 1).