Following the general definition of the Gaussian ensembles, for the $n\times n$ matrix $W$ I take the distribution $P(W)\propto e^{-\tfrac{1}{2}\beta\,{\rm tr}\,W^2}$, with $\beta=1$ if $W$ is real symmetric and $\beta=2$ if $W$ is complex Hermitian. The eigenvalues of $W$ are denoted by $w_1,w_2,\ldots w_n$. I consider two quantities:
$$e^{nF(w_i)}=\prod_{j=1, j\ne i}^n|w_i-w_j|\Rightarrow F(w_i)=n^{-1}\sum_{j=1,j\ne i}^n\ln|w_i-w_j|,$$
$$e^{n^2G}=\prod_{i=1}^n\prod_{j=1, j\ne i}^n|w_i-w_j|\Rightarrow G=n^{-2}\sum_{i=1}^n\sum_{j=1,j\ne i}^n\ln|w_i-w_j|.$$

For $n\gg 1$ we may replace sums by integrals, upon the replacements of $w_i\mapsto n^{1/2}\lambda(i/n)$ and of 
$$\sum_i u(w_i)\mapsto n\int dx\, u(n^{1/2}\lambda(x))=n\int d\lambda \,\rho(\lambda)u(n^{1/2}\lambda),\;\;\text{with}\;\;\rho(\lambda)=(d\lambda/dx)^{-1}.$$
We thus have the expressions
$$F(n^{1/2}\lambda)=\int \rho(\mu)d\mu\,\ln|n^{1/2}\lambda-n^{1/2}\mu|,$$
$$G=\int\rho(\lambda)d\lambda\int\rho(\mu)d\mu\,\ln |n^{1/2}\lambda-n^{1/2}\mu|.$$

The rescaled eigenvalue density $\rho(\lambda)=(d\lambda/dx)^{-1}$ is given in the large-$n$ limit by the <A HREF="https://en.wikipedia.org/wiki/Wigner_semicircle_distribution">Wigner semicircle</A>),
$$\rho(\lambda)=\frac{1}{\pi}\sqrt{2-\lambda^2},\;\;\lambda\in(-\sqrt 2,\sqrt 2).$$
We thus arrive at
$$F(n^{1/2}\lambda)=\tfrac{1}{2}\ln n+\frac{1}{\pi}\int_{-\sqrt{2}}^{\sqrt{2}}d\mu\,\sqrt{2-\mu^2}\ln|\lambda-\mu|,$$
$$G=\tfrac{1}{2}\ln n-\tfrac{1}{4} -\tfrac{1}{2}\ln 2.$$
For $\lambda=0$ the integral over $\mu$ evaluates to
$$F(0)=\tfrac{1}{2}\ln n-\tfrac{1}{2} -\tfrac{1}{2}\ln 2.$$

Now to make contact with the OP, note that
$$f(\lambda_i)^{1/n}=\prod_{j=1, j\ne i}^n|\lambda_i-\lambda_j|^{1/n}=n^{-1/2}e^{F(n^{1/2}\lambda)},$$
hence we conclude that for $n\gg 1$ we have
$$f(\lambda_i)^{1/n}=\exp\left[\frac{1}{\pi}\int_{-\sqrt{2}}^{\sqrt{2}}d\mu\,\sqrt{2-\mu^2}\ln|\lambda-\mu|\right].$$
For $\lambda_i=0$ this reduces to
$$f(0)^{1/n}=e^{-(1+\ln 2)/2}=0.428882.$$
I checked this numerically, and it agrees.