plot[U=-1/sqrt(y^2+(1+x)^2)-1/sqrt(y^2+(1-x)^2)]