Simpson's Rule in Matlab Write a Matlab function which computes the integral by Simpson's rule. Your function should be used by the following command in Matlab command window: >> v = Simpson ('funItg' a, b, n) where funItg. m is the name of the file of the function f(x), and a, b is the interval, and n is the number of sub-intervals (ie., 2n+1 will be the number of points.) Test your function on problem 1. Write a script that computes, for n = 2, 4, 8, 16, 32, 128, the absolute error for each n, and a plot (with loglog) of the absolute error against n. How does the error change when n doubled? Compare the results with trapezoid rule and comment. What to hand in? Your function in Simpson.m, your script file, the plot of error, and your comments