程序如下：
Clear["Global`*"];
V = q/(4 Pi*\[Epsilon]0)*(1/r1 - 1/r2)
cond = {r1 -> R - l/2*Cos[\[Theta]], r2 -> R + l/2*Cos[\[Theta]]}
V = V /. cond // Simplify
Series[V, {l, 0, 6}] // Normal
V1 = pCos[\[Theta]]/(4 \[Pi] (R^2) \[Epsilon]0 )
Solve[V1 == c, R]
到这儿都对，结果是
{{R -> -(Sqrt[pCos[\[Theta]]]/(
2 Sqrt[c] Sqrt[\[Pi]] Sqrt[\[Epsilon]0]))}, {R -> Sqrt[
pCos[\[Theta]]]/(2 Sqrt[c] Sqrt[\[Pi]] Sqrt[\[Epsilon]0])}}
调用适量分析包，错了，提示是Get::noopen: 无法打开 Calculus`VectorAnalysis`. >>
<<Calculus`VectorAnalysis`
Ee=Grad[-V1,Spherical[R,\[Theta],\[Omega]]]
希望大神能够帮忙一下，万分感谢

直接用Grad就行了

Ee = Grad[-V1, {R, \[Theta], \[CurlyPhi]}, "Spherical"]