The existence of multiple positive periodic solutions of the system of difference equations with a parameter x(n + 1) = A(n, x(n))x(n) + λf(n, xn), is studied. In particular, we use the eigenvalue problems of completely continuous operators to obtain our results. We apply our results to a well-known model in population dynamics.