Mosaic crystals of b-type, for example, pyrolytic graphite, is widely used in experimental physics as neutron and X-ray monochromators and provide a higher yield of monochromatic radiation than perfect ones. The technique is proposed and implemented to calculate the reflectivity of such crystals by Monte Carlo modeling, correctly considering the multiple reflections of photons inside the crystal and the geometry of experiment for random distribution of the mosaic
