plot z=1/((x)^2+(y-1)^2)+1/((x-sqrt(3)/2)^2+(y+1/2)^2)+1/((x+sqrt(3)/2)^2+(y+1/2)^2)