A consistent approach to estimating nuclear effect functions RA RvA (x2) and RSA(x2) based on numerical iteration technique is presented in the quark-parton model when taking into account the nonconstancy of quantum chromodynamics correction factor K. ARv (x2) and RsA(x2) correspond respectively to the valence quark distributions for one bound nucleon within the nucleus and to the sea quark ones. Related numerical analysis is given for nuclei 6 12C,20 40Ca, and 26 56Fe. As the basis, it adopts both experimental data of the high energy proton-nucleus Drell-Yan process and of the high energy lepton-nucleus deep inelastic scattering.