clf;

GSS_A = -161.8414350882-459.4479639194;
MP2_A = -161.8414350882-459.5524332776;

load GSS.OUT;
load MP2.OUT;

plot (MP2(:,1),MP2(:,2)-MP2_A);
hold on;

plot (GSS(:,1),GSS(:,2)-GSS_A,'--');
axis ([1 7 -0.155 0.35]);

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

v = axis;
txt = ['\10 Solid line = Moller-Plesset 2nd-order PT'];
stext (v(1)*0.55+v(2)*0.45,v(3)*0.1+v(4)*0.9,txt);
txt = ['\10 Dashed line = 6-31G^{**}  UHF Mult=1 '];
stext (v(1)*0.55+v(2)*0.45,v(3)*0.16+v(4)*0.84,txt);

