% Input: n, vettore
% Output: m, vettore dei n!!
function m = double_factorial(n)
  m = 0*n ;
  for i = 1 : length(n)
    m(i) = prod( n(i):-2:1 ) ;
  end
end