function f=fact(n)
% f=fact(n) implementa il calcolo del fattoriale di n
% f = 1 * 2 * ... * n

f = 1;
for k = 1 : n % k = n: -1: 1
    f = f*k;
end