discover function factorial(n). change factorial.m to test_factorial_limit.m