Numerical methods for eigenvalue problems in the description of drift instabilities in the plasma edge