data power;
file 'z:\stat706\power.txt';
alpha=.05;
a=5;
am1=a-1;
do s02=1 to 5;
  do n=2 to 15;
    nc=n*s02;
    nma=a*(n-1);
    fcrit=finv(1-alpha,am1,nma);
    power=1-probf(fcrit,am1,nma,nc);
    put n s02 power;
  end;
end;
run;