plot 1/sqrt((x-1)^2+(y-1)^2)+1/sqrt((x-1)^2+(y+1)^2)+1/sqrt((x+1)^2+(y-1)^2)+1/sqrt((x+1)^2+(y+1)^2), at -1<x<1, y=0