function E=Efac(x)

N=numel(x);
E=N;
for n=1:N
  E(n)=prod(x(n):-1:0);
end
end