python implementation of the sellmeier formula