options ls=90 nocenter;
%let seed=12345;
data one;
   array vmeans{10} (28,  30,30,30,30,30,  33,33,33,33);
   do block=1 to 4;
      B=3*rannor(&seed);
      do variety = 1 to 10;
         yield = vmeans{variety} + B + rannor(&seed);
         output;
      end;
   end;
run;
proc print data=one noobs;
   var block variety yield;
run;
proc glm data=one;
   class block variety;
   model yield = block variety;
   means variety/tukey lsd;
   means variety;
   ods output means=means;
run;
proc sort data=means;
   by mean_yield;
run;
data means; set means; amper="&"; eol = " \\"; mean_yield=round(mean_yield,.01);
proc print data=means noobs;
   var variety amper mean_yield amper amper eol;
run;