function out=factorial(n);

out=1;
for i=1:n
  out=out*i;
end;  









