A theoretical formulation for the analysis of free vibration of a cylindrical shell with a circular plate attached at an arbitrary axial position of the shell under various kinds of boundary conditions was derived and programed to get the numerical results for natural frequencies and mode shapes of the combined system. The boundary conditions of the shell to be considered here are clamped-free, clamped-simply supported, both ends clamped and both ends simply supported. The frequencies and mode shapes from theoretical calculation were compared with those of commercial finite element code, ANSYS. The results showed good agreement with those of ANSYS in frequencies and mode shapes. The program will contribute to the design optimization of a shell/plate combined system through the analysis of natural frequencies and mode shapes for the system.