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

data params;
  input a alpha ratiovar;
cards;
  5 .05 2.0
;

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

proc print;
var n beta;
run;