Due to schedule problems I am unable to perform your requested tests.
However, I already have an Excel file I had generated demonstrating the problem I have encountered. It is attached.
The call I make to get the results of Sheet2 is:
..
double[] fitVals = NI.CurveFit.PolynomialFit(tempsFixed, sensors, fitOrder, out coeffs, out err);
..
the "temps" column corresponds to the input array "tempsFixed", "sensor" corresponds to the input array "sensors", and fitOrder has the value 4. The column labelled "fit" displays the return values returned in the "fitVals" array. I extracted these data values by cutting and pasting from the Visual Studio Watch window during a breakpoint in program executions just following the statement above.
I tried scaling both the x and y data to each fit in the range -1.0 to +1.0 by subtracting the midrange value of each data set and then scaling each to the min-to-max range. Then the fit is correct as evidenced in the plot of Sheet1. I have derived transformations of the fit coefficients back to the original scaling and tested the results and am using this workaround strategy for now.