A computational method for studying molecular vibrations and spectra for symmetrical systems with many degrees of freedom was developed. The algorithm allows overcoming difficulties on the automation of calculus related to the symmetry determination of such oscillations in complex systems with many degrees of freedom. One can find symmetrized displacements and, consequently, obtain and classify normal oscillations and their frequencies. The problem is therefore reduced to the determination of eigenvectors by common numerical methods, and the algorithm simplifies the procedure of symmetry determination for normal oscillations. The proposed method was applied to studying molecular vibrations and spectra of the fullerene molecule C₆₀, and the comparison of theoretical results with experimental data is drawn. The computational method can be further extended to other problems of group theory in physics with applications in clusters and nanostructured materials.