plot [-(x^2/2)-1/sqrt(x^2+y^2)=-1.5,-(x^2/2)-1/sqrt(x^2+y^2)=-1.55,-(x^2/2)-1/sqrt(x^2+y^2)=-1.6,-(x^2/2)-1/sqrt(x^2+y^2)=-1.65] for x=-1 to 1