Numerical methods based on shifted polynomials in swirling flows stability analysis