La implementación de la relación de recurrencia o la ecuación de diferencia para generar polinomios ortogonales clásicos discretos es vulnerable a la propagación de errores. Este problema se aborda para el caso de las funciones de Krawtchouk, es decir, la base ortonormal derivada de los polinomios de Krawtchouk. Se propone un algoritmo para la determinación estable de estas funciones. Esto se logra definiendo puntos iniciales adecuados para el inicio de las recurrencias, equilibrando el orden de la dirección en la que se ejecutan las recurrencias y restringiendo de manera adaptativa el rango sobre el cual se aplican las ecuaciones. La adaptación está controlada por una desviación especificada por el usuario de la norma unitaria. Se proporciona el fundamento teórico, se explica el concepto algorítmico y se demuestra el efecto de la precisión controlada mediante ejemplos.