function A = build_A(X,Y,n)
    A = []
    for i = 0:1:n
        for j = 0:1:(n-i)
            A = [A,(X.^i) .* (Y.^j)]
        end
    end

end