factorial.m

n = 10;
fact = 1;           
                   
for k = 1:n         
   fact = k * fact;
   disp( [k fact] )
end