#NEWLINE# def ellip_normal( h2 , k2 , n , p ) : #NEWLINE# #INDENT# '\n    Ellipsoidal harmonic normalization constants gamma^p_n\n\n    The normalization constant is defined as\n\n    .. math::\n\n       \\gamma^p_n=8\\int_{0}^{h}dx\\int_{h}^{k}dy\\frac{(y^2-x^2)(E^p_n(y)E^p_n(x))^2}{\\sqrt((k^2-y^2)(y^2-h^2)(h^2-x^2)(k^2-x^2)}\n\n    Parameters\n    ----------\n    h2 : float\n        ``h**2``\n    k2 : float\n        ``k**2``; should be larger than ``h**2``\n    n : int\n        Degree.\n    p : int\n        Order, can range between [1,2n+1].\n\n    Returns\n    -------\n    gamma : float\n        The normalization constant :math:`\\gamma^p_n`\n\n    See Also\n    --------\n    ellip_harm, ellip_harm_2\n\n    Notes\n    -----\n    .. versionadded:: 0.15.0\n\n    Examples\n    --------\n    >>> from scipy.special import ellip_normal\n    >>> w = ellip_normal(5,8,3,7)\n    >>> w\n    1723.38796997\n\n    ' #NEWLINE# with _ellip_lock : #NEWLINE# #INDENT# with np . errstate ( all = 'ignore' ) : #NEWLINE# #INDENT# return _ellip_normal_vec ( h2 , k2 , n , h2 )