We develop theory and provide an algorithm that allows us to find eigenpairs of a class of (unbounded) self adjoint operators that are localised in a user specified subdomain. We show a practical realisation of the algorithm through numerical experiments.