% T=300K, a=10 [b], K/2e=0.1 in a 41960 [b] (10 um) by 10240 [b] cell.

clear;
% K/2e, stress [e/b^2]
cs = [
0.0001 0.45E-4 
0.0005 1.05E-4 
0.001  1.65E-4 
0.002  2.50E-4 
0.003  3.05E-4 
0.005  3.95E-4 
0.01   5.75E-4
];

clf;
loglog (cs(:,1), cs(:,2), 'o');
xlabel('Obstacle volume ratio c');
ylabel('2\sigma_c^{sim}/\mu(T)');
title('Verification of Friedel \sigma_c \propto c^{1/2} scaling law for Cu at T=300K, a=10, K/2\epsilon=0.1');
cc = [0.0001 0.01];
hold on; plot (cc, 2./sqrt(100./cc)*(0.1^1.5), 'r--');
text (0.00011, 0.00013, 'theoretical prediction');
axis ([0.0001 0.01 4E-5 1E-3]);
print -depsc c.eps
