Factor out common spline functions.
Add plot method for ksfit objects.
Add functions for visually evaluating fits.
Add working baseline code for gradient based fitting.
Check in documentation.