Побудовано однокрокову рекурентну схему інтегрування задачі Коші для великих систем звичайних диференціальних рівнянь, які виникають після просторової напівдискретизації початково-крайових задач для систем нелінійних параболічних рівнянь. Основною специфікою побудови цієї схеми є урівноваження порядків похибки частинами лінійної апроксимації за часом та похибки лінеаризації Ньютона.

Математичні методи та фізико-механічні поля [Текст] // Математичні методи та фізико-механічні поля - Львів : Інститут прикладних проблем механіки і математики. — Л. : Інститут прикладних проблем механіки і математики, 2001.

• УДК // Нелінійні рівняння та системи/Нелинейные уравнения и системы/Nonlinear equations and systems
• УДК // Поверхневі явища в контактних системах ( на поверхнях поділу)
• УДК // Числові методи розв'язання диференціальних рівнянь з частинними похідними