Fitting polynomial passing through a fixed point

This function will fit a polynomial passing though a fixed point using least square method

Syntax:

P = polyfix(xi,yi,x0,y0,m)

This function fits a polynomial passing thru a fixed point of (x0, y0).

(xi, yi) are the data points to be fited.

It will give a result vector of P = [P1,P2,P3,...,Pm,Pm+1], corresponding to the polynomial of:

y = P1 x^m + P2 x^(m-1) + ,..., + Pm x + Pm+1

Please note that xi, and yi must be one dimensional vectors. Multidimensional data is not yet supported in this version.