function k=my_fact(n)
if nargin~=1, error('輸入變量個(gè)數(shù)錯(cuò)誤，只能有一個(gè)輸入變量'); end
if nargout>1, error('輸出變量個(gè)數(shù)過多'); end
if abs(n-floor(n))>eps | n<0 % 判定 n 是否為整數(shù)
   error('n 應(yīng)該為非負(fù)整數(shù)');
end   
if n>1      % 如果 n>1, 進(jìn)行遞歸調(diào)用
   k=n*my_fact(n-1);
elseif any([0 1]==n) % 0!=1!=1
   k=1;
end