After introducing the Gaudin spin chain as a quantum integrable system and the Bethe Ansatz as applied to this system, I will describe some recent work by A.Bourget, D.Jakobson, J.Toth and myself about the average distribution of zeros of Heine-Stieltjes polynomials in the thermodynamic limit.