Abstract: Electromagnetic fields within a waveguide can be expressed in terms of the complex amplitudes of the electromagnetic modes it supports. The electromagnetic fields can be shaped by controlling the complex amplitudes of modes. Here, mode-converting metasurfaces are designed to transform a set of incident modes on one side to a different set of desired modes on the opposite side of the metasurface. A mode-converting metasurface comprises multiple inhomogeneous (spatially-varying) reactive electric sheets that are separated by dielectric spacers. The reactance profile of each electric sheet to perform the needed mode conversion is found through optimization. The optimization routine takes advantage of a multimodal solver that uses two main concepts: modal network theory and a discrete Fourier transform algorithm. With modal network theory, the modes can be translated between the electric sheets using matrix multiplication.