This article deals with the existence of multiple positive solutions to the following system of nonlinear equations involving Pucci’s extremal operators: (Formula presented.) where Ω is a smooth and bounded domain in RN and fi:[0,∞)×[0,∞)⋯×[0,∞)→[0,∞) are Cα functions for i=1,2,⋯,n. The multiplicity result in this work is motivated by the work Amann (SIAM Rev 18(4):620–709, 1976), and Shivaji (Nonlinear analysis and applications (Arlington, Tex., 1986), Dekker, New York, 1987), where the three solutions theorem (multiplicity) has been proved for linear equations. Later on, it was extended for a system of equations involving the Laplace operator by Shivaji and Ali (Differ Integr Equ 19(6):669–680, 2006). Thus, the results here can be considered as a nonlinear analog of the results mentioned above. We also have applied the above results to show the existence of three positive solutions to a system of nonlinear elliptic equations having combined sublinear growth by explicitly constructing two ordered pairs of sub and supersolutions.