Verallgemeinertes Potential

U[X_, Y_]:=-  ((1-α)/(Sqrt[Y^2+(X+α)^2]))- (α/(Sqrt[Y^2+(X-1+α)^2]))-1/2  (X^2+Y^2)

α = 1/4;
range = 1.4;
ContourPlot[U[X,Y],{X,-range,range},{Y,-range,range}, Contours->50]

Plot3D[U[X,Y],{X,-range,range},{Y,-range,range}]

Plot[U[X,0],{X,-1.5range,1.5range}]