factorial := proc ( n )
  if (n=0) then
    return(1);
  else
    return(n * factorial(n-1));
  fi;
end: