bc_poly function

Fit and subtract a polynomial to each spectrum in a dataset.