estart; xh := proc (x)
local i,cj;
cj := n;
for i to x do
cj := cj*(n-i)
end do;
print(expand(cj))
end proc;
M := [seq(a,a = 1 ..4)];
for j in M do xh(j) end do;
2
n - n
3 2
n - 3 n + 2 n
4 3 2
n - 6 n + 11 n - 6 n
5 4 3 2
n - 10 n + 35 n - 50 n + 24 n