Modeling the Stark broadening of spectral lines in plasmas is a complex problem. The problem has a long history, since it plays a crucial role in the interpretation of the observed spectral lines in laboratories and astrophysical plasmas. One difficulty is the characterization of the emitter’s environment. Although several models have been proposed over the years, there have been no systematic studies of the results, until now. Here, calculations from stochastic models and numerical simulations are compared for the Atoms 2014, 2 300 Lyman-α and -β lines in neutral hydrogen. Also discussed are results from the Helium-α and -β lines of Ar XVII.