Vectorizing the creation of a matrix of successive powers

Question

Let x=1:100 and N=1:10. I would like to create a matrix x^N so that the i^th column contains the entries

Answer 1

Not sure if it really fits your question.

bsxfun(@power, cumsum(ones(100,10),2), cumsum(ones(100,10),1))

EDIT: As pointed out by Adrien, my first attempt was not compliant with the OP question.

xn = 100;
N=10;
solution = [ones(1,xn); bsxfun(@power, cumsum(ones(N,xn),2), cumsum(ones(N,xn),1))];