%
%
%
% Optical Imaging and Spectroscopy
%
% David J. Brady
% Duke University
% www.opticalimaging.org
%
% 
%
function qc = QQ(u,m)
%Riesz denominator 
qc=u;
for pip=1:length(qc)
    qc(pip)=1./sum(sinc(qc(pip)+((-1000):1000)).^(2*(m+1)));
end
return