En este trabajo se introduce un procedimiento iterativo para encontrar la solución de un modelo constitutivo no lineal para refuerzos de acero embebidos. El modelo presenta diferentes multiplicidades, donde los parámetros son seleccionados aleatoriamente dentro de una región de solubilidad. Para lograr esto, se modifica una clase de esquemas iterativos de punto fijo multipunto para raíces simples para encontrar raíces múltiples, logrando el cuarto orden de convergencia. Se emplean técnicas complejas de dinámica discreta para seleccionar los miembros con el rendimiento más estable. El problema mecánico mencionado anteriormente, así como algunos problemas académicos que involucran raíces múltiples, se resuelven numéricamente para verificar el análisis teórico, la robustez y la aplicabilidad del esquema propuesto.