function [Me,v] = BM_Average(m)
tic
data = zeros([m,1]);


for i = 1:m
    data(i) = BM(1000,100000,0);
end

Me=mean(data)
v=var(data)/m
toc