This work connected with a micro stress and micro deformation determination problem in matrix composites with random inclusions. Stochastic structure of the composites is described with set of moment functions, with help of which deformation fields statistical characteristics are searching. The boundary problem is reduced to integro-differential equation in pulsations of displacements. Solution of the problem is searched in second approach. Analytical expressions for stress and deformation fields statistical characteristic are obtained. Numerical results for average quantities for simple shear macroscopic mode of deformation are presented.