A method is presented to study the free vibrations of rectangular laminated composite plates with general layups and arbitrary boundary conditions. Based on the first-order shear deformation theory, the governing differential equations and boundary conditions are deduced via Hamilton’s principle. Generalised displacements are expanded as series with Legendre polynomials as the base functions. A generalised eigenvalue problem is obtained by following a variational approach, where energy functional is extremised and boundary conditions are introduced by means of Lagrange multipliers. In order to overcome some difficulties in obtaining the natural frequencies and corresponding mode shapes, a new numerical strategy is proposed.