Use a better interpolation method in the function 'interstellar_extinction'

Instead of linear interpolation in wavelength, Akima cubic spline
interpolation in inverse wavelength is now used, which yields a
smoother extinction curve.