A set of two-dimensional (2-D), second-order approximate equations for extensional, thickness-stretch and symmetric thickness-shear vibrations of piezoelectric ceramic plates with electroded faces is extracted from the infinite system of 2-D equations deduced previously. The new truncation procedure developed recently is used for it improves the accuracy of calculated dispersion curves. Closed-form solutions are obtained for free vibrations of circular disks of barium titanate. Dispersion curves calculated from the present approximate 2-D equations are compared with those obtained from the 3-D equations, and the predicted resonance frequencies are compared with experimental data. Both comparisons show good agreement without any corrections. The frequencies of the edge modes calculated from the present 2-D equations are very close to the experimental data. Furthermore, mode shapes at various frequencies are calculated in order to identify the frequency segments of the spectrum at which one of the coupled modes--i.e., the radial extension (R), edge mode (Eg), thickness-stretch (TSt), and symmetric thickness-shear (s.TSh)--is predominant.