Periodic solutions for systems of coupled nonlinear Schrödinger equations with three and four components

Periodic solutions for systems of coupled nonlinear Schr\"odinger equations (CNLS) are established by the Hirota bilinear method and elliptic functions. The interesting feature is the choice of theta functions in the formulation. The sum of moduli of the components or the total intensity of the beam in physical terms, will now be a rational function, instead of a polynomial, of elliptic functions. Each component of the CNLS may have multiple peaks within one period.