This paper presents an algorithm for constructing mixture designs based on orthogonal arrays of index unity containing the least number of runs for a specified strength and a specified number of levels using pairwise orthogonal Latin squares. The algorithm permits the generation of economical and efficient mixture designs for fitting Scheffe's canonical polynomials.