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