function y = fact(N)

% File: FACT.M
% CALL: y = fact(N)
% This file computes N! (factorial).


y = 1;
while (N > 0)
  y = y*N;
  N = N-1;
end;
end;





