In this work, the model of adiabatic waveguide modes is studied by means of computer algebra. Within the model, the solution of the system of Maxwell’s equations is reduced to a form expressed via the solution of a system of four ordinary differential equations and two algebraic equations for six components of the electromagnetic field. In the case of multilayer waveguides, by means of a computer algebra system, the equations are reduced to a homogeneous system of linear algebraic equations, which is studied symbolically. The condition for non-trivial solvability of the system defines a dispersion relation, which is solved by the symbolic-numerical method, while the system is solved symbolically. The paper presents solutions that describe adiabatic waveguide modes in the zeroth approximation, taking into account the small slope of the interface of the waveguide layer, which are qualitatively different from solutions that do not take into account this slope.