<p>Application of higher order finite element to the analysis of thin solid circular plates on elastic foundation using the modified Vlasov model in polar coordinate system is attempted. A Matlab code has been written for present formulation that can handle any type of loading and boundary conditions. Validation of the code is carried out after the convergence study. The results are compared to those obtained by other researchers. It is concluded that the present formulation behaves extremely well for thin solid circular plate resting on elastic foundations with good convergence rate and accuracy. The present plate bending element resting on elastic foundation using a higher order finite element is simple and requires less computational efforts, resources and time.</p>