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