
V[x_,y_]:=Module[{r},r=Sqrt[x^2+y^2];If[r>1,Return[(1/r^2-r)*Re[x+I y]/r],Return[0]]]
			

Ex[x_,y_]:=Module[{r,cos,sin,Er,Etheta},
	r=Sqrt[x^2+y^2];cos=Re[x+I y]/r;sin=Im[x+I y]/r;
	Er=(2/r^3+1)*cos;
	Etheta=(1/r^3-1)*sin;
	If [r>1, Return[Er*cos-Etheta*sin],Return[0]];
]
Ey[x_,y_]:=Module[{r,cos,sin,Er,Etheta},
	r=Sqrt[x^2+y^2];cos=Re[x+I y]/r;sin=Im[x+I y]/r;
	Er=(2/r^3+1)*cos;
	Etheta=(1/r^3-1)*sin;
	If [r>1, Return[Er*sin+Etheta*cos],Return[0]];
]
			

p1=ContourPlot[V[x,y],{x,-5,5},{y,-5,5},PlotPoints->200,Contours->30,ContourShading->None,ContourStyle->Red]
			

p2=Show[StreamPlot[{Ex[x,y],Ey[x,y]},{x,-5,5},{y,-5,5},StreamPoints->100],Graphics[Circle[{0,0}]]]
			
Show[p1,p2]
