gm = zeros(K,N); % gaussian component in the nominator