En este trabajo especial de grado se implementó el algoritmo del backward/ forward sweep (BFS) como método para resolver redes radiales o débilmente malladas, como es normalmente la estructura de los sistemas de distribución. Dicho método utiliza las características típicas de estas topologías de redes para simplificar las ecuaciones del sistema y así optimizar el tiempo consumido para encontrar la solución de la red. Sin embargo, este tipo de métodos presenta particulares problemas con los límites de convergencia que han sido estudiados para el análisis de sistemas eléctricos sin presencia de componentes armónicas. En este trabajo se amplió el estudio y se identificaron los límites de convergencia del método en presencia de componentes armónicas, mediante la implementación de dos variantes del BFS aplicadas a una red trifásica desbalanceada en presencia de cargas no lineales. La implementación se realizó mediante la programación en Matlab® de los algoritmos estudiados y su aplicación a tres redes de prueba, con la posterior verificación de los resultados obtenidos.