
clear;

load GROO.out;
load pam.out;
A = pam(1);
B = pam(2);
C = pam(3);

START = 35;
clf;

aa = size(GROO);
bb = aa(1);
for i=START:bb
r1 = GROO(i,1);
r2 = r1*r1;
r6 = r2*r2*r2;
V(i-START+1) = A*exp(-B*r1)+C/r6;
V1(i-START+1) = 1895.294/r6/r6-1.88354/r6;
x(i-START+1) = r1;
end;

plot(x,V1);
hold on;
plot(x,V,'--');
axis ([1 10 -0.001 0.001])


