Evaluation of the polynomials should be done by normalizing the data
as indicated above.
Something like:
const real8 const real8 minL = master.originalwindow.linelo; const real8 maxL = master.originalwindow.linehi; const real8 minL = master.originalwindow.pixlo; const real8 maxL = master.originalwindow.pixhi; matrixreal4(N,1) l\_axis = linenumbers; normalize(l\_axis); matrixreal4 f = polyval(l\_axis, p\_axis, coefficients, [degree]);
NOTE: It is faster to evaluate a polynomial on a grid than point by point.