In[27]:= << VectorAnalysis`
In[28]:= Laplacian[f[r, th, phi], Spherical[r, th, phi]]
$\frac{\csc (\text{th}) \left(r^2 \sin (\text{th}) f^{(2,0,0)}(r,\text{th},\phi )+2 r \sin (\text{th})
f^{(1,0,0)}(r,\text{th},\phi )+\sin (\text{th}) f^{(0,2,0)}(r,\text{th},\phi )+\cos (\text{th})
f^{(0,1,0)}(r,\text{th},\phi )+\csc (\text{th}) f^{(0,0,2)}(r,\text{th},\phi )\right)}{r^2}$
Regards
Hans Dolhaine