Symplectic aspects of Gaudin integrable systems
and Szegö kernel variational method