Este trabalho descreve um estudo teórico e sua aplicação computacional sobre a operação de usinas hidráulicas num sistema elétrico de potência. Parte-se das equações dinâmicas do sistema, considerando os parâmetros que influenciam um sistema multireservatório (afluências e potência demandada, aleatória, tempo de percurso da água entre dois reservatórios). Essas equações são um conjunto de equações diferenciais de 1º ordem, quais são aproximadas por um conjunto de diferentes finitas, com a finalidade de que o sistema em um instante t só dependa do instante anterior (t-1). Quando apresenta-se o problema do tempo de percurso da água entre dois reservatórios, t(i), o sistema em um instante t vai depender de (t-l) e de (t-t(i)). Então as equações recorrentes com transmissão não instantânea são transformadas em um conjunto de equações de estado, que representam estados intermediários entre os reservatórios, assim aumentando o número de equações de estado, faz-se com que o retardo de tempo t(i) não apareça explicitamente. Tendo definido o modelo matemático do sistema, o objetivo é maximizar a quantidade de água armazenada, satisfazendo uma demanda de potência aleatória e as restrições em hidráulicidade (características hidráulicas do sistema) e em instalações. São usados os métodos de relaxação e de intervalos iguais. Os programas foram feitos para o sistema IBM 370-165 da PUC/RJ.