A Multiphase Shape Optimization Problem for Eigenvalues: Qualitative Study and Numerical Results