Numerical calculations and Monte Carlo simulations have been performed to explore the dynamic scattering behavior of solutions of copolymer chains composed by two nonsegregated blocks in an athermal solvent. The experimental investigation of this property for real copolymers has shown the presence of a variety of modes, with a complex variation of their locations and intensities, which we try to partially understand with the present simulation work. For this purpose, we use chain samples with a small amount of polydispersity in the block lengths. We have introduced opposite scattering contrast factors for the A and B units. In this way, zero-averaged contrast conditions have been set for the overall scattering, avoiding the presence of the collective mode that otherwise would manifest for any nondilute system of polymer chains. Static and dynamic scattering functions have been obtained for the different samples. For the dynamic scattering functions, three modes are observed in agreement with existing experiments. The first two modes, predominant for low values of the scattering variable, q, are assigned to correspond to the compositional heterogeneity and the first internal relaxation of the chains, as predicted by the theory. The variation of the mode intensities and positions with polymer concentration is analyzed. It semiquantitatively agrees with the existing experimental data.