function poly=poly_fit(X,Y,p)
A=zeros(length(X),p+1);
for h=1:length(X)
    A(h,:)=X(h).^(0:p);
end;  
[r,c]=size(Y);
if r==1
   poly=(A'*A)\(A'*Y');
else
   poly=(A'*A)\(A'*Y);  
end;    
return;