In this work multivariate experiments were conducted to optimize the operating conditions for inductively coupled plasma optical emission spectrometry (ICP OES) for multielemental determinations in chocolate drink powder. The operating conditions were investigated using a 2(3) central composite design, where the variables studied were radio frequency power, nebulization flow rate, and auxiliary argon flow rate. The effects of these parameters on plasma robustness and on signal to background ratio (SBR) were considered in parallel, allowing the evaluation of robustness and detectability using few and fast experiments to select the best conditions for the determination of the analytes. In this case, the proposed experiments were applied to the optimization of a method aimed at the determination of Al, Ba, Cd, Co, Cr, Cu, Fe, Mg, Mn, Mo, Ni, P, Pb, V, and Zn in chocolate drink powder. The compromise conditions that allowed obtaining a robust and sensitive analytical method were radio frequency power of 1200 W, nebulization flow rate of 0.6 L/min, and auxiliary argon flow rate of 0.3 L/min. Using these conditions, recoveries between 95 and 105% and relative standard deviations lower than 5% were obtained for the majority of the analytes. The proposed method was successfully applied to the analysis of 15 samples of chocolate drink powder. The highest concentrations of metallic species were found in diet and light products.