options ls=100;

data one;
   array alpha a1-a8 (0.2,0.15,0.1,0.05,0.025,0.01,0.005,0.001);
   array tstar t1-t8; 
   do df=1 to 30,40 to 100 by 10,200,100000;
      do i=1 to 8;
         tstar{i}=tinv(1-alpha{i},df);
      end;
      output;
   end;
run;

proc print;
   var df t1-t8;
run;