N = [10 100 1000 10^4 10^6 10^6 10^6];
[one n] = size(N);
for i=1:n
  NN = N(i);
  r = rand( NN, 1);
  rsq = r.*r;
  Av = sum( rsq(1:NN,1) );
  Av = 3*Av/(2*NN);
  printLine = sprintf(' %e   %e  %e', NN, Av, Av - .5);
  disp(printLine)
  clear r
  clear rsq
end