В течениях вязкой несжимаемой жидкости имеются пограничные слои с быстрым изменением компонент решения, для успешного численного исследования которых представляется полезным повышать порядок точности разностных схем или использовать неравномерную сетку, сгущающуюся в зонах больших градиентов. В дипломной работе совмещаются оба этих подхода в комплексе, т. е. применяется компактная схема второго порядка точности на неравномерной сетке. Разработан итерационный алгоритм реализации схемы с граничными условиями того же порядка точности и программа расчета краевых задач в операционной среде МАТЛАБ. Проведены расчеты тестовых задач в сравнении с другими методами.