clf;

STO_A = -126.6045250888;
GSS_A = -128.5226687496;
MP2_A = -128.7392465733;

load GSS.OUT;
load STO.OUT;
load MP2.OUT;
load CI3.OUT;

plot (MP2(:,1),MP2(:,2)-2*MP2_A);
hold on;
axis ([2.2,3.6,-0.00026 0.0012]);

plot (GSS(:,1),GSS(:,2)-2*GSS_A,'--');

plot (STO(:,1),STO(:,2)-2*STO_A,':');

a = size(CI3);
plot (CI3(:,1),CI3(:,2)-CI3(a(1,1),2),'o');

ylabel ('a.u.');
xlabel ('Angstrom');
title ('Neon-Neon Binding Curve');

v = axis;
txt = ['\10 Dotted line = STO-3G Minimal Basis'];
stext (v(1)*0.55+v(2)*0.45,v(3)*0.1+v(4)*0.9,txt);
txt = ['\10 Dashed line = 6-311G^{**}  RHF '];
stext (v(1)*0.55+v(2)*0.45,v(3)*0.16+v(4)*0.84,txt);
txt = ['\10 Solid line = 6-311G^{**} MP2'];
stext (v(1)*0.55+v(2)*0.45,v(3)*0.22+v(4)*0.78,txt);
txt = ['\10 Circle = TCI using 6-311G^{**} basis'];
stext (v(1)*0.55+v(2)*0.45,v(3)*0.28+v(4)*0.72,txt);

ss=input('Continue... ','s');
clf;
S = MP2(:,2)-2*MP2_A-GSS(:,2)+2*GSS_A;
plot (MP2(:,1),S,'--');
axis ([2.2,3.6,-0.00075 0.0015]);

txt = ['\10 Dashed line = MP2-RHF (6-311G^{**})'];
stext (v(1)*0.55+v(2)*0.45,v(3)*0.16+v(4)*0.84,txt);

