clear all
clc
B=[2000];
n=10;
m = size(B);

x = rand(1,n);


mean1=mean(x);
median1 = median(x);


for i=1:m
	for j=1:B(i)
		u1 = randint(1,n,[1,n]);	
		for k=1:n
			x2(k,:) = x(u1(k));
		end
		mean2(j)=mean(x2);
		median2(j)=median(x2);

    end
    mean1
    
%     mean2 = mean2(end)
%     median2 = median2(end)
	mean_bt=mean(x2)
    mean_bais = mean_bt - mean1
    median1
	median_bt=median(x2)
	median_bais = median_bt - median1
end