factorial

// Make a table of factorial
n = (0:30)';
[n factorial(n)]

// See the limits of factorial: f(171)=%inf
factorial(170) // 7.257415615307998967e306
factorial(171) // %inf

// Plot the function on all its range.
scf();
plot ( 1:170 , factorial , "b-o" )
h = gcf();
h.children.log_flags="nln";