
load GROO.out;
load pam.out;
A = pam(1);
B = pam(2);
C = pam(3);
	
clf;
plot (GROO(:,1),GROO(:,2),'--');
hold on;

beta = 1./(305.*1.380658E-23/2.3071E-18);

aa = size(GROO);
bb = aa(1);
for i=40:bb
r1 = GROO(i,1);
r2 = r1*r1;
r6 = r2*r2*r2;
V(i-39) = A*exp(-B*r1)+C/r6; 
x(i-39) = r1;
y(i-39) = exp(-beta*V(i-39));
end;

plot(x,y);


