Numerical methods for a singular eigenvalue problem with eigenparameter in the boundary conditions