The quadratic term of a Taylor expansion of a multivariate scalar valued function can be expressed in terms of the Hessian. Is there a similar form for vector valued functions, in the sense that all partial derivatives can be arranged in a matrix or tensor?