Robust algorithm for computing quasi-static stark broadening of spectral lines

A method previously developed to solve large-scale linear systems is applied to quasi-static Stark line broadening calculations. The method is formally exact, numerically stable, and allows optimization of the integration over the quasi-static field to assure numerical accuracy. Furthermore, the method is computationally efficient compared to the conventional approach, particularly for complex problems involving large matrices.