function y= f(n)  % f.m factorial function
 if (n==1)
   y=1;
 else
   y=n*f(n-1);
 end
end
