human
Complete this python program:
def _norm_factor(p, k):
assistant
#Here is the Python code you requested.

    """Numerically find frequency shift to apply to delay-normalized filter such
    that -3 dB point is at 1 rad/sec.
    `p` is an array_like of polynomial poles
    `k` is a float gain
    First 10 values are listed in "Bessel Scale Factors" table,
    "Bessel Filters Polynomials, Poles and Circuit Elements 2003, C. Bond."
    """
    p = asarray(p, dtype=complex)
    def G(w):
        '\n                                Gain    of    filter\n                                '
        return abs((k / prod(((1j * w) - p))))
    def cutoff(w):
        '\n                                When    gain    =    -3    dB,    return    0\n                                '
        return (G(w) - (1 / np.sqrt(2)))
    return optimize.newton(cutoff, 1.5)