sumGM = zeros(N,1); % denominator