code: def _roots_hermite_asy(n):
      r"""Gauss-Hermite (physicst's) quadrature for large n.
  
      Computes the sample points and weights for Gauss-Hermite quadrature.
      The sample points are the roots of the n-th degree Hermite polynomial,
      :math:`H_n(x)`.  These sample points and weights correctly integrate
      polynomials of degree :math:`2n - 1` or less over the interval
      :math:`[-\infty, \infty]` with weight function :math:`f(x) = e^{-x^2}`.
  
      This method relies on asymptotic expansions which work best for n > 150.
      The algorithm has linear runtime making computation for very large n
      feasible.
  
      Parameters
      ----------
      n : int
          quadrature order
  
      Returns
      -------
      nodes : ndarray
          Quadrature nodes
      weights : ndarray
          Quadrature weights
  
      See Also
      --------
      roots_hermite
  
      References
      ----------
      .. [townsend.trogdon.olver-2014]
         Townsend, A. and Trogdon, T. and Olver, S. (2014)
         *Fast computation of Gauss quadrature nodes and
         weights on the whole real line*. :arXiv:`1410.5286`.
  
      .. [townsend.trogdon.olver-2015]
         Townsend, A. and Trogdon, T. and Olver, S. (2015)
         *Fast computation of Gauss quadrature nodes and
         weights on the whole real line*.
         IMA Journal of Numerical Analysis
         :doi:`10.1093/imanum/drv002`.
      """
      iv = _initial_nodes(n)
      nodes, weights = _newton(n, iv)
      # Combine with negative parts
      if n % 2 == 0:
          nodes = hstack([-nodes[::-1], nodes])
          weights = hstack([weights[::-1], weights])
      else:
          nodes = hstack([-nodes[-1:0:-1], nodes])
          weights = hstack([weights[-1:0:-1], weights])
      # Scale weights
      weights *= sqrt(pi) / sum(weights)
      return nodes, weights
  
imports: from __future__ import division, print_function, absolute_import
  import numpy
  from numpy import exp, inf, pi, sqrt, floor, sin, cos, around, int, hstack, arccos, arange
  from scipy import linalg
  from scipy.special import airy
  from .None import _ufuncs
  from .None import _ufuncs
  from .None import specfun
  from ._ufuncs import binom, eval_jacobi, eval_sh_jacobi, eval_gegenbauer, eval_chebyt, eval_chebyu, eval_chebys, eval_chebyc, eval_sh_chebyt, eval_sh_chebyu, eval_legendre, eval_sh_legendre, eval_genlaguerre, eval_laguerre, eval_hermite, eval_hermitenorm
label_desc: Use image node_id 17 for calling a global function with example usage: _roots_hermite_asy(n) and returns: nodes,  weights