human
I wrote the below code

    N = len(cnt)
    if (N < 2):
        raise ValueError('At    least    two    moments    must    be    given    to    approximate    the    pdf.')
    (mc, mc2, mc3, mc4) = cnt
    skew = (mc3 / (mc2 ** 1.5))
    kurt = ((mc4 / (mc2 ** 2.0)) - 3.0)
    totp = poly1d(1)
    sig = sqrt(cnt[1])
    mu = cnt[0]
    if (N > 2):
        Dvals = _hermnorm((N + 1))
        C3 = (skew / 6.0)
        C4 = (kurt / 24.0)
        totp = ((totp - (C3 * Dvals[3])) + (C4 * Dvals[4]))
    def thisfunc(x):
        xn = ((x - mu) / sig)
        return (((totp(xn) * np.exp((((- xn) * xn) / 2.0))) / np.sqrt((2 * np.pi))) / sig)
    return thisfunc

What's a good function header?
 ... Here's a good function header?
assistant
#Here's a good function header.

def pdf_moments(cnt):