def hermpow(c, pow, maxpower=<NUM_LIT:16>):
<EOL>[c] = pu.as_series([c])<EOL>power = int(pow)<EOL>if power != pow or power < <NUM_LIT:0>:<EOL><INDENT>raise ValueError("<STR_LIT>")<EOL><DEDENT>elif maxpower is not None and power > maxpower:<EOL><INDENT>raise ValueError("<STR_LIT>")<EOL><DEDENT>elif power == <NUM_LIT:0>:<EOL><INDENT>return np.array([<NUM_LIT:1>], dtype=c.dtype)<EOL><DEDENT>elif power == <NUM_LIT:1>:<EOL><INDENT>return c<EOL><DEDENT>else:<EOL><INDENT>prd = c<EOL>for i in range(<NUM_LIT:2>, power + <NUM_LIT:1>):<EOL><INDENT>prd = hermmul(prd, c)<EOL><DEDENT>return prd<EOL><DEDENT>
Raise a Hermite series to a power.

    Returns the Hermite series `c` raised to the power `pow`. The
    argument `c` is a sequence of coefficients ordered from low to high.
    i.e., [1,2,3] is the series  ``P_0 + 2*P_1 + 3*P_2.``

    Parameters
    ----------
    c : array_like
        1-D array of Hermite series coefficients ordered from low to
        high.
    pow : integer
        Power to which the series will be raised
    maxpower : integer, optional
        Maximum power allowed. This is mainly to limit growth of the series
        to unmanageable size. Default is 16

    Returns
    -------
    coef : ndarray
        Hermite series of power.

    See Also
    --------
    hermadd, hermsub, hermmul, hermdiv

    Examples
    --------
    >>> from numpy.polynomial.hermite import hermpow
    >>> hermpow([1, 2, 3], 2)
    array([ 81.,  52.,  82.,  12.,   9.])