Write a MATLAB function named estimate. m that takes as input the number of points to use (call the variable npts) and returns the estimate for π using the Monte Carlo method described above. Next, write a script named distribution. m that calls estimate. m 10,000 times with the number of points used npts=100 and plots a histogram of the estimates for π (use the command histogram() for that). What do you notice about the distributions if you run distribution. m with npts=1000 or npts=10,000 compared to npts=100? (You do not have to report that in your homework; just change npts and see what happens). Whe