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