Integrable model for interacting electrons in metallic grains

We find an integrable generalization of the BCS model with nonuniform Coulomb and pairing interaction. The Hamiltonian is integrable by construction since it is a functional of commuting operators; these operators, which therefore are constants of motion of the model, contain the anisotropic Gaudin Hamiltonians. The exact solution is obtained diagonalizing them by means of Bethe ansatz. Uniform pairing and Coulomb interaction are obtained as the "isotropic limit" of the Gaudin Hamiltonians. We discuss possible applications of this model to a single grain and to a system of few interacting grains.