(*Three body problem*)

(* PRIMARY OBJECTS WITH EQUAL MASSES *)

v[x_,y_,α_]=-((1-α)/(Sqrt[(((x-α))^2)+(y^2)]))-(α/(Sqrt[(((x+1-α))^2)+(y^2)]))-(1/2)((x^2)+(y^2));

Plot3D[v[x,y,0.5],{x,-1.5,1.5},{y,-1.5,1.5}]

ContourPlot[v[x,y,0.5],{x,-1.5,1.5},{y,-1.5,1.5},Contours-> 50]

(* SUN-JUPITER SYSTEM: α=0.000953875 *)

Plot3D[v[x,y,0.1],{x,-1.5,1.5},{y,-1.5,1.5}]

ContourPlot[v[x,y,0.1],{x,-1.5,1.5},{y,-1.5,1.5},Contours-> 40]

ContourPlot[v[x,y,0.1],{x,-1,0},{y,0,1},Contours-> 300]