An iteractive Kantorovich method is presented for the vibration analysis of rectangular orthotropic thick plates. Mindlin plate characteristic functions are derived in general forms using the Kantorovich method. Initially, Timoshenko beam functions consistent with the boundary conditions of the plate were used. Through numerical calculations of natural fairs of appropriate models, it has been confirmed that the method presented is superior to the Rayleigh-Ritz analysis or the finite element analysis in both accuracy and computational efficiency.