The propagator method (PM) belongs to a class of subspace based methods for direction-of-arrival estimation which only requires linear operations but does not involve any eigendecomposition or singular value decomposition as in common subspace techniques. In this paper, we apply the PM for estimating the frequencies of multiple real sinusoids in noise and a computationally simple as well as high resolution multiple frequency estimation algorithm is developed. The estimation accuracy of the proposed method is contrasted with the conventional MUSIC and Cramer-Rao lower bound under different noise conditions.