PolyFit functions - should it return up to order + 1 actual values?

Question

I have been using the PolyFit function in CVI 2010 (version 10.0.0 (360)) and I am trying to figure out if I am implementing incorrectly, there is a bug for higher order or if I am getting the exepected behavior.
The function always return zero for coefficients above 26 (order equals between 45 to 90 ). I could understand this if the fit were always good but in cases where the fit is obviously less than optimum the first twenty seven coefficients (o through 26) are populated but the remainder never are. Is this an expected behavior?
The application is the removal of baseline roll / drift from mass spectrometry data. The file sizes are large (4 MPoint and up) which must be tackled in smaller chunks (at least currently). One may can correct misfit data by limiting the chunk size manaually but the end goal is routine,  robust automation across multiple blind spectra.

Answer 1

In our case the baseline roll is relatively slow compared with the real signals. However, the correction using the current polyfit has an effective  limit of 500 to 1500 Hz (chi by eye). This requires several regions to be fit and then the  "joints" to be smoothed to cover the approxiamtely 1 MHz bandpass that we are looking for. Such  a piecewise fit is doable but it requires a lot more knowledge of the system if the analysis is to be performed in batch mode (no human intervention.) 
The main question of what should the return values be suggest that the function should be upward bounded to a value less than 90 (the current limit as described in the  documentation and reported by the debugger ).