No way. In spherical coordinates the field is simply [itex]\hat r/r^2[/itex]. There is no phi or theta dependence in there, since it's radial remember? All you have to do is lookup the divergence in spherical coordinates and do a simply derivative.