The differential equations describing stellar oscillations are transformed into an algebraic eigenvalue problem. Frequencies of adiabatic oscillations are obtained as the eigenvalues of a banded real symmetric matrix. We employ the Cowling-approximation, i.e. neglect the Eulerian perturbation of the gravitational potential, and, in order to preserve selfadjointness, require that the Eulerian pressure perturbation vanishes at the outer boundary. For a solar model, comparison of first results with results obtained from a Henyey method shows that the matrix method is convenient, accurate, and fast.