function f=my_factorial(n)
%
f=1;
for k=1:n
    f=f*k;
end
