
new;
!las gr.out
!ri temp.out
!pr rep

load Data/GROO.out;
load Data/GRHHO.out;
load gr.out;

con = ' at T=305K Pressure = 1 bar ';
MAX = 9.15;

clf;
plot (GROO(:,1),GROO(:,2),'--');
hold on;
plot(gr(:,2),gr(:,3));
v = axis;
v(2) = MAX;
axis(v);
title(['Goo(r)' con]);
[Y,I] = max(gr(:,3));
app = ['(' num2str(gr(I,2)) ', ' num2str(Y) ')'];
xlabel(['dashed line: experiment (2.875, 3.092)   solid line: MD  ' app]);
PRINT;

clf;
plot (GRHHO(:,1),GRHHO(:,2),'--');
hold on;
plot(gr(:,2),gr(:,4));
v = axis;
v(2) = MAX;
axis(v);
title(['Goh(r)' con]);
xlabel('dashed line: experiment(1986)     solid line: new potential ')
PRINT;

clf;
plot (GRHHO(:,1),GRHHO(:,3),'--');
hold on;
plot(gr(:,2),gr(:,5));
v = axis;
v(2) = MAX;
axis(v);
title(['Ghh(r)' con]);
xlabel('dashed line: experiment(1986)     solid line: new potential ')
PRINT;
