A class of problems of wave propagation in waveguides consisting of one or several layers that are characterized by linear variation of the squared refractive index along the normal to the interfaces between them is considered in this paper. In various problems arising in practical applications, it is necessary to efficiently solve the dispersion relations for such waveguides in order to compute horizontal wavenumbers for different frequencies. Such relations are transcendental equations written in terms of Airy functions, and their numerical solutions may require significant computational effort. A procedure that allows one to reduce a dispersion relation to an ordinary differential equation written in terms of elementary functions exclusively is proposed. The proposed technique is illustrated on two cases of waveguides with both compact and non-compact cross-sections. The developed reduction method can be used in applications such as geoacoustic inversion.