Free vibration analysis using the method of NDIF (non-dimensional dynamic influence function), which was developed by the author, is extended to arbitrarily shaped polygonal plates with free edges. Local Cartesian coordinate systems are employed to apply the free boundary condition to nodes distributed along the edges of the plate of interest. Furthermore, a new way for applying the free boundary condition to nodes located at corners of the plate is for the first time introduced by considering the phenomenon of stress concentration at the corners. Two case studies show that the proposed method is valid and accurate when the eigenvalues by the proposed method are compared to those by FEM(ANSYS).