A review about the application of response surface methodology (RSM) in the optimization of analytical methods is presented. The theoretical principles of RSM and steps for its application are described to introduce readers to this multivariate statistical technique. Symmetrical experimental designs (three-level factorial, Box-Behnken, central composite, and Doehlert designs) are compared in terms of characteristics and efficiency. Furthermore, recent references of their uses in analytical chemistry are presented. Multiple response optimization applying desirability functions in RSM and the use of artificial neural networks for modeling are also discussed.