# subfactorial(n) is n!;
# subfactorial(x_):=gamma(x+1,-1)/e;


subfactorial(0):= 1;
subfactorial(1):= 0;
subfactorial(2):= 1;
subfactorial(3):= 2;
subfactorial(4):= 9;
subfactorial(5):= 44;