We propose simple Monte Carlo method for calculating parton distribution in nuclei. Only events satisfying the exact kinematical constrains of the corresponding deep-inelastic reaction probing given nuclear distribution are selected to form the final distribution we are looking for. The EMC effect is automatically included by means of two parameters, which characterize the change of the nuclear pion field. Good agreement with experimental data in the broad range of variable x is obtained.

PACS numbers: 12.38.Aw, 12.38.Lg, 12.39.--x