Orthogonal and symplectic matrix models with real analytic potentials and multi interval supports of the equilibrium measures will be discussed. For these models universality of local eigenvalue statistics and bounds for the rate of convergence of linear eigenvalue statistics and for the variance of linear eigenvalue statistics are obtained.