In this paper, a non-linear programming model to operate distribution systems considering energy storage devices and distributed generation is presented. Mathematical formulation is made taking in to account four terms for minimization: The cost operation of the electrical grid, reducing greenhouse emissions, reduction of electrical losses in conductors and voltage profile improve. As model constraints are employed the active a reactive nodal power balance, maximum capacities from distributed generation and energy storage devices, and the voltage profile regulation. To solve the proposed mathematical model a commercial optimization software GAMS and CONOPT solver are employed. To verify the efficiency and applicability of the model developed are used two test nodes 9 and 10, which have been adapted to the typical operating conditions of Colombian distribution systems. The results show the ability to adapt to various conditions and operating instructions of the utility.