function a = poly_fit( x, y, n )

[~, m] = size(x);
A = repmat(x', [1, n+1]);
[X,~] = meshgrid(0:n, 0:m);
size(A)
size(X)
%A = A .^ X

end

