(* Li2 with external clouds, Kimball_Loebl ekin, 13.08.2012 *)
Clear[u1,u2,u3,u4, c0, c1,c2,c3,T,Vnn,Vne,Vee]; Z = 3.0;
(* Kimball's parameter, k2 lowered *)
c1={k1->1.0, k2-> 0.8, s1->0.4, s2->0.4};

Vne=-2*3*Z/Q-4*Z/(P+Q)-2*Z/(P+Q);
Vee=2*4/(P+Q)+2/(P+Q)+3*s2/P+2*3*s1/Q;
Vnn=Z*Z/(P+Q)/2;
T=2.25*k2/P^2+2*2.25*k1/Q^2;

func=T+Vne+Vee+Vnn /. c1
N[t=FindMinimum[func, {P,2.5},{Q,0.6}],10]
u2=t[[2]];
N[-(Vne+Vee+Vnn)/T /.u2 /.c1,10]
N[2*(P+Q)*0.529177 /.u2/.c1,10]
Vne /.u2/.c1
Vee /.u2 /.c1
Vnn /.u2 /.c1

plot1=Graphics[{Circle[{-(P+Q),0},Q],Circle[{P+Q,0},Q],Circle[{0,0},P],Disk[{-(P+Q),0},0.08],Disk[{(P+Q),0},0.08],{Thickness[0.01],Line[{{-(P+Q),0},{(P+Q),0}}]}}]/.u2
Show[plot1,{AspectRatio->Automatic,Axes->True,GridLines->Automatic,PlotRange->{{-5,5},{-4,4}},Frame->True}]