Rewrite Voigt profile calculation

Adapted from the Armstrong 1967 code