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 (fun, a, b, n) where fun is a function handle of f(x), and a, b is the interval, and n is the number of sub-intervals (i.e., 2n+l 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. Compute the convergence order numerically (compare the error change when n doubled). What to hand in? Print out your function in the file Simpson's, your script file, the plot of error, and your comments