function B = dctscale(A)
% get the dct scale factors we hope
  
  sf(1) = 1;
  for i=2:8
    sf(i) = cos((i-1)*pi/16) * sqrt(2);
  end
  
  
  for i=1:8
    for j=1:8
      B(i,j) = sf(i) * sf(j);
      
    end
  end
  