Este trabalho apresenta o resultado obtido da implementação de um algoritmo de cálculo do fluxo de potência em redes de distribuição radiais trifásicas desequilibradas levando em consideração o condutor neutro e o efeito da terra. Esta implementação é parte do projeto de desenvolvimento de uma ferramenta para ajuste e dimensionamento de proteções em redes de distribuição. O algoritmo iterativo foi implementado no Matlab a partir de uma adaptação do método do somatório de correntes. O algoritmo foi testado com os dados de uma rede real de baixa tensão de 29 barras. Os valores de tensão calculados foram comparados com os valores obtidos pelo software OpenDSS e o resultado indica que o algoritmo funciona satisfatoriamente.