Flexural vibration in a finite piezoelectric circular solid cylinder of crystal class 6mm is studied for traction-free and electric potential-free boundaries. The solution is obtained in the form of a series, in which the boundary conditions for the shear stresses on curved and plane end surfaces and electric potential on the plane ends are satisfied term-by-term. Remaining other boundary conditions are satisfied by an orthogonalization procedure. The series converges rapidly within a few terms. Numerical computation is carried out for the piezoelectric ceramic PZT4. For various ratios of half-height to radius of the cylinder, the natural frequencies are calculated for the first and second circular modes. The stress distribution over the cylinder for the first circular mode are presented graphically. © 1994.