In this paper, a novel analytical coupled trajectory optimization of a seven degrees-of-freedom (7DOF) Baxter manipulator utilizing extremum seeking (ES) approach is presented. The robotic manipulators are used in network-based industrial units, and even homes, by expending a significant lumped amount of energy, and therefore, optimal trajectories need to be generated to address efficiency issues. These robots are typically operated for thousands of cycles resulting in a considerable cost of operation. First, coupled dynamic equations are derived using the Lagrangian method and experimentally validated to examine the accuracy of the model. Then, global design sensitivity analysis is performed to investigate the effects of changes of optimization variables on the cost function leading to select the most effective ones. We examine a discrete-time multivariable gradient-based ES scheme enforcing operational time and torque saturation constraints in order to minimize the lumped amount of energy consumed in a path given; therefore, time-energy optimization would not be the immediate focus of this research effort. The results are compared with those of a global heuristic genetic algorithm (GA) to discuss the locality/globality of optimal solutions. Finally, the optimal trajectory is experimentally implemented to be thoroughly compared with the inefficient one. The results reveal that the proposed scheme yields the minimum energy consumption in addition to overcoming the robot's jerky motion observed in an inefficient path.