%test

%SCREEN OUTPUT: Plot of PVal against D

plot(D,PVAL,'k.','MarkerSize',20);

hold on;

plot(D,PVAL,'k-','LineWidth',2);

hold on;

SIG = .05*ones(length(D),1);

plot(D,SIG,'r--');

plot(D,1 - SIG,'r--');


xlabel('Radius','Fontsize',12, 'FontWeight', 'bold');

ylabel('attraction                   P-Value                    repulsion',...
   'Fontsize',12, 'FontWeight', 'bold');