Question

I am trying to compute x1^i * x2^j * x3^k * ......

This is my code so far:

for l = 1:N
f = 1;
for i = 0:2
    for j = 0:2-i
        for k = 0:2-j
            for m = 0:2-k
                g(l,f) = x1(l)^i*x2(l)^j*x3(l)^k*x4(l)^m;
                f = f+1;
            end
        end
    end
end
end

How can I do this easier or without a loop?

Answer 1

I do not have MATLAB on hand here, but what I’d do is make a vector X = [x1, x2, …, xn] of bases and a vector P = [i, j, k, …, z] of powers, and then compute prod(power(X, P)).

power() does an element-wise power function, and prod takes the product of every element in the vector.