/* module 8b: */

/* 8-b: 1 */

data one;
do j = 1 to 8;
  if j = 1 then samples = 2;
  if j = 2 then samples = 5;
  if j = 3 then samples = 10;
  if j = 4 then samples = 20;
  if j = 5 then samples = 50;
  if j = 6 then samples = 70;
  if j = 7 then samples = 100;
  if j = 8 then samples = 500;
  do i = 1 to 1000;
    z = 0;
    do k = 1 to samples;
      x = ranexp(37);

      z = z + (x / samples);
      z2 = z - 1;
    end;

    output;
  end;
end;
run;
proc chart data=one;
  by samples;
*  vbar z;
  vbar z2;
run;
proc means data=one;
  by samples;
run;