Предложен итерационный метод определения напряженно-деформированного состояния в линейно упругих механически неоднородных системах, содержащих произвольно расположенные упругие включения. Эффективность предложенного метода проиллюстрирована на решении задач, в которых используются разные математические аппараты (теория функций комплексного переменного, представление Папковича-Нейбера, граничные интегральные уравнения).