Abstract: The theory and methodology of a four-level orthogonal array design were developed for the optimization of analytical procedures The construction and characteristics of an OA16(45) matrix are described in detail  and the orthogonality of the matrix is shown by means of a third-order regression model