In this work, it is shown how the shapes of surface plasmon dispersion curves can be engineered by manipulating the distribution of the electromagnetic fields in multilayer structures, which themselves are controlled by the free electron density in metal-like materials, such as doped semiconductors in the THz spectral range. By having a nonuniform free electron density profile, reduced relative to that in typical bulk metals, the electromagnetic fields of surface plasmons are distributed in different metallic materials that have different complex dielectric permittivities. As the in-plane component of surface plasmon's wave-vector increases, they become more confined to a particular layer of the multilayer structure and have energies that are predictable by considering the permittivity of the layer in which the fields are most concentrated. Unusual and arbitrary shapes of surface plasmon dispersion curves can be designed, including stair steps and dovetails shapes.