We present a vectorial formalism to determine the approximate solutions to the problem of a composite body made of L homogeneous, rigidly rotating layers bounded by spheroidal surfaces. The method is based on the 1st-order expansion of the gravitational potential over confocal parameters, thereby generalizing the method described in Paper I for L = 2. For a given relative geometry of the ellipses and a given set of mass-density jumps at the interfaces, the sequence of rotation rates and interface pressures is obtained analytically by recursion. A wide range of equilibria result when layers rotate in an asynchronous manner, although configurations with a negative oblateness gradient are more favorable. In contrast, states of global rotation (all layers move at the same rate), found by solving a linear system of L − 1 equations, are much more constrained. In this case, we mathematically demonstrate that confocal and coelliptical configurations are not permitted. Approximate formula for small ellipticities are derived. These results reinforce and prolongate known results and classical theorems restricted to small elliptiticities. Comparisons with the numerical solutions computed from the Self-Consistent-Field method are successful.