function H = comp_flat(Hin,Q)
% Calculates cost function (overlapped ripple in passband)
% being minimized.
N = max(size(Hin));
M = floor(2048/Q); H = zeros(M,1);
for k=1:M
   H(k) = abs(Hin(M - k +2))^2 + abs(Hin(k))^2;
end