options nocenter ps=35 ls=72;
/* This is how you can calculate the power in one data set */

data params;
  input a alpha d var;
cards;
  5 .05 3.0 1.8
;

data new;
 set params;
 do n=2 to 10;
    df = a*(n-1);
    nc = n*d*d/(2*var);
    fcut  = finv(1-alpha,a-1,df);
    beta=probf(fcut,a-1,df,nc);
    output;
 end;

proc print;
var n nc beta;
run;