Este trabajo presenta un método iterativo para resolver las ecuaciones de flujo de potencia a partir del establecimiento de funciones no lineales en cada nodo. Con el método propuesto se generan N-1 funciones no lineales, las cuales se resuelven iterativamente hasta lograr la convergencia en los voltajes de cada nodo del sistema eléctrico de potencia. El análisis de convergencia demostró que el método iterativo de dos pasos es de cuarto orden. En este trabajo se utiliza la fórmula original del método de Newton-Raphson en el primer paso y un método tipo Halley en el segundo paso. Los sistemas de prueba de 14, 30 y 57 nodos son simulados para conocer las características del método propuesto. Entre las fortalezas del método se puede mencionar que no requiere de una programación sofisticada y que representa una opción más para calcular flujos de potencia.