% T=300K, a=10 [b], c=0.001 in a 41960 [b] (10 um) by 10240 [b] cell.

clear;
% K/2e, stress [e/b^2]
ks = [
0.05 0.57E-4 
0.1  1.65E-4 
0.2  4.85E-4
0.3  8.45E-4
0.4  1.50E-3
];

clf;
loglog (ks(:,1), ks(:,2), 'o');
xlabel('Obstacle strength K/2\epsilon');
ylabel('2\sigma_c^{sim}/\mu(T)');
title('Verification of Friedel \sigma_c \propto (K/2\epsilon)^{3/2} scaling law for Cu at T=300K, a=10, c=0.001');
L = sqrt(41960*10240/4194);
kk = [0.01 1];
hold on; plot (kk, 2/L*kk.^1.5, 'r--');
text (0.0105, 3.5E-5, 'theoretical prediction');
axis ([0.01 1 0.1E-4 1E-2]);
print -depsc K.eps
