電磁気学の分野の問題になってしまうのですが、

r=(x^2+y^2+z^2)^(1/2)

E(ベクトル)=(r(ベクトル)/r^3)(1+kr)e^(-kr) kは定数

のとき

∇・E(ベクトル)を求めてください

模範解答にどうしてもならず、ここから先の問いに答えられません
普通の偏微分でやってみたのですが、何回やってもできません
よろしくお願いします