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