function f=factorial(n);
f=1;
for i=1:n
    f=f*i;
end