wolffd@0: function x = factorial(n)
wolffd@0: % FACTORIAL Compute n!
wolffd@0: % x = factorial(n)
wolffd@0: 
wolffd@0: if n == 0
wolffd@0:   x = 1;
wolffd@0: else
wolffd@0:   x = n*factorial(n-1);
wolffd@0: end