|MATLAB Function Reference|
is the product of all the integers from
prod(1:n). Since double pricision numbers only have about 15 digits, the answer is only accurate for
n <= 21. For larger
n, the answer will have the right magnitute, and is accurate for the first 15 digits.