Nstat=10;
bpd=2;

srand=dir_srand;
k1=full(sum(srand'));
k1_max=max(k1);

k2=full(sum(srand));
k2_max=max(k2);

b1_max=1+floor(bpd.*log10(k1_max));
b2_max=1+floor(bpd.*log10(k2_max));

n_1_2_asym=zeros(b1_max+1,b2_max+1);
binned_srand;

n_1_2_asym_orig=n_1_2;
dir_generate_srand;

binned_srand;

aver_n_1_2_asym=n_1_2;
aver_sq_n_1_2_asym=n_1_2.^2;

for k=2:Nstat;
   k
   dir_generate_srand;
   binned_srand;
   aver_n_1_2_asym=aver_n_1_2_asym+n_1_2;
   aver_sq_n_1_2_asym=aver_sq_n_1_2_asym+n_1_2.^2;
end;
aver_n_1_2_asym=aver_n_1_2_asym./Nstat;
aver_sq_n_1_2_asym=aver_sq_n_1_2_asym./Nstat;
err_n_1_2_asym=sqrt(aver_sq_n_1_2_asym-aver_n_1_2_asym.^2);


dir_ratio_1_2_asym=n_1_2_asym_orig./(aver_n_1_2_asym+0.0001.*(aver_n_1_2_asym==0));
dev_n_1_2_asym_orig=(n_1_2_asym_orig-aver_n_1_2_asym)./(err_n_1_2_asym+0.0001.*(aver_n_1_2_asym==0));


dir_ratio_1_2_asym(end,:)=dir_ratio_1_2_asym(end-1,:);
dir_ratio_1_2_asym(:,end)=dir_ratio_1_2_asym(:,end-1);
dev_n_1_2_asym_orig(end,:)=dev_n_1_2_asym_orig(end-1,:);
dev_n_1_2_asym_orig(:,end)=dev_n_1_2_asym_orig(:,end-1);

figure;  pcolor(dir_ratio_1_2_asym);  colorbar


shading interp; 
colorbar;

figure;  pcolor(dev_n_1_2_asym_orig);  



shading interp; 
colorbar;

%shading interp; colormap hot; 
num2str(dir_ratio_1_2_asym(1:end-1,1:end-1))
num2str(dev_n_1_2_asym_orig(1:end-1,1:end-1))