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.