The finite element method is adopted for computation of axisymmetrical conical frusta. The normal disp1acement function ω and circumferential displacement function υ are given by two cubic po1ynomials，the longitudinal displacement function u is given by linear polynomial. Nodal degrees of freedom can extract three displacements，one rotation and one shear strain on the nodal circ1e. The computed resu1ts had compared with referencies[l]，[2] and [3]. This method and computer program are applicable to find vibration characteristics of conical frustum shell.