A general framework is developed for constructing higher order spectral ref
inement schemes for a simple eigenvalue. Well-known techniques for ordinary
spectral refinement are carried over to higher order spectral refinement y
ielding faster rate's of convergence. Numerical examples are given by consi
dering an integral operator.