Let G be an N×N real matrix whose entries are independent identically distributed standard normal random variables Gij∼𝒩(0,1). The eigenvalues of such matrices are known to form a two-component system consisting of purely real and complex conjugated points. The purpose of this paper is to show that by appropriately adapting the methods of [E. Kanzieper, M. Poplavskyi, C. Timm, R. Tribe and O. Zaboronski, Annals of Applied Probability26(5) (2016) 2733–2753], we can prove a central limit theorem of the following form: if λ1,…,λNℝ are the real eigenvalues of G, then for any even polynomial function P(x) and even N=2n, we have the convergence in distribution to a normal random variable
1√𝔼(Nℝ)(Nℝ∑j=1P(λj/√2n)−𝔼Nℝ∑j=1P(λj/√2n))→𝒩(0,σ2(P))
as n→∞, where σ2(P)=2−√22∫1−1P(x)2dx.