The modeling of edge plasma of present and future magnetic confinement devices requires the knowledge of accurate line shapes. For instance, spectroscopic diagnostics using detailed line profiles may be used for obtaining the main plasma parameters. Detailed lines shapes also enter in the study of the radiation transfer in the divertor plasma by using edge plasma simulations [1,2]. For Lyman and Balmer series of hydrogen isotopes, the structure of several lines is affected by the Zeeman-Stark effect due to the simultaneous action of the plasma microfield, and of the externally applied magnetic field. In the future international fusion machine ITER [2], edge plasma conditions with a density around 5x10 14 cm-3 , and a temperature of the order or lower than 1 eV are expected. For such conditions, lower Lyman and Balmer lines may be affected by ion dynamics effects, i.e. a standard Stark broadening model with static ions may not be accurate enough. We have studied such cases with numerical integrations of the Schrödinger equation of emitters in a fixed magnetic field, and submitted to different microfield simulations [4,5,6]. Besides such benchmark simulations, a more computationally efficient modelling using the Fluctuation Frequency Model may provide fast calculations for systematic calculation of line shapes in the large range of conditions foreseen in the ITER divertor. For this wide density and temperature range, we will show comparisons between simulations and the Frequency Fluctuation Model, for Lyman α and Balmer α line.