The elastic interaction energy between several precipitates is of interest since it may induce ordering of precipitates in many metallurgical systems. Most of the works on this subject assumed homogeneous systems, namely, the elastic constants of the matrix and the precipitates are identical. In this study, the elastic fields, and self and interaction energies of inhomogeneous anisotropic precipitates have been solved and assessed, based on a new iterative approach using the quasi-analytic Fourier transform method. This approach allows good approximation for problems of several inhomogeneous precipitates in solid matrix. We illustrate the calculation approach on γ′-Ni3Ti precipitates in A-286 steel and demonstrate that the influence of elastic inhomogeneity is in some incidences only quantitative, while in others it has essential effect. Assuming homogeneous system, disk shape precipitate is associated with minimum elastic energy. Only by taking into account different elastic constants of the precipitate, the minimum self-energy is found to be associated with spherical shape, and indeed, this is the observed shape of the precipitates in A-286 steel. The elastic interaction energy between two precipitates was calculated for several configurations. Significant differences between the interactions in homogeneous and inhomogeneous were found for disk shape morphologies. Only quantitative differences (9% higher interaction between inhomogeneous precipitates) were found between two spherical precipitates, which are the actual shape.