Dynamic programming with successive approximation has been used in the past for optimizing multi-reservoir water resources systems. In this study, the State Incremental Dynamic Programming (SIDP) model is developed for energy optimization of multi-reservoir systems. A random file access method is used for reaching initial and intermediate data to cope with the curse of dimensionality of dynamic programming. A conventional dynamic programming method is used for each single reservoir to find the initial trajectory of the reservoirs. Then, the computer program developed in the study is applied to the multipurpose-multi-reservoir system in Lower Seyhan Basin, which has six reservoirs, some of which are serial and some parallel. First, extended historical flows were used to maximize firm energy in the critical period, and then total energy in the total flows. The program was run with 50-year long segments (20 flow scenarios) of the synthetic flow data generated by using the HEC-4 generalized computer program to take into account the stochastic nature of streamflows. An increment of approximately 20 percent in total energy was obtained by using the model for the Lower Seyhan System, as compared to that calculated previously by conventional methods..