The aim of the present study was to optimize a chromatographic method for the analysis of atorvastatin (acid and lactone forms), ortho- and para-hydroxyatorvastatin by using an experimental design approach. Optimization experiments were conducted through a process of screening and optimization. The purpose of a screening design is to identify the factors that have significant effects on the selected chromatographic responses, and for this purpose a full 2³ factorial design was used. The location of the true optimum was established by applying Derringer’s desirability function, which provides simultaneously optimization of all seven responses. The ranges of the independent variables used for the optimization were content of acetonitrile in mobile phase (60–70%), temperature of column (30–40 °C) and flow rate (0.8–1.2 mL min⁻¹). The influences of these independent variables were evaluated for the output responses: retention time of first peak (p-hydroxyatorvastatin) and of last peak (atorvastatin, lactone form), symmetries of all four peaks and relative retention time of p-hydroxyatorvastatin. The primary goal of this investigation was establishing a new simple and sensitive method that could be used in analysis of biological samples. The method was validated and successfully applied for determination of atorvastatin (acid and lactone forms) and its metabolites in plasma.