This modifies an input parameter. ``k = k * real(prod(2*fs - z) / prod(2*fs - p))`` would be better.