This is simple in a rectangular basis. It can be found in arbitrary basis by playing with partials. This is pretty messy.
Another method for orthogonal coordinate systems:
Let and . Then the coordinates are orthogonal if satisfy .
The length element in the new coordinate system is . So .
For spherical coordinates,
. Define .
For volume integrals, where .
For the gradient: The inverse matrix is needed so, so gradient =
For divergence: