Compute the length and the direction of the gradient of the function u = 1/r where r =sqrt(x^2+y^2+z^2) at a point M(x0, y0, z0). (Remark: the direction should be described by the unit vector having the same direction as ∇f)