Mars Wavelength Correction
Continuum removal with natural spline

Gives results that match idl, but is terribly slow because it work sin
decimal. Need to convert the decimal spline fitting to work with regular
arrays.