data set1;
  do i = 1 to 300000;
    score = 1000*ranexp(290204920);
	output;
  end;
run;


proc univariate;
  var score;
  histogram /midpoints =(500 to 7500 by 1000) cfill=red vaxis = axis1;
  histogram /midpoints =(250 to 7750 by 500) cfill=red vaxis = axis1;
  histogram /midpoints =(100 to 7900 by 200)  cfill = red vaxis = axis1;
  histogram /midpoints =(10 to 7990 by 20)  cfill = red vaxis = axis1;
  histogram /midpoints =(250 to 7750 by 500)   normal vaxis = axis1;
  histogram /midpoints =(250 to 7750 by 500)   normal(fill percents = 10 to 90 by 10) nobars vaxis = axis1;
  axis1 value = (height = 2) label = (height=1.5);
  axis2 value = (height = 2) label = (height=1.5) major=none;
run;

* Select a subset.  Choosing only 300;
data set1;
  set set1;
  if i = 1000*(int(i/1000));
run;


proc univariate;
  var score;
  histogram /midpoints =(500 to 7500 by 1000) cfill=red vaxis = axis1;
  histogram /midpoints =(250 to 7750 by 500) cfill=red vaxis = axis1;
  histogram /midpoints =(100 to 7900 by 200)  cfill = red vaxis = axis1;
  histogram /midpoints =(10 to 7990 by 20)  cfill = red vaxis = axis1;
  histogram /midpoints =(250 to 7750 by 500)   normal vaxis = axis1;
  histogram /midpoints =(250 to 7750 by 500)   exponential vaxis = axis1;
  histogram /midpoints =(250 to 7750 by 500)   normal(fill percents = 10 to 90 by 10) nobars vaxis = axis1;
  histogram /midpoints =(250 to 7750 by 500)   exponential(fill percents = 10 to 90 by 10) nobars vaxis = axis1;
  axis1 value = (height = 2) label = (height=1.5);
  axis2 value = (height = 2) label = (height=1.5) major=none;
run;