Codes for equivalent inhomogeneous inclusion method

This code determines the stress and strain fields of multiple spherical /circular inhomogeneities by Equivalent Inhomogeneous Inclusion Method (EIIM). EIIM is an inverse of Equivalent Inclusion Method and simulates the subdomains of matrix with known strains with equivalent inhomogeneous inclusions. The stress and strain fields of multiple inhomogeneities are decomposed into the superposition of matrix under applied load and each solitary inhomogeneous inclusion with polynomial eigenstrains by the iteration of EIIM.