This paper presents a numerical solution for vibration analysis of cantilevered non-uniform trapezoidal thick plates. Based on the first shear deformation theory, kinetic and strain energies of the plate are derived and using Hamilton's principle, governing equations and boundary conditions are derived. A transformation of coordinates is used to convert the equations and boundary conditions from the original coordinates into a new computational coordinates. Using Differential quadrature method (DQM), natural frequencies and corresponding modes are derived numerically. Convergence and accuracy of the proposed solution are confirmed using results presented by other authors and also results obtained based on the finite element method using ANSYS software. Finally, as the case studies, two cases for variation of thickness are considered and the effects of angles, aspect ratio and thickness of the plate on the natural frequencies are studied. It is concluded that two angles of the trapezoid have opposite effect on the natural frequencies. Also, it is shown that all frequencies rise as value of thickness increases or value of the aspect ratio of the plate decreases. The most advantage of the proposed solution is its applicability for plates with variable thickness.