Un problema importante relacionado con el estudio de la estabilidad robusta de un sistema lineal que presenta variación en términos de un parámetro incierto consiste en comprender la variación en las raíces del polinomio característico del sistema en función del parámetro incierto. En esta contribución, proponemos un algoritmo para proporcionar condiciones suficientes sobre el parámetro incierto de tal manera que una familia robustamente estable de polinomios tenga todas sus raíces dentro de un subconjunto específico de su región de estabilidad. Nuestro método se basa en el teorema de Rouché y utiliza polinomios robustamente estables construidos mediante propiedades básicas de polinomios ortogonales.