This paper relates to the flexural vibration analysis of slender spindle systems with multiple thin disks, supported by the ball bearings by means of the finite element method. Each system component is analytically modeled taking into account its flexibility and also the centrifugal effect especially for the disk. In order to show the rapid convergence rate and accuracy of the proposed approach, an experimental set-up is built to be versatile. In two distinct cases, its natural modes are numerically computed using only a small number of total element meshes as the shaft rotational speed is varied, and verified through experimental frequency response function obtained by the impact test.