The simulation of polymer-dispersed cholesteric liquid crystals is carried out using a model, which takes into account the domain structure of the layer. The model is based on the 4x4 matrix calculation method of Berreman. One considers different orientations of the helical axis in multi- and single-domain configurations. The distribution of the helical axis orientations in the multi- domain model allows the description of optical properties of multi-stable director configurations.