На основе адиабатического представления собственных мод интегрально-оптического многослойного волновода для вертикального распределения электромагнитного поля в волноводе получены обыкновенные дифференциальные уравнения и граничные условия. Для плавно-нерегулярных волноводов применен асимптотический метод и выделены вклады нулевого порядка малости в дифференциальных уравнениях и граничных условиях. Получены явные выражения для вертикального распределения электромагнитного поля в волноводе и граничных условий на границах слоев. В итоге задача сведена к решению однородной системы линейных алгебраических уравнений, зависящей от спектрального параметра, и поиску значений параметра, при которых система нетривиально разрешима. В заключение рассмотрены метод и алгоритмы вычисления вертикального распределения адиабатических мод.