Warning: Polynomial regression via the normal equations tends to generate the badly conditioned Hilbert matrix and, thus, numerical problems. But can also do

Warning: Polynomial regression via the normal equations tends to generate the badly conditioned Hilbert matrix and, thus, numerical problems. But can also do the fitting using some orthogonal polynomials that are much more stable numerically.