Abstract: A stepping scheme was described for analyzing the effective properties of composites containing multiple kinds of inclusions with different shapes and sizes. In the stepping scheme the inclusions are treated kind by kind and the effective stiffness coefficients of the resultant composites in each step are calculated. Numerical results indicate that the stepping scheme can precisely deal with the problem with high volume fraction of inclusions and multiple inclusions. For the binary composites consisting of matrix and one kind of inclusion, the present result is in agreement with that of the differential scheme.