(* Model V(r) = -1/r + Lambda[1-exp(-Delta r)]/r *)
(*SetDirectory["~/math/screened"];*)
lamin=0;
lamax=2;
emin=-0.05;
emax=0.5;
tabd={0.01,0.2,0.4,0.6,0.8,1,1.2,1.5,2};
o="sum.dat";
elist=ReadList[o,{Number,Number,Number,Number,Number}];
Dimensions[elist];
npoints=Dimensions[elist][[1]];
md=Length[tabd];
Clear[redat,imdat,f,la,pl,plo];
Do[
delta=tabd[[nd]];
redat[nd]=imdat[nd]={};
Do[{lambdan,deltan,re,im,acc}=elist[[n]];
   If[deltan!=delta,Goto[1]];
   If[acc<3,Return[]];
   lamaxn=lambdan;
   redat[nd]=Append[redat[nd],{lambdan,-re}];
   Label[1],{n,npoints}];
f[nd]=Interpolation[redat[nd]];
pl[nd]=Plot[f[nd][la],{la,lamin,lamaxn},
PlotStyle->{{Thickness[0.003]}}]
,{nd,md}];
plinf=Plot[(1-la)^2/2,{la,lamin,1},
PlotStyle->{{Thickness[0.005]}}];
plo=Apply[Sequence,Array[pl,md]];
   (* plotting together *)                                                     

Clear[unit,leng,thic,xmin,xmax];
ticks[unit_,leng_,thic_][xmin_,xmax_]:=
   Block[{x,x0},
   x0=unit Ceiling[xmin/unit];
   Table[{x,x,{leng,0},{GrayLevel[0.],Thickness[thic]}},{x,x0,xmax,unit}]];
plt = 
   Show[plo, plinf,Graphics[{
Text[StyleForm[
    "\!\(V \((r)\)\  = 
    \(-\ \(1\/r\)\) + \ \(\[Lambda]\/r\) \((1 - e\^\(-\[Delta]r\))\)\)", 
    FontSize -> 14],{1.5,0.45},{0,0}],
Text[StyleForm[
    SequenceForm[\[Delta]," = ",tabd], 
    FontSize -> 12],{1.5,0.38},{0,0}]
}],
    TextStyle -> {FontFamily -> "Times", FontSize -> 24}, 
    (*FormatType -> TraditionalForm, *)
    AxesLabel -> {\[Lambda], Subscript[E, I]}, 
    AxesStyle -> {Thickness[0.004]}, AspectRatio -> 10/7.5, 
    PlotRange -> {{lamin, lamax}, {emin, emax}},
    Ticks->{ticks[0.5,0.015,0.004],ticks[0.1,0.015,0.004]}];
Display["plot2.gif",plt,"GIF",ImageSize->72 {15/2,10}];