function out = newmapFeature(X, poly)
%NEWMAPFEATURE Summary of this function goes here
%   Detailed explanation goes here


out = ones(size(X(:,1)));
[m, n] = size(X);


out(:,end+1) = X(:,1);
out(:,end+1) = X(:,2);

%n is the number we need
for i = 1:n%1:n
    for j = 0:n

        if (i + j) > n
        continue
        end

        fprintf('X%f * X%f\n',i, (i+j));
        out(:, end+1) = X(:,i).*X(:,i+j);

        if i == n
        break
        end

        %if i == n break  if (i + j) > n break
    end
end

end