function m=facto(n)
% m=facto(n), factorielle
if n==0,m=1;
elseif n==1,m=1;
elseif n==2,m=2;
elseif (floor(n)-abs(n))==0,
m=n*facto(n-1);
else error('n doit etre un entier positif');
end;