modify fit_function

[jmzobitz]

fit function only does linear interpolation. We want to modify that so it can do:

• spline interpolation w/ more than 3 data points.
• linear interpolation with 2 or 3 data points.