n = 50;
r = rand(n,1);
plot(r)

 m = mean(r);
hold on
plot([0,n],[m, m])
hold off
title('Mean of Random Uniform Data')

nsamples = 5;
npoints = 50;

for k = 1:nsamples
   iterationString = ['Iteration #',int2str(k)];
   disp(iterationString)
   currentData = rand(npoints,1);
   sampleMean(k) = mean(currentData)
end
overallMean = mean(sampleMean);

if overallMean <.49
    disp('Mean is less than expected');
elseif overallMean > 0.51
    disp('Mean is greater thean expected');
else
    disp('Mean is within the expected range')
end
Iteration #1

sampleMean =

    0.5010    0.5636    0.4783    0.5249    0.4814

Iteration #2

sampleMean =

    0.5010    0.5130    0.4783    0.5249    0.4814

Iteration #3

sampleMean =

    0.5010    0.5130    0.4967    0.5249    0.4814

Iteration #4

sampleMean =

    0.5010    0.5130    0.4967    0.4758    0.4814

Iteration #5

sampleMean =

    0.5010    0.5130    0.4967    0.4758    0.5272

Mean is within the expected range