A modelagem dos mercados de energia modernos exige a representação das seguintes características principais: (i) um processo de decisão dinâmico estocástico , com incertezas relacionadas aos os custos de produção e dos combustíveis renováveis, entre outros; e (ii) teoria dos jogos que representa o comportamento estratégico de múltiplos agentes , por exemplo, em propostas de preços diárias. Esses recursos podem ser, em teoria, representados como uma recursão de programação dinâmica estocástica, onde temos um equilíbrio de Nash para múltiplos agentes. No entanto, o problema resultante é muito difícil de resolver. Este trabalho apresenta um processo iterativo para resolver o problema acima para sistemas de energia realistas. O algoritmo proposto é composto de um algoritmo de ponto fixo, no qual, cada passo é resolvido através do método de programação dinâmica dual estocástica. A aplicação do algoritmo proposto está ilustrado em estudos de caso com os sistemas de energia Panamá.