code: def qnwgamma(n, a=1.0, b=1.0, tol=3e-14):
      
      return _make_multidim_func(_qnwgamma1, n, a, b, tol)
docs: Computes nodes and weights for gamma distribution
  
      Parameters
      ----------
      n : int or array_like(float)
          A length-d iterable of the number of nodes in each dimension
  
      a : scalar or array_like(float) : optional(default=ones(d))
          Shape parameter of the gamma distribution parameter. Must be positive
  
      b : scalar or array_like(float) : optional(default=ones(d))
          Scale parameter of the gamma distribution parameter. Must be positive
  
      tol : scalar or array_like(float) : optional(default=ones(d) * 3e-14)
          Tolerance parameter for newton iterations for each node
  
      Returns
      -------
      nodes : np.ndarray(dtype=float)
          Quadrature nodes
  
      weights : np.ndarray(dtype=float)
          Weights for quadrature nodes
  
      Notes
      -----
      Based of original function ``qnwgamma`` in CompEcon toolbox by
      Miranda and Fackler
  
      References
      ----------
      Miranda, Mario J, and Paul L Fackler. Applied Computational
      Economics and Finance, MIT Press, 2002.