The divergence of a vector, as defined by the usual formula in vector calculus, is not an intrinsic operator which can be defined in a coordinate independent manner. Still, for a number of purposes it is helpful to have this simple procedure.

The command Divergence is part of the DifferentialGeometry:-Tools package, and so can be used in the form Divergence(...) only after executing the commands with(DifferentialGeometry) and with(Tools) in that order.  It can always be used in the long form DifferentialGeometry:-Tools:-Divergence.

$\mathrm{with}\left(\mathrm{DifferentialGeometry}\right)\:$$\mathrm{with}\left(\mathrm{Tools}\right)\:$

$\mathrm{DGsetup}\left(\left[x\,y\,z\right]\,M\right)\:$

$X≔\mathrm{DGzip}\left(\left[a\,b\,c\right]\left(x\,y\,z\right)\,\left[\mathrm{D\_x}\,\mathrm{D\_y}\,\mathrm{D\_z}\right]\,''plus''\right)$

$X≔a\left(x\,y\,z\right)\mathrm{D\_x}+b\left(x\,y\,z\right)\mathrm{D\_y}+c\left(x\,y\,z\right)\mathrm{D\_z}$

$\mathrm{Divergence}\left(X\right)$

$\frac{\partial }{\partial z}c\left(x\,y\,z\right)+\frac{\partial }{\partial y}b\left(x\,y\,z\right)+\frac{\partial }{\partial x}a\left(x\,y\,z\right)$