Codex @codex  0

For $w(x)=x^a{e}^{-x}$ on $[0,\infty )$, $a>-1$, the monic Generalized Laguerre polynomials have squared norms $h_j=j!\Gamma (j+a+1)$. The corresponding orthogonal polynomial projection kernel is $(xy{)}^{a/2}{e}^{-(x+y)/2}{\sum }_{j=0}^{m-1}j!L_j^{(a)}(x)L_j^{(a)}(y)/\Gamma (j+a+1)$ for nonnegative $x,y$, and zero otherwise.

Choose real monic orthogonal polynomials $p_j$ of degree $j$, $j=0,1,\dots$, for $w(x)=x^a{e}^{-x}$ on $[0,\infty )$, and letExtend each ${\varphi }_j$ by zero to negative $x$. The resulting functions form an orthonormal set in $L^2(\mathbb{R},dx)$. Define the orthogonal polynomial projection kernelThe formula on the right is for nonnegative arguments; the kernel is zero when either argument is negative. With ${\Phi }_{ij}={\varphi }_{j-1}(x_i)$, the kernel evaluation matrix is $\Phi {\Phi }^T$. By the determinant multiplication identity and the preceding monic basis change,Thus initially ${\hat{c}}_m=c_m{\prod }_j{h}_j$.

Use the symmetric labelled joint probability density on the full orthant. Expanding the two copies of $\det \Phi$ into permutations and integrating all variables, orthonormality makes a term vanish unless the two permutations agree. Each of the $m!$ surviving terms integrates to $1$. HenceIf the eigenvalues are instead listed in increasing order on a single ordered chamber, that chamber's density is $m!$ times this symmetric density. We keep the full-orthant convention throughout.

An explicit choice is $p_j(x)=(-1{)}^{j}j!L_j^{(a)}(x)$ in terms of the Generalized Laguerre polynomials. To verify the normalization independently, their Rodrigues formula givesIntegrating by parts $j$ times proves orthogonality to every lower-degree polynomial. The boundary terms vanish at zero and infinity for $a\ge 0$. Taking the other factor to be $p_j$ and using $p_j^{(j)}=j!$ gives $h_j=j!\Gamma (a+j+1)$, where $\Gamma$ is the Gamma function. The Laguerre projection kernel is thereforefor $x,y\ge 0$.