Write a program in MATLAB to determine the value e accurate to $3$ digits after the decimal point. The main program should call a function that returns an estimate to the main program. Use an initial estimate of $1$ and $10000$ randomly generated points. The number of tries should be progressively increased until the desired accuracy is obtained. It should be noted that e is treated as an unknown, and comparison with the known value is not permitted. The program should have an upper limit on the number of iterations.