An iterative method for computing robustness of polynomial stability