The purpose of this study is to investigate both the fundamental and some higher natural frequencies and mode shapes of compressive members resting on the linear elastic foundation. The model of compressive member is based on the classical Bernoulli-Euler beam theory. The differential equation governing free vibrations of such members subjected to an axial load is derived and solved numerically for calculating the natural frequencies and mode shapes. The Improved Euler method is used to integrate the differential equation and the Determinant Search method combined with the Regula-Falsi method to determine the natural frequencies, respectively. In numerical examples, the hinged-hinged, hinged-clamped, clamped-hinged and clamped-clamped end constraints are considered. The convergence analysis is conducted for determining the available step size in the Improved Euler method. The validation of theories developed herein is also conducted by comparing the numerical results between this study and SAP 90. The non-dimensional frequency parameters are presented as the non-dimensional system parameters: section ratio, modulus parameter and load parameter. Also typical mode shapes are presented.