This paper describes some new finite difference methods of order 2 and 4 for computing eigenvalues of a two-point boundary value problem associated with a fourth
order differential equation of the form (py″)′′+(q−λr)y=0. Numerical results for two typical eigenvalue problems are tabulated to demonstrate practical usefulness of our methods.