The spline function technique is used to analyze the vibration of multi-layered circular cylindrical shells with cross-ply walls including first-order shear deformation theory. Both antisymmetric and symmetric cross-ply laminations are considered in this analysis. The governing equilibrium equations are obtained in terms of displacement and rotational functions. A system of coupled ordinary differential equations in terms of displacement and rotational functions are obtained by assuming the solution in a separable form. These functions are approximated by using Bickley-type splines of suitable order to obtain the generalized eigenvalue problem by applying point collocation techniques with appropriate boundary conditions. Parametric studies are performed to analyze the frequency response of the shell with reference to the material properties, number of layers, fiber orientation, thickness to radius ratio, length to radius ratio and circumferential node number. Reasonable agreement is found with existing results obtained by FEM and other methods. Valuable results are presented as graphs and discussed.