An eigenvalue analysis of the circular Mindlin plates with free boundary conditions is presented. The analysis is based on the Chebyshev-Fourier pseudospectral method. Even though the eigenvalues of lower vibration modes tend to convergence more slowly than those of higher vibration modes, the eigenvalues converge for sufficiently fine pseudospectral grid resolutions. The eigenvalues of the axisymmetric modes are computed separately. Numerical results are provided for different grid resolutions and for different thickness-to-radius ratios.