Computes the least squares fitting for a polynomial of order k, for a data-set in a file with columns for x and y.
The condition k< n should be satisfied.