Para cualquier suficientemente función derivable
$F:(x, y, z): \Bbb R^3 \to \Bbb R \tag{1}$
que puede escribirse como una función de $r = \sqrt{x^2 + y^2 + z^2}$ solamente, tenemos la fórmula general
$\nabla F(r) = \dfrac{dF(r)}{dr} \nabla r \tag{2}$
que se desprende directamente de la regla de la cadena, a saber.
$\dfrac{\partial F}{\partial x} = \dfrac{dF(r)}{dr} \dfrac{\partial r}{\partial x}, \tag{3}$
con el correspondiente epxpressions para el $y$ $z$ derivados. También,
$\dfrac{\partial r}{\partial x} = \dfrac{\partial \sqrt{x^2 + y^2 + z^2}}{\partial x} = \dfrac{1}{2 \sqrt{x^2 + y^2 + z^2}}(2x) = \dfrac{x}{\sqrt{x^2 + y^2 + z^2}} =\dfrac{x}{r}, \tag{4}$
y del mismo modo
$\dfrac{\partial r}{\partial y} =\dfrac{y}{r}, \tag{5}$
$\dfrac{\partial r}{\partial z} =\dfrac{z}{r}. \tag{6}$
Así tenemos
$\nabla r = (\dfrac{\partial r}{\partial x}, \dfrac{\partial r}{\partial y}, \dfrac{\partial r}{\partial z}) = \dfrac{1}{r}(x, y, z). \tag{7}$
En la realización de los cálculos anteriores que, por supuesto, suponemos que permanecer lejos del punto de $(0, 0, 0)$ donde $r$ no es diferenciable.
A continuación, se deduce de (2) y (7) que
$\nabla F(r) = \dfrac{dF(r)}{dr} \dfrac{1}{r}(x, y, z). \tag{8}$
Tomando
$F(r) = f(r) = \dfrac{1}{\cosh^3 (kr)} = \cosh^{-3}(kr) \tag{9}$
tenemos
$f'(r) = \dfrac{df(r)}{dr} = -3k \cosh^{-4}(kr) \sinh(kr)$
$ = -3k \cosh^{-3}(kr) \dfrac{\sinh (kr)}{\cosh (kr)} = -3k \cosh^{-3}(kr) \tanh (kr) = -3k \tanh (kr) \mathrm{sech^3}(kr); \tag{10}$
ahora sigue de (8)-(10) que
$\nabla f(r) = -\dfrac{3k \tanh (kr) \mathrm{sech^3}(kr)}{r}(x, y, z) = -3k \tanh (kr) \mathrm{sech^3}(kr) \dfrac{(x, y, z)}{r}; \tag{11}$
vemos que el $x$-componente de $\nabla f$ contiene un factor de $x$ etc.
Quizás vale la pena señalar que la expresión
$\dfrac{(x, y, z)}{r} = \dfrac{(x, y, z)}{\sqrt{x^2 + y^2 + z^2}} \tag{12}$
ocurre en (11) es, en realidad,$\vec{e}_r$, el vector unitario en la dirección radial:
$\vec{e}_r = \dfrac{(x, y, z)}{\sqrt{x^2 + y^2 + z^2}}; \tag{13}$
por lo tanto (8) puede escribirse
$\nabla F(r) = \dfrac{dF(r)}{dr} \vec{e}_r, \tag{14}$
un útil fórmula general.
En el cierre, se nota que no hay nada en especial en tres dimensiones sobre estas fórmulas. De hecho, con $F:\Bbb R^n \to \Bbb R$ y coordina $x_1, x_2, \dots, x_n$ tenemos $r = \sqrt{\sum_{i = 1}^n x_i^2}$ y
$\vec{e}_r = \dfrac{(x_1, x_2, \ldots, x_n)}{r}; \tag{15}$
entonces (14) aún se une.