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: 1 >> 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+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. Compute the convergence order numerically (compare the error change when n doubled).