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