An experimental design and response surface methodologies using Plackett–Burman and Box-Behnken designs were applied for selecting and optimizing the most appropriate parameters which significantly affect the separation and quantitative estimation of five skeletal muscle relaxants and four analgesic drugs (baclofen, methocarbamol, dantrolene sodium, orphenadrine citrate, cyclobenzaprine hydrochloride, ketoprofen, etoricoxib, ibuprofen, and mefenamic acid) with a relatively short duration of analysis in a single run. For the separation of the nine drugs, an INERTSIL ODS-V3-5 µm C18 column (250 × 4.6 mm I.D.) was used with the optimum mobile phase conditions (45.15 mM ammonium acetate buffer pH 5.56 adjusted with acetic acid, acetonitrile, and methanol in a ratio of 30.5:29.5:40, v/v/v with a flow rate of 1.5 mL/min) and UV-detection at 220 nm. The optimized method was successfully subjected to the validation steps as described in ICH guidelines for linearity, precision, accuracy, robustness, and sensitivity. The optimized and validated method was effectively applied to determine the content of the studied drugs in their pharmaceutical preparations and to expand its applicability to the counterfeit estimation of etoricoxib in different brands of tablet dosage forms.