wilk:= proc(n)
  product('(x-i)','i'=1..n);
end: