A very simple and computationally efficient numerical method is developed for the free vibration of arbitrary shape plates. A set of two- dimensional simple series functions is used as an admissible displacement functions in the Rayleigh-Ritz method to obtain the natural frequencies for the arbitrary shape plates. From the prescribed starting function satisfying only the geometric boundary conditions, the higher terms in the series function are constructed with adding order of polynomial. Natural frequencies are obtained for the arbitrary shape plates, with combinational boundary conditions. The obtained numerical results are presented, some cases are verified with other numerical methods in the literature.