A new numerical method—Green quasifunction is proposed.The idea of Green quasifunction method is clarified in detail by considering a vibration problem of simply-supported thin polygonic plates on Pasternak foundation.A Green quasifunction is established by using the fundamental solution and boundary equation of the problem. This function satisfies the homogeneous boundary condition of the problem.The mode shape differential equation of the vibration problem of simply-supported thin plates on Pasternak foundation is reduced to two simultaneous Fredholm integral equations of the second kind by Green formula.There are multiple choices for the normalized boundary equation.Based on a chosen normalized boundary equation,a new normalized boundary equation can be established such that the irregularity of the kernel of integral equations is overcome.Finally,natural frequency is obtained by the condition that there exists a nontrivial solution in the numerically discrete algebraic equations derived from the integral equations.Numerical results show high accuracy of the Green quasifunction method. A new numerical method-Green quasifunction is proposed. The idea of ​​Green quasifunction method is clarified in detail by considering a vibration problem of simply-supported thin polygonic plates on Pasternak foundation. A Green quasifunction is established by using the fundamental solution and boundary equation of the problem. This function satisfies the homogeneous boundary condition of the problem. The mode shape differential equation of the vibration problem of simply-supported thin plates on Pasternak foundation is reduced to two simultaneous Fredholm integral equations of the second kind by Green formula. there are multiple choices for the normalized boundary equation, a new normalized boundary equation can be established such that the irregularity of the kernel of integral equations is overcome. Finally, natural frequency is obtained by the condition that there exists a nontrivial solution in the numerically discrete algebraic equations d erived from the integral equations. Numerical results show high accuracy of the Green quasifunction method.