function f=fact(n)
% f=FACT(n) calcola il fattoriale di n
f=1;
for k=1:n
    f=f*k
end
end