1.  A gamma variate X with shape parameter ALPHA and scale BETA can be generated:

 x=beta*rangam(seed,alpha);

2. you can use the following program to  generate  random variables for gamma distribution.

data test;
do i=1 to 10000;
      X=rangam(20,2);        /*20 means random seed*/   /*2 is a shape parameter*/
      output;                           
end;
run;


3.  goodness of fit test for gamma distribution

proc capability data=test;
      var x;
      histogram /          gamma   ;
run;