sub factorial
{
  my $n = @_;