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