Using a Monte Carlo numerical method, we perform a complementary analysis of the hyperfine field distribution in Magnetic/Non-Magnetic (M=Fe/NM=Cu) multilayers and a comparison with Mössbauer experiment results of [Fe(40)=Cu(20)]₂₀ multilayers. From differing relaxation rates for spins with few near neighbor spins and particularly those located at the interface, the Mössbauer result can be simulated. Our model consists of an alternate stacking of magnetic and non-magnetic layers (M_nM/NM_nNM) with disordered interface. The simulation results confirm that the concentration of interface alloys (M_1-xNM_x) and magnetic layer thickness modify systemically the magnetization distribution in the M/NM multilayers. The result is in agreement with Mössbauer experimental analysis.